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PAPERS 
FLOW INA CHANNEL OF DEFINITE 


“ 


he effect of anliatit roughness of the wall upon the resistance to 7 
es has been thoroughly : studied, there seems no equally — 
satisfactory study of the effect of known wall roughness on the flow in open 
channels. . This paper Teports t the ‘results of some e experiments bearing on this - 


- question w which were performed by ‘the writer during the summer of 1940 in the © 
laboratory. of the Institute of Hydraulic Research at Iowa City, Iowa. AD 
“total c of two hundred runs on a » rectangular flumé, with e eleven different rough- 


“nesses an and four different slopes, indicated that Chézy’ could be 

as 40 logio R/e + cs, in which R is the hydraulic ra radius, eisa measure of the 

“Toughness, and ¢, is a small shape correction. Ss ‘The mean of the deviations of 


C=4 40 R/C. + Ge +4 
The mean deviation of this Riedie from 


value of C. Shooting flow in the same gave the formula 
¢ 60 R/C + 30 F + — 175, in which F is the Froude number. ‘The 


mean deviotion of formula from the data was. 2. 44, or about 2.1% of the 


iid these pon are cosieat, ‘two ‘surprising results follow: (1) The ‘ universal 

constant,” "b= 2. 5, changes 2. 5 for pipes to ut 3.0 in open channels; 

and 1 (2) the law of ‘resistance to shooting flow i is qu 


‘of resistance | to tranquil 


_ Nors.—Written comments are invited for r immediate publication; to to insure publication the > last dis- 
= should be submitted by May 11,1945. 


Prof. of Mechanics, Ohio State te Univ., Columbus, Ohio. 
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December, CHANNEL ROUGHNE SS 


The pioneer sr research on flow in 1 open 
artificial | roughness was made by H. Bazi n?3 between (1855 and 1860. These 
tests w ere made in a channel 2 m wide and nearly | 600 m long, and the measure= is 
“ments were quite accurate; but only a few of the  Toughnesses used were SO 
definitely described that ‘they - could be accurately duplicated today, and appa-_ 
“rently Bazin made no attempt to reduce his results to a formula that included — a 


= definite > roughness figure. In fact, the idea of a definite linear unit of rough- — 


ness seems not to have originated until the twentieth cer century. 


——_ Measurements of the resistance to flow in pipes of ‘definite roughness | 
seem: to have been first made by L. Hopf! and K. Fromm‘ in (1923. A ‘rational — 
b. agd of velocity distribution and, from it it, formul: as f for resistance to flow Ww ere 
later developed by Theodor von Karman, ™M. ‘Am. Soe: C. * and Ludwig 


"Prandtl. In 1933, J. Nikuradse® published his” epoch. th- paper com- 
bining Prandtl’s with definite experimental me: surements of the 
in pipes lined with sand | gr ains of size Sine gh a 


ing," C. F. Colebrook": and C. M.W White. 
= ‘The corresponding problem of the « effect of definite roughness on on flow in 


has published a the problem, but with ‘no experimental 
values other than Bazin’s. _ The experiments: described in this paper are a a 


— attempt to do for open channels what Nikuradse did for pipes. 


_ 2“*Recherches hydrauliques,” by H. Bazin, Paris, 1865. His results reduced to aes units are given 
in ‘Flow of Water in Rivers and Other Channels,” by E. Ganguillet and W. R. Kutter, translated into 
English by Rudolph Hering and John C. Trautwine, John Wiley & Sons, Inc., New York, N. Y., 1893. q 
(The temperatures are not given in the translation, but they are evidently in the original as they are given 
3**Laws of Turbulent Flow in Open Channels,” by Garbis H. Keulegan, Research . Paper RP 1161, 
National Bureau of Standards (U. Vol. 21, December, 1938, pp. 707-741. 
eo ‘Elemente der technischen Hydromechanik,”’ by Richard von Mises, Leipzig, | 1914, p. 50, et seq. se, 
5‘*Die Messung der hydraulisc hen Rauhigkeit,” L. Hopf, fiir angewa andte Mathem 
Mechanik, Vol. 3, 1923, pp. 329-339. 
©“Strémungswiderstand in rauhen ‘Rohren,” by K. Fromm, ibid., >. 339- 358. 
ths ““Mechanische Ahnlichkeit und Turbulenz,”’ by Theodor von Karman, Proceedings, 3 3d International 
Cong. for Applied Mechanics, Stockholm, 1930, pp. 85-93. 
8“‘Zur turbulenten Strémungen in Rohren und lings Platten,” by Ludwig Prandtl, E rgebnisse der 
_ Aerodynamischen Versuchsanstalt zu Gdttingen, Vol. IV (1932), pp. 18-29; or ‘ ‘Neuere_ Ergebnisse der 
Turbulenzforschung,” by L. Prandtl, Zeitschrift des Vereines Deutscher Ingenieure, Bd. 77, No. 5 (1933), rs 
_ Strémungsgesetze in rauhen Rohren,” by J. Nikuradse, Forschungsarbeiten auf dem gebiete des inge- 
‘Frictional Resistance in Artificially Roughened Pipes,’ ” by Victor Streeter, Transactions, Am. 
_ “‘Experimentelle Untersuchungen zu zum Rauhigskeitsproblem,”’ by H. Schlichting, 
Vol. VII, No. 1 (February, 1936), pp. 1-34. An abridged translation of the paper is given in Proceedings, — 
Am. Soc. C. E., November, 1937, pp. 16-31 (jollowing 
“Turbulent Flow in Pipes, with Particular Reference to the Transition Region between the Smooth © 
Rough Pipe Laws,” by C. F. Colebrook, Journal, Inst. of Civ. Engrs., Vol. 11, 1939, pp. 133-156. 
‘Experiments with Fluid Friction in Roughened Pipes,”’ by C. F. White, 
Proceedings, Royal Soc. London, Series A, Vol. 161, 1937, p pp. 
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: were of lucite and the base was a structural steel channel with its back upper-_ 4 ae 
al 
| 
a 
a 
— 
= 


vith paint. The of the offered attach- 
ment for the lucite side panels and for ~ vertical light steel angles which | 
supported these panels at 20-in. intervals. | YS The lucite panels were also stiffened — 
— along the upper edge by a lucite flange. © A light steel strap was placed ACTOSS 
the top at the middle of each panel and bolted to these flanges. 7 At six points. 
along its, length the flume ‘Tested on on blocks tom trestles, 


ols ~The definite v wall roughness was provided by square steel strips which 
extended down the sides and across the bottom so as to act as cleats or battens. 
The arrangements of the battens to create the various Toughnesses a are shown 


in Fi ig. larger | of th the tw two ‘being in. square and the smaller, in 


7 


. 


- — 2 ox 10" =20" 


13 ov 


4 5=20" 


The water ‘was measured by ‘a diaphragm orifice in the 6-in. supply line 
a 4 just visible in the upper ‘ right corner of Fig. 1. bor the larger flows the ¢ differ 


ential head on this orifice was measured by the mercury manometer shown at 


‘a the right of Fig. oa, 4 For most of the flows for which the differential head would | 
re have been less than 2 in. of mercury, it was read directly i in inches of water in 
the vertical tubes shown to the left of the mercury manometer. 
A He At the downstream end of th the flume a control gate was arranged as § shown 
in Fig. . 4(). _ This allowed part of the flow to. pass under and Part over, but 
by adjusting the wing m nut (shown. at the: upper right) the depth at the down- 

stream gage could be adjusted sO that the flow throughout the flume was 
oe “ practically uniform. A. thermometer i is shown at the left in Fig. A(e). The 


a ey depth was read by hook gages in wells as — at the right of Fig. 1 ool: 


ad 
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“let of Fig. 2. ‘The term ‘ ‘second slope” ” in the caption of Fig. 4 means that 
‘the bottom slope was 0.0080. The wells were connected by rubber hose with 
the static pressure tubes shown in Figs. 4(a), 4(c), 4(d), 4(f), 5(a), and 6(d). 

‘These w ere made of copper tubing, 0.187 in. outside diameter, with the up- : 


stream end _ plugged with solder and filed to a “‘streamlined’ ia contour. .. 


PRESSURE TUBE 3 


(d) ROUGHNESS RUN 22.3, 
SHOWING PRESSURE TUBE 3 


ify} NO BATTENS, RUN. 23-2, 
SHOWING PRESSURE TUBE 1 


holes were ‘drilled with | a No. 54 drill, 1 in. the | nose—one ‘on on 
: the bottom, and one on each side. It was felt that some such arrangement was — 
_zRecessary because the pressure at an ordinary flush -peizometer opening would 
undoubtedly | have been affected by the proximity of a batten, and, with flows — a 
such as those shown i in n Fig. 5(a) (which are > similar to the flow i in run ee 


‘it woul 


OWING PRESSURE TUBE 3 > y 
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the 
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draft of this p paper on file in the E ngineering ‘Library. Wo 


() OBLIQUE VIEW 
FACING UPSTREAM 


T THE STEEPEST Store 


—EXAMPLes OF FLow at THE STEEPEST SLOPE 


4, W When the w ork was begun the flume was set with its floor at a slope of 
This was enough to shooting flow in the 


CHANNEL ROUGHNESS Dee 
alls. The location of these pressure tubes is shown in Fig. Furth 
gim 
sim 
— | Re. / 
the 
smé 
Aw. : sg wit 
nin 


set u up. After runs h: id been made at this slope with the smooth ‘flume and 
with each of ten roughnesses (with about five different depths. for each rough- 
ness), the slope was ch anged successively to 0.0080, 0. 0020, and @ 0.0005, and 
similar runs were made at each slope. At the flattest slope many - of the flows 
were toi too small to be measured by the orifice m« meter, even with the water-air 


manometer. q For these a a deflector w as arranged at the dow nstream end to 


tN 


Zero 


, in Feet Above 


i Distance, in Feet from End of Flume 
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- direct the flow into a w eighing tank. Ae he time required to run 300 lb of w ater. 
mas: measured with a stop- watch, suet for a few of the sm: allest flows where 
f the: ‘quantity was 1 reduced to 200 i oF even on only 100 Ib; but the time for these. 
‘s smaller Ww veights was never less than 64 sec and was r read to fifths of a second, so — 


that these flows v w ere probably determined more accurately than those measured 


the orifice, 

Before any Tuns were take width of the channel was. measured at 
ions and at two depths in each section with a micrometer 


by. gage niet to thousandths of an inch, The w width s at about: 0.3 in, from the 


4 
SQ 
— 
| 
= 
— 
a — 
— 


= 


7 


CHANNEL ROUGHNESS 


* height of about 4in in. above the floor the width varied ied from 8. 8. eveave to 8. le ) in. ; 
with a mean of 8.166 in. _ In general, “each 20-in. panel was sprung outward 
so that the channel was about 0.03 in. wider at the middle of the panel ‘than 
at the junction with adjoining panels; and these junctions sometimes involved 
an offset of 0. 01 in. to 0.02 in. la Therefore, the so-called “smooth channel” | 
was not absolutely smooth but had a slight | amount of roughness and an 


the orifice meter had been calibrated a 1 few ‘months | before these ie. om 
ments by a a student as result, oft this ‘ealibr: ation on was put into the form of a 


curve of values of Q/V HH plotted against. VH H to o such a a scale that the ec constant | 
could be read to four significant figures. 


described i in detail i in the ‘original draft of this. paper on file in the ean 
Societies ree al 130 _ Them most important data and results are given in Table el 
(Appendix). The slope wa w as determined from plotted energy, gradients, a 
sample of which i is shown in Fig. 7. For the steeper slopes there the 
: ‘possibility that the resistance might depend on the Froude number, - 


i The most fundamental treatment of the resistance to flow i in open ¢ channels 


seems to be that of Mr. t. Keulegan,* w whose formula fc fors smooth walled channels 


which: as and b are constants; v is the kinematic of the 

flowing; B is a function of the shape of the \ cross section; é measures the re- 

- sistance of the free surface; and Vs, the so-called friction velocity, is defined 

n which R= 4RV/v oe the Reynolds number. Mr. Keulegan shows that 


« 


oe 14 ‘‘Laws of Turbulent Flow in Open Channels,” by Garbis H. Keulegan, Research Paper RP br 
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for a a rectangular cl channel ab ab inthe off 
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has the. same value | as that which for pipes, an 
zero, the shape correction term is equal to 
4 


> 


Values of 1 of this term were calculated for various depths and are plotted against 7 
the corresponding hydraulic radius in Fig. 8. - Since this term represents ‘uch 


a small part of the total C, it is felt that a0 


the approximation of making ¢=0 and | 
= 2.5 will introduce | appreciable 


error. By grouping the. terms, rms, § 
| a= ets... | 
in which the logarithm is now to 
The same assumptions as those in Eq. 6 
“would mi makes m = 32. 65. DefineJas: 0 005 010 0.15 0.20 


“Using Nikuradse’ 8 value of a, = 5.5 in addition to these assumptions, J would 


be 32.16 (5.5 2.5 loge 2.01 01) = = 21.96 .96 1 hott 
It may be noted that for a Cross Mr. ‘Keulegan that 
8 = 0.259, so that c, = 14.18 X 0.259 = 3.67, and © | epee nadia ae 


which i is almost identical with the formula for flow i in smooth channels derived — 
“gal However, when C — c, was plotted against logo R/( C, the data for this smooth © 


channel failed to agree with da the first the for 
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the s steepest slope, where the flow was definitely shooting (with F = 2. er - spe 
more), fall in a distinct region indicating an entirely different law. The others Nil 
lie a straight line, but its equation, fitted by squares, es, is sta 
C — = 39.7 logio (10) 4 the 


x 


ing 


Equation 


‘a 


“a 
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‘ Fie. 9.—Vatues or Cufzy Corrricrent C ror a Smooru CHANNEL 


The probable error of m, also least squares, is 1. 4, it was 


Table in the gives the of 
peewee by Eq. 11 (w ith the help of a table of values of C + 40 logio C), and 
also the deviations from the observed C. | ‘The m mean deviation was 2.32 and 
ae the ‘ ‘probable error’ was 2. 0. Since C averaged 117.4, the raven error was 
about 2.0% and the probable « error about 1.7%. 

If the walls are wavy, the value of J in Eq. 7 would be expected to be less. 
- than that for straight Ww alls, but the value of m would be expected to remain 
the same. . Howe rever, the value of m found in these ¢ experiments is definitely 


more e than that found by Nikuradse. e. The value of m= 39.7 + 1.4 “corre: 

= 


hits 


— 

— 

= 

— “Sin 
— pro 
of J from 5.6 to 4.8, so that Eq. 10 anc 
coe 
det 
4 (Ay 
el, 
bs 1939 
— 
’ 


3.04 + ife 
s value of = 2.5 is somew (aotwith- 
standing the papa formulas that have been based on ‘it), but i it can scarcel yo 
be this much in error. Therefore, there are either serious systematic errors in- 


or b ad has a for channels than for The 


this of m. “it may also be that the of velocity. 
- distribution i in open channels by E. H. Taylor, 16 Jun. Am. Soe. C. E., following 
a modification of a method proposed by M. P. O’Brien,” M. Am. Soc. C. E., : a 
gave an average value for the “universal constant” of 0.306, which corresponds | 7 
b= 3. 27 and which is even more than the writer’s experiments show. 
As stated, Eq. ll does not represent the runs at the steepest slope (setups a, 
12 2, 3, and 8, in Table 1(b) in the Appendix). When the discrepancies from 
Eq. 11 were plotted against the Froude number F, there was no apparent 
relationship for the flatter slopes | (setups 23, 24, and 46, in Table 1(a), i in which 
F varied from 0.26 to 1. 69) ; but on the steepest slope where F varied from 2.49 
“to 3.28, there was a very definite dependence on F,w which it seemed n might be 


well expressed bya a formula of the 


olving for the most probable values of J, 2, and least ‘squares gave: 
= 30. re 2. 2;2= 58.7 + 3.1; and = 171.5. In view of the - 
: a probable errors it seemed best to round these values tox = 30andz = 6 / 


— which for the best fit numa changing J to — - 175 and gave the formula 


The C computed by Eq. (13 (with aid of a table of values of 
C+ 60 logis C’) and the deviations from the observed C are given in Table ry 
in the Appendix. The mean deviation Was 2 2.44 and the probable error 

ince C ‘averaged 117. 9, this r makes the enone about 2 tien ae 


robable error | about 1. 


case (no was taken, but the was similar to that, 
™ 4(@), it is not ; surprising g that F enters into this case although it did not 
enter v when the flow wa as more ‘tranquil. It is surprising, howev er, that the | 
. coefficient of logi R/C should change. it is also surprising that the change 7 
= Eq. ll to Eq. 13 came between F = 1.69 and F = 2.49 instead of at. 
F »'? 00. (Tests at one or more intermediate slopes would be required to ¥ 
determine y Ww here the change occurs. The sis six runs of setup 28, Table 1(a) 
(Appendix), had Froude numbers: F varying 1g from 1.45 to (1.69 ‘and agreed 


AE 16 “‘Velocity-Distribution in Channels,” by E. H. Taylor, Transactions, Am. Geophy sical Union, 
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with Eq. 11 with an average discrepancy of of 0. 87, whereas the the average dis- 


‘Mr. Keulegan’ for walled in the notation of this 


2 and 


essrs. 
and Keulegan ré represented by k). "Substituting the value V; and [rearranging 


€ 


Asi in the case of danslnith chieunis (see Eq. 6), can be replaced by ¢, and 
the other terms combined to make 


‘For roughness values of C — c, were plotted against Since 
: ‘roughness IX (Fig. 3) constituted a half-seale model of roughness I except for 


penalties IX, with hydraulic radii twice those observed, were pense plotted on 
the same sheet. © Similarly, roughness VIII was combined with roughness | II, 


roughness VII with ‘roughness Ill, and roughness’ X with roughness XI. 
_ each case the slope of the line m= 40 seemed to give as good a fit as any. 


‘This value was obtained visually without a least squares solution. 


taken as 0.25 in. = = 0.02083 ft, = 48, and 40 logi 1/e 


would be 72.2 — 67. 25 = 4.95. On the other hand, 4. 95/40 = 0.12375 
= logo 1.327 and the reciprocal of 1.327 is 0.754, s o that if « had been taken 
as about three quarters the actual height, J would have been zero. Since the 


measure of roughness depends not only upon the height of the particle, but 
also. on the spacing, it is obvious that the unit of measurement of roughness i is 


7 quite arbitrary, and that it would be quite allowable to take it such that J = 0. 
Therefore, the following formula was used: 


a! = 40 logio logio R — 40 log . (17) 


a cg ns From the a aforementioned plottings, and similar ones for roughnesses IV, 
nd VI, was obtained and then Cc was computed for each run and the dis-_ 


— 

Tab 
— betw 
— juste 
Lan. 
= 
proc 
= ft 
were 
the } 
erro: 
ee velo 
1.99 
This 
— 
— char 
— Prot 

la 
No: 
crep 
plot 
| relat 
— valu 
22.5 
“othe 
valu 
max 
to re 
hii ] 
aig ‘Teco 
— and ODServed Values 


— 


ii 


December, 1044 CHANNEL ROUGHNESS 
Table 1(c) to Table 1(m) i in the: Appendix. ot If there v was any | serious difference 
between the sums of the plus and minus s discrepancies, the value of « was ad- 
_ justed t to make the discrepancies balance. _ The final values are listed in n Table 
2(a) in the Appendix. Table 2 also gives es the values of h, s, s/h, in which 
his the actual dimension of the batten, and s is the spacing of the battens 
center to center. Since roughnesses IV, V, and VI consisted of both g-in. and 


tin. no value of h was assigned tothem. 


ling to probability theory, ‘the percentage error” of the 
product | or quotient of two observed quantities i is the square root , of the sum 
a the squares of the percentage probable errors of the factors. ; Estimates | 
were made of the probable error in each of the various kinds of readings, and, _ 
from them, of the probable errors in the computed quantities, based on the 
aforementioned probability theory. In the case of the slope, the deviation of 
the plotted points from a straight profile also gave an indication of the probable | 
error. These estimates gave percentage probable errors, as follows: For the 

velocity, 1.5%; for the : slope, 1.7%; for the hydraulic radius, 1.3%; and for C, 

1. 9%. _ These « are only estimates an and may be too low. © At any rate the aver-— 
age deviation of the 200 runs on rough channels from Eq. 17 was 4. 37%, , which | 
by C. A. F. Peters ’ formula (probable error is approximately 0.845 times the 
average deviation) would indicate a probable | error of 3.7% in one observation. — 
This includes any error due | to the choice of a wrong general form for Eq. 17, 

" any error in the value m = 40, and any error due to the wrong choice of € for 
each particular roughness. | As to the value of m for tranquil flow in the smooth 

channel, it already has been noted that the least squares solution gave a 
probable value of 39.7 + 1 A, and as the a average value of logio R/C was 3.441 . 
the probable error in the term 40 logio R/C might be estimated at 1. 4X 3. A4l 
= 4.8; but the error in the computed value of C would not be this much as it. 
is largely compensated for by choosing a value. of J to o correspond t to the m used. ; 

No study was made of the accuracy of the value m = 40 in Eq. 17, hae from 
‘the plotting the writer would judge the probable error to be about 1.5. 

_ The shape correction, ¢,, is of the same order of magnitude as the dis- 
-crepancies, so° that the data are insufficient to check its correctness; but a 
‘Plotting of the discrepancies against values of the depth y : showed no ‘systematic: _ 
relationship, so that there is no evidence against its correctness. it may be 
mentioned again t that the relation c, = 14. 188 was based on ‘Nikuradse’s — = 
value of b. The writer’s value of m = 40, instead of 32. 65, is an an increase of = 
22. 5% and would increase the coefficient of B from | 14.18 to 17. the 


other hand, é, the free surface effect, was taken as zero. Ifit has an appreciable 


maximum change i iN Cy could be. only 6, > that it was. thought unnecessary 

general, it must be admitted that the > accuracy of the 
--Tecorded herein is only of about the same order as those of Osborne Reynolds"* — 
Experimental Investigation of the Circumstances Which Determine Whether the Motion of 


P Water Shall Be Direct or Sinuous, and of the Law of Resistance in Parallel Channels,”’ by Osborne or J 
Philosophical Transactions, Soc. of London, Vol. ‘174, Ill 
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CHANNEL ROUGHNESS 
on pipes in the 1880's se it can only be hoped that, in this case too, they may 
lead to more accurate measurements by observers w up 
7 mM %. It might be urged that, since quite large standing waves formed on 1 the 
surface at the higher velocities in the ‘rough channel, a term involving the 
_ Froude number ‘ should be included i in E Eq. .17 . However, a a plotting of the dis- 
-erepancies from Eq. 17 against F showed no no_systematic relationship. This 
as to be expected because, since | =V gy and C = V R RS, and since 
a Rva varies nearly as Y, Tesistance depending on F could hardly be distinguished 
other resistance except by the difference in slopes, and in these tests 
“there as only one ‘slope a at which F was large ge enough to much effect. 
12 of Bazin’ teste" ere on plank channels with 


that the were not square. it seems that this 


_ relatively little difference, Table 1(n) (Appendix), has been computed. The 
of which gave the best fit in Eq. 


= was selected (and listed i in Table 2(b) 


in the Appendix), ‘and ‘computed 
and the deviations from the observed “i 


are giv en in Table 1(n). mean dis | 

erepancy 1 was only 1.6, as compared to 

The ratios of the rough 
ness to the actual height of the ‘ ‘rough- 
particle” and the ratio of spacing 
to » height, as. given in Table’ 2 in the 

Appendix, a are plotted in Fig. 10. The 
dashed portion of this curve is ad- Its 
mittedly uneertain, but it indicates half 
that the spacing of square battens 
J} woul ould give the the maximum resistance is. oe 
< 02 06 about twelve times the height of the 


10.—RELATION Rovas- roughness” is about 0.83 of the > height. whe 
It must be remembered that the scale of Fig 
“measurement o of “equivalent ro t roughness” was s taken ai arbitrarily to o make J in J in Eq. bis 
x - 15 equal to zero, and that therefore this last result would not necessarily a agree 2 
with Nikuradse’s sand grain size and his formula for pipes. The facts that: an 
a _ Schlichting” | obtained a value of e/h = 11.9 for the flanges of a small angle s sec- and 
9 ate a tion placed across his pipe, and that Mr. Streeter » found values of ¢/h as high Rey 


as 4, .57 for square circumferential grooves, both investigators ‘using € € equal | to 


Nikuradse’ sk, indicate | that the unit of € used in Eq. VW may be, say, from 


five to ten tes as large : as that used by Nikuradse. Although a standard 
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value wold be desirable, it would not seem feasible to a dopt one until ae 


one geometrical type of roughness (such as the one used in this study or some 
other that could be reproduced as as easily) has been carefully calibrated by the — 
‘Schlichting m¢ method ii in a pipe, and by some method comparable t to that of the = 


writer, in open 1 channels, ‘Then all other roughnesses could be referred to this — 


ay The discrepancies of Eq. 17 from the observed data a average - nearly twice 


as much as would be expected from. the probable i inaccuracies | of the measure- 
ments a as estimated by the observer; but these are only estimates and it may 
be that Eq. 17 is entirely correct and the > discrepancies entirely due to errors in 
observation or to errors of observation plus errors in the chosen values of 
It seems more probable that the discrepancies are partly due to these som 
factors and partly due to wrong values in the constants of Eq. 17, but that the | 


form of the equation is satisfactory, at least for tranquil flow. ws air 


For ' the ca case of the depth half the width, c, is 2.74. on resulting formulas" : 


For channels, 


= 


and, for channels, 


For flow, that of the hydraulic gradient i is the same in 


C= = 
57.0 
lt should be 1 be noted that Eqs. 18a, 18, ‘and | 18¢ apply. only when the depth is 
“half the width. h. For smaller ‘depths the constant in 18¢ is 


ties the width is infinite compared with the depth ogc 18 are ® plotted i in 
— “Rig. 11, which sums up most of the results of this paper. I 
stood that, as in the corresponding case for pipes, the same may act 
as rough orsmooth. F or example, a channel 4 ft ft wide flowing 2 ft deep, with a. 
an absolute roughness (e) of 0.0010 ft, would have R= 1.0 ft and R/e = 1,000, q 
and C would be 122.7 for any Reynolds number greater than 92,900; but, for 
Reynolds ‘numbers less than this, it would behave as a smooth channel with Cc. 
decreasing as R decreased down to about ¢ 67 at R=2 ,000, when Eq. ‘ibe 


om Review of Flow in Pipes and Channels,” by S. J. i and C. M. White, Engineering, Vol. 23 
» Pp. 71. 
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uncertain, has been studied by L.. Hopf,” H. E. Herten, 
_M. Am. Soe. C. E., the late H. R. Leach, Assoc. M. Am. Soe. C. E., and R. 
Vliet,” and by L. G. ‘Straub, M. Am. Soe. E. Professor Straub, 
however, used a different definition for Reynolds number and also seems § to 


. have used the constant for a round pipe in his equation n for laminar flow. 
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11. é AND R For Recriwovtsn CHANNELS 
WITH THE DeptH To THE WIDTH Va! 


By analogy with pipes one may suppose also that the region just to the 
oy right of the smooth channel line is probably ‘not given correctly by the lines. 
re as drawn. One would expect that for uniform roughness there | would be an 
— ; upward bend (just as Nikuradse® found a downward bend for the friction 
eo coefficient fin pipes) and that for irregular roughnesses the bend would be in 
the opposite direction , as f found by Messrs. Colebrook and White. This 
transition range in pipes has ‘been discussed quite fully by. Mr. Colebrook” 

| by Hunter Rouse M. Am. Soc. C. E. i Unfortunately, the w riter’s tests 


PA: did not include any case in which the rough channel acted as smooth. it 


20 **Turbulenz bei einem Flusse,” by | L. Hopf, Annalen der Physik (4) ) xxxii (1910), pp pp. 777-808. 


= 21‘*The Flow of Water in an Inclined Channel of Rectangular Section,” by Harold Jeffreys, Philo- 
am ae _ sophical Magazine and Journal of Science, London, Edinburgh, and Dublin, Vol. 49 (May, 1925), pp. 


i 2 "Laminar Sheet-Flow,”’ by Robert E. Horton, H. R. Leach, and R. V an Vliet, Transactions, Am. 


23 ‘Studies of the Transition Region Between | Laminar and Turbulent Flow in Open Channels,” by 
6. Straub, tbid., 1939, pp. 649-653. ky) Pole 


%“Evaluation of Boundary Roughness,’ by Hunter Rouse, Proceedings, 2d Hydraulic Conference, 
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would b be a ‘difficult, condition to arrange fo for W water, _ but very easy | with a liquid 


& In all the computations of this ‘paper the area and hydraulic radius are 
“based on the ; gross section of the channel without the battens. . ae the battens 7 
were quite close together there certainly would be a question as to whether 7 #8 
- computation should not be based on the net section, or on some kind of average — - 
section. This knotty question is left for further study. ; 
—o«d—f Fig. 11 correctly represents the flow of liquids in channels, the implica- a 
tions for river model studies are evident. if the prototype | acts as a rough 
© tue it is not sufficient to give th the model the same ‘telative e roughness, 
because, at its smaller Reynolds number, it might be acting | ‘as a smooth © 
The only solution seems ms to be to use a ‘sufficiently large model to 
bring iti into the rough channel 1 range. 
The writer wishes to express his grateful appreciation of the | generosity of | 
7 low: a Institute of Hydraulic Research w hich placed its facilities at his 
disposal, and met the incidental expenses of the study; especially to F. M. 
Dawson ,M - Am. Soc. C. E., Director of the > Institute, KE. W. Lane, M. Am. 
Soe. C. E., , then Associate Director in Charge, and to Hunter | Rouse, M. Am. 
Soc. C.E. and ALA. Kalinske, Assoc. M. Am. Soc. C. E., and to A. McLeod. 
He also whkes -to thank Paul Nemenyi and James Robertson, Jun. Am. Soc. 
CE, at that time | both connected with the Institute, and to express his appre- 
ciation of the work of E. W. Thomas and Dale Harris, technicians, who assisted 
in the construction of the apparatus. x Several Gudea, working under pro- 
visions of the National Youth Administration, assisted with the computations. 
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symbols, used in the paper, conform to American 
Letter Symbols for Hydraulics,* prepared by a Committee of ‘the 
American Standards Association, with Society representation, and approved 


4 a = aconstant ‘used by Mr. Keulegan,? subscript r denoting rough channels 


“ti and subscript s denoting smooth channels; 
8 2 constant, the reciprocal of von Kaérman’s “universal constant”’ ; 
Ba = the Chézy coefficient i in the formula Zz = Cy RS; 
=a correction; c, = a small all shape correction; 


= coefficient of pipe friction in the formula S = {_— 


= head on the measuring orifice; 
he height of roughness particles, or battens; 
a constant defined by Eq. 8b; Pata 
k aid symbol for for roughness used by Messrs. and Keulegan 


a 


R= 


R = Reynolds number; 
slope. of the energy gradient; 
= of Toughness battens (center to center); 
velocity as defined by Eq. 


: — the unknown coefficient of F in Eq. 12; 


of flow; ; that i is, vertical or measurements 


= an unknown coefficient i in Eq. 


a function of the shape of the cross ‘itil as defined by Eq. 5; 


: absolute or equivalent roughness; that is, a measure of the Toughness 4 


of the channel walls (see k); 


és = = 8 measure of the resistance of the free surface Bs », 


= - kinematic viscosity of the flowing liquid; 
7 = shearing stress; tr. = shearing stress at the — aia 
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CHANNEL ROUGHNESS 
TABLE 1.—Osservep Test Data Aan 
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0.148 | 0.103 
0.187 | 0.121 
5 | 0.239 | 0.140 
| 0.298 | 0.159 
7 | 0.349 | 0.172 
0.410 | 0.186 
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0.109 | 0.083 
2 | 0.159 | 0.108 
0.217 | 0.133 
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46- 98 104.1] 0.5 
106.4 | 0.30 | 107.2} .... | 0.8 
1209 _ | 0.32 116.4} 4.5 | .... 4 
6 129.0] .... 1.4 
132.9| 3.2 |.... — 
5 1369) ....1 08 4 
34 | 74 | 1370 | 1.57 139.9 “3.9 
154 | 143.1] 0.0 |<... — 
— 
127.9) 12 ].... | 
1276|....| 02 — 
| 128.0| .. 0.9 — 
| 0.271 | 0.151 | 8.24 128.0} 4.9 | 
31 | 0.100] 0.077] 5.22. | 291 107.0| 
3-2 |0.130 0.094 | | 208 | 2.9 .... 
10.137] 0.098 | | 308 2:9 41 
| 0.107} 6.40 | 281 7 | 2.8 117.3] 2...) 06 
35 | 0.180] 0.118 121.0] .... | 21 
3-6 | 0.195 | 0.124 124.1) .... | 5.2 
3-7 | 0.246 | 0.143 123.3) 3.9 | .. an, — 
3-8 | 0.257 | 0.146 123.0] 3.8 — 
0.279 | 0.153 
83 | 0.1 
8&5 [0.1] 0.111 — 
8-6 | 0.173 | 0.116 
| 0.184 | 0.119 
$8 | 0.193 | 0.123 
| 0.204 | 0.128 
| 0.233 | 0.138 6.4 
8-13 | 0.3 120. | 
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TABLE 1.—(Continued) 
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0.247 | 0.143 
0.259 | 0.147 
0.272 | 0.151 
0.282 | 0.154 
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0.093 | 0.073 | 0.198 
0.130 | 0.094 | 0.243 
0.199 | 0.126 | 0.368 
0.240 | 0.141 | 0.345 
0.317 | 0.164 | 0.509 
0.391 | 0.182 | 0.556 
0.092 | 0.072 | 0.410 
| 0.134 | 0.096 | 0.612 
0.168 | 0.112 | 0.718 
0.238 | 0.140 | 0.874 
0.277 | 0.153 | 0.911 
0.339 | 0.170} 1.03 
0.108 | 0.082 
0.134 | 0.096 
0.150 | 0.104 
0.178 | 0.117 
0.190 | 0.122 
0.247 | 0.143 | 
0.255 | 0.146 
0.313 | 0.163 
0.319 | 0.165} 1.93 
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| —(C ontinued) 


R Ss , | (degrees 


0.171 | 0.114 | 2. 45.6 |1.06| 28.4 
0.236 | 0.139 | 3. 310 1.10 4 : 
0.287 | 0.156 | 3. «46.3 | 1.06} 
0.331 | 0.168 | 3.41 | 310 3 | 104] 28.6 


0.067 27.38 10.10] 21.3 
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0.077 | 0.346 21.0 27.8 
25] 0.6231 20.4 | 39.0 
0.153 | 0.756 | 218 | 418 
1941 0.095, 1.03 | 75 | 
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0.101 | 43.3 | 0.54] 30 | 2 43.5) 
0.53 | 12:1 | 2:7 | 460] 06 = 
26/256 | 309 | 411 26 | 43.1] .... 
| 309 $3.0 f 25.3 | 2.7 | 46.4 
VI (Fic. 3); Discrepancies BY EQ. 17 2 
Roveuness VI (Fre, 
39-2 | 0.127 | 0.05 5.2 | 448 | 0.14] 22. 210 | 26 | 466 
30-8 0247 |O1a3 | O380-+ | | Oud | 316 | 27 | 492 36 
| 0.247 | 0.143 | 5.4 221 | | "0 | 271 36- 
39-5 | 0.321 | 0.165 | 0.494 | 6.9 | 50.7 | 0.16 | 22.1 36- 
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(62 | 52.7 1016} 2.9 | 38 
0.386 22.6 1.29 1.6 411 7.9 
0.412 20.2 33.2 | 0.23 22.6 | 19 | 45.3 
0.133 | 0.096 0.595 | 20.2 42.8 2.4 | 50.8 
.133 | 0.096 | 0. 1 | 0.32 26 | 534.1261 538 3.0 
4 he O31) 226 | | 26 | 538 
0.268 | 0.15 909 | 526 | 031 27 
0.103 | 0.079 | 1.20 475 | 0.06 | 25.9 6.39 | 2.0 
|0.147| 0.103} 1.47 | BO 26.0 9.05 | 2.4 
| 0:196 | 0.124 | «| 200 | 25 
0.306 | 209 | 83 | 572 |oe7| 200 | 
0.306 | 0.161} 2.09 | ; 8 
1.09 27.6 | 8.23 | 1. 
0.119 | 0.088 | 2.12 1.09 2:1 
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"Rune (ft per served | F puted 
(k) Roveuness VIII 1G. 3); DISCREPANCIES BY Ea. 17 
37-1 | 0.083 | 0.067 | 0.254 2 0.640 45.6| 0.2 |....— 
37-2 | 0.130 | 0.094 | 0.278 0.992 | 19 | 52.0] .... | 
37-3 | 0.181 | 0.118 | 0.422 8 | 189 | 231 564 
37-4 | 0.257 | 0.146 | 0.543 2.6 | 60.4 0.5 
37-5 | 0.322 | 0.165 | 0.628 4 | 395 | 2.7 | 62.6 340 
37-6 0.388 | 0.181 | 0.724 4 4.99 | 2.7 | 64.2 2.5 
0.097 | 0.076 | 0.561. 1.60 | 1.5 | 47.8 3.0 
| 0.133 | 0.096 | 0.706 2.55 | 19 52. Ls 
0.202 | 0.127 | 0.926 | 446 | 24 B77 
0.245 | 0.142 | 0.962 | 5.26 | 26 | 59.9 3.6. 
0.320 | 0.165} 114 9 735 2.7 = 
0.372 | 0.178, 7 18.1 | 2.7 0.9 
10.084} 2.60 | | O48 117 | 
(10.1511 0.105] 3.20 | : 1.24. | 
3 | 0.195] 0.124 | 3.72 | 20.0 | $3183) 
0.246 | 0.143 | 4.12 | 2 | 256 | 2.6 | 60.0| 17 
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0.268 | 0.150 | 4.60 | 312 | 1.56 
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q ENERGY FLOW EXPANSIONS 


Information on the energy changes that occur when fluids flow i 


in expanding 


a 


conduits is presented in in this paper. Specific data are presented on me mean 
6 “velocity distribution, turbulenes, and pr pressure changes for 7.5°, 15°, 30°, a and 
180° expansions from 2.75-in. to 4.75-in. pipes. The velocity and turbulence 
measurements were made photographing, on motion-picture film, the flow 


in ‘transparent conduits. The water therein had “suspended in it small, 


af 


miscible droplets w which formed streaks on the pictures, from which the direction 


A Deer 


and magnitude of the velocity at an any y instant in two dimensions could be 


Da The analysis of the data for mean velocity distribution, variation of turbu- 
lence, and pressure changes gives considerable information on the /mechanics 
the energy conversion : and dissipation Processes. Iti is found that the total 
loss of energy is considerably greater | r than the total « energy of turbulence that oo : 


is produced, thus indicating that a major portion of the > energy is is lost by direct — = 


conversion into heat at the regions of ‘high, local shear in the fluid. The 


ane 


: 


@ 


toto: 


turbulence created is a by-produet of the energy dissipation process; however, 

it plays an impor ant réle in the shaping of ‘the mean velocity distribution ~_. 


—— ore in and be ond the expansion. Data are resented on the ener + loss 
a) and efficiency of the energy conversion for two rates of flow through the four ye: 


‘ Examples of fluid flow in conduits where the cross- sectional area increases, F a 
with, a resultant decrease in velocity and increase in pressure, are common 
“wherever fluids are carried and controlled. Specific examples Tange all the 
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EXPANSIONS 


and the loss of head in any such may ‘be only of ‘miner 
importance to the engineer; or it may be of considerable importance. — Through- 
out the years, considerable data have been accumulated, relating to losses in 
various types of such expansions, s, and in general hydraulic engineers have a 
fair idea of the pressure losses involved. Therefore, this paper is not concerned 
with the mere presentation of further pressure loss dete, but rather with supply- 
ing data that w will shed some light on the internal mechanics of the energy 
“* changes that Occur in the fluid flowing i in and beyond expansions. eye. 
a The experimental data presented herein have been | an outgrowt: th of a study 
relating to the general problem of the hydraulic ¢ conversion of kinetic to poten- 
tial energy, which was started by the Society’s Committee on Hydraulic Re- 
‘search in 1935 at the Iowa Institute of Hydraulic Research. In order that 
quantitative information might be available relating to the conversion ‘process, 
it was necessary that all the elements entering into the process be measured, 
‘Therefor for each expansion studied the mean velocity distribution “was 
determined at several ‘Sections in and beyond the expansion. — Since the con- 
version process involves transformation of a part of the kinetic energy into 
turbulence energy, it was desirable to have some measurements re relating to the 
intensity of turbulence at various ‘sections. This, of course, proved to be the 
most difficult part of the experimental aie: however, some very significant 
data were secured. _ The change in the pressure along the conduit in which 
the expansion was installed was also determined. 


a Before the experiments | and results are discussed i ‘in detail, a brief review 


and pertinent past work on this subject. Also, since certain items relating to 
fluid turbulence are basic to this particular study, a short. description will be 


: given . of the relation of turbulence to the energy changes that occur in flowing 


ais letter symbols, used in this paper, conform essentially to American 

‘Standard Letter Symbols for ‘Hydraulies,? prepared by a Committee of the 

American Standards Association, with Society ‘representation, approved 


B, Mecuantcs or Fiurp Frow 1n Expansions ty 


“ad From a practical standpoint the principal item of interest is the pressure 
that occurs W fluids flow through expansions. However, in order 


along the conduit, it is necessary that the flow mechanism be 
in detail. This p paper is concerned only with incompressible. fluids; therefore, 


& ‘density and temperature changes will be 2 neglected. 
ee Visual: observations of flow in conical expansions reveal that there is a1 


3 appreciable i increase in in the angularity of the velocity of fluid particles near the 


boundaries as they p pass through an expansion. In other words, ‘the flow 


a 


will be g given of some of the fundamental facts known about Save in expansions 


as 
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becomes more eddying, ‘and the eddies created near the boundaries ; gradually 
Pe into the center stream which for a while after leaving the smaller conduit _ 


is relatively y free of turbulence. — The eddies coming from the wall naturally 
E have relatively low forward velocities, and, as they intermingle with the rapidly — 
28 town central part of the stream, the fast-moving fluid is gradually slowed 


1ed down. As the mean velocity ‘decreases there is a gradual pressure rise, and 
ly- the efficiency of the conversion process of kinetic to potential energy depends 
ey ~~ how - much of the initial kinetic energy is used in forming eddies and in 
idy 2 The loss in. any y expanding. section can be divided into two parte—that due 


en- [— to friction on the boundaries and 1 the extra loss due to the shape and form of © 
the conduit. Except for slow expansions (those having total angles less than 


about 10°) the frictional loss is negligible compared to what might be called | 


hat 
4“ form loss.” ” Both analytical work and experimental data indicate that the 


| total expansion angle giving the least loss varies between 5° and 10°, the exact 
was value being dependent on the velocity ratio and the initial velocity distribu- - 
oie os As the expansion angle is increased, the form loss increases rapidly, as 

it is associated with the flow separation that occurs at the boundaries; that ‘is, 
the as the expansion angle becomes large, the fluid near the boundaries — 
the have sufficient ki kinetic energy to move against the pressure gradient; instead 
cant a backward flow occurs. There is thus created a . discontinuity inside the fluid 


hich | with the Tesult that high shearing : stresses occur in the fluid and violent turbu- 


view Some extensive data on the loss in expansions varying from 3° to 180° 
expansion) betwe een circular pipes of 1.5 3.0 in. ‘in diameter were 
- obtained by A. H. Gibson? He expressed the loss in terms of the Borda lo 
“for sudden expansions, and showed that the total loss was a minimum for an 
expansion. of 5.5°, being about 13.5% of the Borda loss; it reached a 4 maximum | 


sions: 
to 
ll be 


wing 
a of 121% when “the an le was 63°, and then diminished to (102% for a a 180° : 
expansion. Extensive experiments on the loss in sudden ext ‘expansions “were 
in which Um are the mean velocities in the upstr ream and 
In his work conversion of energy cross-sectional divergences, 
1 a Peters® * emphasized the influence of the conditions of inflow on the efficiency 
takes of the conversion process. He stated that complete conversion requires a 
aie discharge length beyond the expansion which depends on the angle and the ~ 
otaah velocity distribution at the entrance to the expansion. | He defined the effi- "i 
i: sca of conversion as the ratio of the actual rise in pressure, corrected for wall 
flow by W. H. Am. Soc. C. E., Vol. LXXVI (1913), p.999. 
_5“Conversion of Energy in Cross-Sectional Divergences Under of ” by 


Technical Memorandum No. 787 (Translation), National Advisory Committee for 


rs 
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ene: was calculated by net of the actual velocity distribution ei 

any section. — ‘Mr. Peters showed that superimposing a spiral motion \ on the 

— flow entering any expansion. increased the efficiency markedly. It a appears 
that the rotational motion delays separation, and thus decreases 1 the energy loss. % 

ae . Kratz a and J. R. Fellows’ conducted detailed experiments on pressure 

; ducts, and | among their findings i is the conclusion | that diverging s sections having 
angles than 60° give losses at least as as those in a 


sections angles of 5°, and it is apparent from his data that the 
shape of the velocity curves depends on the angle of divergence, the curves 
‘becoming more pointed for the larger angles. He noted that separation 
occurred at a 5° divergence and that, for slightly | larger angles, there was a 
__ tendency for asymmetrical velocity distribution and a we of the nal 


_ Some analytical work has been done i in regard to coimpiied and pea 
flows by extension of the knowledge relating to turbulent boundary layers to 
s where a pressure gradient exists. Calculations: for such conditions 

‘indicate whether a given velocity distribution and p pressure » gradient are possible . 
“without separation of the boundary layer. _ A detailed di discussion of this - 


Certain Concerts RELATING TO TURBULENCE 


“guid ts turbulence, an ana of some of the basic ideas relating rn 
turbulence is enmutial. _ Certain aspects of the statistical theory of turbulence 

have been ‘presented in a previous paper. The important items relating to 
turbulence are its origin, designation of its intensity and energy, its diffusing 

& } In general, for various dynamic conditions, re regions s of high shearing stress, 


“eddies” are ‘created. “Such regions can. solid in 

ing a slower-moving or still fluid, and 

- surfaces between two fluids of ¢ different density when a sufficient relative 
exists. As the vortices or eddies develop, they move away and enter 


other parts of the fluid. d. Although the motion of an isolated vortex, such as of 
the loos helical type, might | be: analyzed, when large numbers of such vor- 


“Pressure Losses a from icles es in Cross-Sectional Area in Air Ducts,” by A. P. Kratz and 

R. Fellows, Bulletin No. 52, Univ. of Illinois, Urbana, 1938. 

Bc’ “‘Untersuchungen fiber die Strémungen des Wassers in konvergenten und divergenten Kanilen,” 

by J. Nikuradse, Forschungsarbeiten des Vereines Deutscher Ingenieure, Heft 289,1929. | 

**Aerodynamic Theory,” W. F. Durand, Editor, Berlin, 3, 1935, 155 
#**Relation of the Statistical of Turbulence Hydraulics,” by J 
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. that it becomes impossible to predict the behavior of a pestienlen: eddy. When 
a mass of fluid has a large number of eddies of varying size and strength, itis — 
aa to be turbulent. tt has been accepted by va various workers in the field of 
- fluid turbulence that the term would be used only to describe ‘fluid motion 
that is entirely haphazard i in character and which can only be described quan- - 
titatively i in the statistical sense. _ Thus, the w vell- ordered system of vortices r 
one sees immediately behind a cylinder immersed in n a moving fluid ‘stream is oll : 
rey ntensity and Energy of Turbulence- —The notation to be used for velocity 
is as follows: The velocity at any instant will be designated by the three rece 
tangular components U, V, and W. The direction of U will be in the direction Ls x 
of mean flow, or, in this paper, along the conduit axis. The components 
W will ill be at right angles to U . The value of U can be considered as equal o 
U0 +u, u, in which U is the mean axial velocity with respect | to time at any _— 
r point, and u is the fluctuating part due > to the turbulence. Since the mean 
~ flow is in the direction of U, V and W are zero. _ The turbulence ‘components 
at any instant perpendicular to the ‘wcomponent a are v and w. Borrowing the 
concept of ‘standard deviation” from statistics, the e intensity of the turbulence 
in each of the three directions will be defined as Vy uw, Vv, and Vw BS) omit 
ot The total kinetic energy of the fluid at any point in the section, —- 
. in foot-pounds per pound of fiuid, will then be equal to Ve citer a5 


in hydraulics a as s velocity head), and will be designated bye The s sum of 
the other terms is the kinetic energy of the turbulence and will be —r 7 


Wty 


tr item to. a in mind is that similar values of v , and w? for owe flow 
s ‘conditions do not necessarily mean that the turbulence is identical. The ae 
intensity of the turbulence only characterizes the rotational velocity of the ba 
eddies. For a complete description of the turbulence, some measure of the 


eddy size is necessary. diffusing power of ‘the turbulence is 


directly proportional to both the intensity and the scale. 
‘Turbulence Energy Diffusion Dissipation—For : a fluid moving ur under 
os _ shear, the energy producing the flow i is continually being transformed into heat. 


_ The work done by the shear, per unit volume of fluid, per unit time, for parallel Eat 
flow, is equal to —— , in which 7 is the unit shear ‘stress, and dU /dy is the 


y 
mean gradient transverse to the flow. For viscous flow this rate is 


to , in which is the of dynamic viscosity. In tur- 


‘equivalent to the rate of. energy y conversion into heat. Part of the work 


: by the Shear results i in the creation of turbulence; that i is, the energy producing 


flow is in part transformed into energy of turbulence, which i is ‘eventually con- sees 


4 


Bae flow the rate at which the mean shear does work is not necessarily — te a 
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verted into 


heat at through the action of viscosit course, in the 
3 


4 acter. A general dissipation of energy i 


ae has been derived by G. I. Taylor” and it contains terms: of the ty ype t 
(2 By , ete. These terms are the mean of the “squares of the 


e 
‘ instantaneous velocity due to the fluctuating velocity components. 

For isotropic turbulence, @ =v = ‘such a as that generated i in the ‘throat 


of a Ww ‘ind tunnel due to a regular system of grids, the rate of turbulence energ 4 
a dissipation, as shown by Mr. Taylor, reduces to a simple expression because —# 
of the ‘simple relation that exists ative een the fluctuating velocity g gradients i in 7 


the different directions. of dissipation per unit volume is then 
OF course, in conduit flow the turbulence is far from. isotropic. The prac- 
a thea significance of Eq. | 3 is that it gives a clue as to what factors affect the 
, “rate e of turbulence energy dissipation. ' The term (0u/dy)? can be set propor- 
tional to to u2/d2, in which d is a scale of turbulence associated with the dissipation - 
process. ‘Thus, as the size of the | eddies decreases, for a given intensity of 


turbulence, rate of dissipation increases. Therefore, can be 


> 


> = 


= 
ono 


the average scale tends to increase, since the sa satien are dissipated 1 more 
‘rapidly. The larger eddies are “ground down” in a process that is not as yet 
clearly understood. fe There seems to be a , stretching of the vortex filaments; 
for instance, Mr. Taylor shows that the rate of dissipation by 


the vortex with a resultant decrease in size of 
but an increase in the angular velocity, which follows from the law of conserva- 


ofangularmomentum. 


oa The foregoing | ideas will serve as a background | for analyzing the ee “ 


“EXPERIMENTAL AND ANALYTICAL PROCEDURES 

Ai ‘complete study of the energy conversion and dissipation processes in the 

- ie in expanding conduits involves a study of the mean pressure change 


_: 19 “Statistieal Theory of Turbulence,” by G. I. Taylor, Proceedings, Royal Soc. of London, Vol. 1514, 
Di ipation of V. Vortic in a’ Turbulent Fluid, by I. Taylor, ibid., » Vol. 164A, 
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ag directly converted into heat since the value of dU/dy is very high in this J 
region. Since turbulence is so diffusive in character, the turbulence energy 
ereated at one point is not necessarily dissipated into heat at that point. a 
dissipation of turbulence energy is an extremely complicated phe- 
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of the turbulence. So far as the turbulence is the 


obtaining of complete information regarding it presents formidable experi- 
- mental difficulties. It was decided that, if its intensity could be measured at 
various points in several sections along the this w could 


‘Since this entire investigation extended over several y years, changes | in 
= to general experimental procedures and apparatus occurred as informa- 
tion was obtained relating to experimental difficulties, and also as the > general 
~ emphasis: of the. study was s altered. The first studies were e made on a 15° 
expansion between pipes of 3-in. and in. diameters. The conical section was 
- constructed of transparent pyralin and the straight pipe of pyralin sheets. 


Studies on this. expansion perfected the experimental techniques, and some 


useful data were secured. However, after this expansion was in use about a 
“year, it warped and had to be | abandoned, thus preventing the collection of a 
‘complete set of data for several rates of flow. Since at that time transparent 
lucite pipe ‘and sheet material had become available, it was decided to com- 
sin, the studies using various e expansions between lucite pipes of 2. 75-i1 in. and 
4.75-in. .. diameters. Time permitted ‘obtaining velocity and turbulence data 
for at least two rates of flow for expansions of 5°, 30°, and 180°. Pressure 
— data were obtained on these expansions and also 0 on a ne between 


oo The general setup used in ie: experimental work consisted ais a constant- 
ba supply tank, calibrated orifice, approach piping, and regulating valve. 
After considerable preliminary work, it was decided to obtain the velocity data 


photographically by use of motion pictures. The technique was as follows: 
"Into th the flowing Ww ater, about 30 ft upstream from the expansion, a mixture of © 


‘ carbon tetrachloride and benzine was injected, the pr proportions being x adjusted. 
so that the mixture had a _ density equal to that of the Water. _ Powdered | 


J anthracene was Wi as suspended i in the mixture to o re it a a milky appearance. ; The 


. and number of w hich could be controlled by adjusting the rate of injection 
~ and by varying the water pressure on the injector. The size of the particles 
considered suitable for this varied between o1 and two millimeters. 
- The water temperature in these experiments varied from 50° to 65° A ben 
sections where v velocity data were ‘desired, the flow with the 
droplets w was illuminated Ww ith a carbon- -are light by focusing an intense beam 


rie 


a double- -slit barrier, which gave a light in. thick. 


A 16 mm, motion- -picture camera was mounted above the conduit and focused 
on the plane of light in the conduit. — The shutter opening and speed of a ; 
- camera were » adjusted for a any ny run, , 8o that for the given ¥ water speed the particles . iN a 
would, on : an average, produce streaks about 4 din. . long when the film was ea 
_ ‘projected to the actual size of the conduit. — A els was inserted next to the 
conduit i uit in the field of v view to give = linear ' dimension. . The field of view 
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Papers” 


only large enough ‘to take in 1 the width of ‘the conduit sit and the | 
i each run, about 15 ft of motion- -picture film was taken, giving about eight | 
hundred individual frames. The time interval for this footage varied with the 

; _ camera | speed, but on the average w was about 1.5 min. For each shutter opening 


and camera speed the camera was calibrated exposure 


4 75" Diameter 


4.75" 


Distance Along Pipe, in Inches oy 

(Computed Q=0.157 Cu Ft per Sec) Gu Ft per See) 


ave 


COCR 
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— Diameter 
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i pictures were analyzed by) projecting th the film to somewhat larger than 


full scale and measuring the length of the streaks i in n the following manner: Ate 
the ‘section of the conduit where velocity data were requir ed, some seven to : 


: nine points were selected across the « conduit. t. At these points on tk the tae 


squares were drawn about 1/25th of the conduit diameter in size. “Streaks st 


falling taken as the velocity at that at that 


or 
with the speed, but was between 1 /20th and 1 /80th a Three 
; were taken each streak—the length | along 
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mean axial velocity at any point. This is designated by the symbol 
‘ Subtracting U from the individual velocity values at various instances es gave u,— ra 
and from these the value of u could be obtained. on 


Q=0.147 Cu Ft per Sec 4.75" Diameter (6) 


— 


(Computed Q=0.148 cu Ft per Sec) |) Q=0.137 Cu Ft per Sec) Cu Ft per Sec) 
- 


4 


4a Diameter 


Square Velocity Fluctuation, in Feet per Second? 
JUARE @ERELOCITY FLUCTUATIONS FOR 7.5° EXPANSION (Q = 0.147 Cu Fr Pur 
- The e length component of the streak perpendicular to the conduit gave the 


"velocity © component, V, and | then, from Lead V), values of the turbulent com- 
n ponent v could be and finally the streak data 
it calculation, at each p oint, of the values of U, and The number of 
to m motion-picture frames varied from 500 to 1, 000. The number of 


Teadings obtained averaged about 300. can be > imagined, t the taking of the | 
data from the pictures and their evaluation were very very time consuming. = 


a ever, since such d data had ‘never been o obtained before, it thought ‘worth 
; 5 while to do the necessary work , since it seemed that some significant information — 


ee The pressure-variation data were obtained by installing wall piezometers — 


along the conduit and | connecting these toa tube bank. Colored water was > 


~ put into ‘0 the top | op of the columns s so as t to make them more visible, and a motion — 


uit minus the streak width (droplet size) gav oa 
rection Of the conduit minus the st - | 
# that A: atinn at that ing ly 4 
7 
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“picture: was taken of. the water levels in all the tubes. % The average value of — 
the pressure was determined by projecting the picture and reading the water- 


column heights from several hundred frames distributed over a time interval : 
E. Mean VELOCITY np VELocITY FLUCTUATION D 


_ several sections of a conduit was first obtained for a flow of 0. 082 cu ft per sec 
through the 15° conical expansion between 3-in. and 5-in. pipes. detailed 

_ Plotting or tabulation of the values of U, u®, and v* for all the tests would require © 
= space; therefore only typical detaited data will be shown for a few 


of the runs. Figs. 1 and 2 contain data obtained for the 7.5° and 180° expan- 
sions for a flow of 0. 147 eu ft per sec. _ From visual observations it appeared 
that on the average the flow was ‘symmetrical; therefore, the actual values 
obtained on one side of the conduit | were also reflected to the other s side to 
assist i in drawing in the a average curve. The actual measured values ar are shown 
gs solid symbols, and the points reflected from the opposite s side are open 
Similar data were obtained for the 30° expansion also for a rate of 

flow of 0.210 cu ft per sec in the 7.5°, 30°, and 180° expansions, and in addition — 


fora flow of 0.082 in the 180° expansion. ll these data will be used in calcu-_ 


. 


_ iin the turbulence gat and i in studying its variation along the conduits. 


apparent that, ata distance of 100 in. beyond the start of ‘the expansion, : 
normal flow ‘conditions had definitely been attained. All the data indicated 
that the u*-values are always larger than the v*-values. The highest u®-values 


v-values occur beyond the sudden expansion. Relative changes from 
section _ to section of u?-values and of »®-values for any expansion a are not com- 
parable; in other words, one of the turbulence components may change dif- 
- ferently” from the other. The highest values of turbulence intensity occur 4 
near the boundaries in the expansion, whereas beyond the expansion the 

turbulence gradually diffuses into the center of the | stream. For the sudden 
expansion (Fig. 2), backflow near the ws alls i is apparent. It w as al also observ ed - 
at the end of the 30° expansion. Only a a slight very nes near the walls. 
was visible for the 15° "expansion, and none was noted for the 7. 7.5° expansion. 
Values of the discharge were computed by. y numerical integration of the 
An mean velocity curves and are shown i in Figs. 1 and 2 2. r In analyzing the large F 
. "amount of data it was necessary to use different people at various times, and, 4 
since considerable personal error can enter into the ‘readings of the streak 
m lengths from the motion- -picture film, it is believed that this explains the dis- = 
crepancy between the actual and alues of In comparing values 


- 4 of ua and e from section to to section, ‘it would be well if these values were cor- 
——-reeted by the square of the ratio of calculated discharge to actual discharge. 


This will be done when the 1 values 1 for the 1 turbulence energy are calculated. — i 
Whether the values of and that were obtained are correct in absolute 4 


is problematical; however, in any case the relative values should be 


| 


as indicated. The actual for normal flow in the straight 
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a (b) 30° Expansion; From 2.75-In. To 4.75-In. Pie 4. 7 ft to. 2. 5 ft per ‘sec. ec. Also, it 
- aN was noted that at this distance the 
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order of wo obtained by H. C. H. Townend" 


tended to be large than that of in Mr. Townend’s studies. He measured 
the w-component and found it to 
1.—Maximum VALUES of the order of the v-component, 
(perpendicular to the 
of | Q, Couns all) tended toward zero, whereas 
the w- -component ¢ did not. 
gs = 
2) 0. 0.082 Significant items to be noted from 
eins aaa ——————_ Fig. 2 are the manner in which the 
(a) 7.5° EXPANSION; From 2. 75-IN. To 4.75-IN. Pire ‘mean center ¥ velocity oni 
0008 .. the turbulence gets into the 


0.094 
Se : ——_ stream i in the 180° expansion. 


esos 


the exp expansion, the center velocity 
reduced from a maximum of 


4 turbulence created at the boundaries 
the expansion already per- 
meated into the center of the stream. 

foregoing reduction in center 
aes |. __ velocity and permeation of the tur- 


ea 
() 1 180° EXPANSION; Frou 2. 75-IN To 4, 75-In. Pree _ bulence d do not. occur, for the 180° ex- 
pansion, until a ¢ distance of 1.5 ft 
| 0.110 3 110 beyond the expansion (see Fig. 2). 
--+ 
0.284 Evidently” the turbulence created in 
0.150 . _ the 30° expansion is of a larger scale 
has greater diffusive power, and 
thereby penetrates and slows down 
(d):15° Expansion; From 3.00-In. To SIN. ‘Pie 
the center portion of the stream in 


Of the relative intensity of the tur- 

—- in the conduits” with the 

location of the section i is given as dis- 

= aol in inches from the beginning of the expansion. various expansions, Table 1 1 contains 

_ The values given at zero distance are in the small values of the ratio of maximum value 
a - conduit, although not necessarily just at the end of the soak 


small conduit. Ga Vu u? to the center velocity i in the 


= cross section, Data shown in this 


__ table are for all the measurements made | in these experiments. It is noted 
Aza — the maximum value of this ratio occurs” for the 30° expansion and is 


a about 0.50 for the flow of 0.210 cu ft per sec. This ratio seems to have 


order to indicate thes variation 


12“*Statistical Measurements of Turbulence in the Flow of Air Through a Pipe,” H.C. 
» Royal Soc. of London, Vol. aA, p. 180. 
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- about the same value in the 2.75-in. pipe and in the 4. 75-in. pipe, when | 


normal fi flow is again developed. _ It seems that it took a distance of about 125 


‘in. or about twenty-five diameters from the start of the expansions for this 
ratio. to attain a constant value in the large pipe. » For any expansion and at 


any , given en section this r: ratio, in general, has a smaller value for the —. 


pa An examination of ‘the typical detailed data shown it in Figs. 1 and 2 wii 


al values in Table 1 offers a fairly complete picture of the variation of the turbu- ad 


lence in conduit expansions. 


, in Foot-Pounds per Second 


ounds per 


‘oun 


0.1 
8 
J 
if 


Foot-P 


in 
E,, in Foot-Pounds per Second 


lues of 
E 


Values o 


Values of 


Peek from Start of “Expansion. tal 


OF Flow OF FoR 7.5° EXPANSION Gin 


To indicate the changes occurring in the kinetic energy of 1 the ‘mean f flow 


and of the turbulence e energy along the conduit for the several expansions, ‘to 
total value of these energies: was calculated by numerical integration of data 


such a as those shown i ‘in ‘Figs. ‘Land 2. kinetic energy of the. mean 
DA. 
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ction, y is the the center. the of the mean 
velocity was uniform, EZ» could be e expressed as - t's , in w vhich U m aie 


| velocity in a cross section. For normal turbulent flow in a 


T= 


© Q=0.210 Cu Ft p per Sec | 
—s Q=0. 147 Cu Ft per Sec 
Q= =0. 082 Cu Ft per Sec 


n Foot-Pounds per Second 


9 
” 
a 
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Et 
Em’ 


Values of 


4. ENERGY OF MEAN FLOW AMD ¢ oF | FOR 180° EXPANSION 


the mean velocity i in the cross ever, for’ velocity 
distribution in and immediately, beyond expansions, the correct vé of En 


“may , be several times the value obtained by using ‘the mean velocity in the 

section. kinetic « energy of the turbulence per of fluid 


The total kinetic energy of the xy of the turbulence i in the Ww vater ater flowing per er unit time 


flow in trom the following expression: pas' 
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past any section is then given by: — 


near the w alls of the conduit. The values of E, were 
lant by numerical integr ation of data such as those shown i in Fi ‘Figs. 1 land 2. - 

_ The ratio of calculated to actual discharge ata any section was assumed to 7 
‘indicate the error in the streak length or in the length o of velocity : scale used. 7 


were therefore corrected by the cube TABLE 2.—Kinetic ENERGY oF 


of this ratio. In Figs. 3 and 4and ‘Mean FL 


in ‘Table 2 are shown | the values of | = 
Location | Krneric ENnerecy, Foor- 
Ratio, 


Et, and the ratio E, t/Em for all “of | Pounps Per Szconp 
section? |__ 


the data obtained in these "studies. ‘or 
( 


mentioned regarding these data: The 
(a) 15° Expansion From 3-In. To 5-In. Pipe; 
maximum value of occurs right at 
$ the end of the 7.5° expansion; for 
the 15° expansion it occurs 0067 a 0.030 
3 beyond; and, for the 30°. and 180° 0110 | 0.088 
expansions, the maximum value of 0. 
110 
a E, is about 0.60 ft ‘beyond the ex- 
pansion. The maximum value 0 of the (b) 30°. EXPANSION FROM 2.75-IN. TO 4. 75-IN. Pree; 
4 ratios E,/Em always occurs beyond Q@ = 0.147 Cu Pr Pern Seo 
tance ) ~ 12° 0.078 
were 1.4 ft for 180°, 0.65 ft for 30°, 0.221 | 
and 0.5( 0.50 ft, for 7.5 5°. tt turbulence 
energy appears to have ceased chang- 


30° E 4.75-In. Pree; 4 
= 0.210 Cu Fr Per 
ing in the 4.75-in. ‘conduit: at a dis- 
tance of 9 ft or about twenty- three 7.56 0.014 


3.85 0.099 
diameters from the start of the 7 oss | 140 | | 
0.765. 0.248 0.325 


+ 


— 


is fairly small beyond about 6.75 | 0570 | 0.018 | 


ten diameters from the : start of eg « The location of the section is given as the dis- 
expansions. tance, in feet, from the beginning of the expansion. 
Th Values given at zero distance are in the small conduit. 
approximate maximum val- 
pansions te tested are as follows: 7.5°, 0. 17; 15°, 0. 0. 28; 30 ° and 180°, 0.50. The. 
values appear to be somewhat higher for the lower rates of flow. The actual 
values of E, for the 30° expansion are slightly greater than those for the 180° 
expansion. The rate of decrease of E, beyond the expansion is 
most rapid for | the 180° expansion, indicating that the turbulence is probably — 
finer grained for the 180° than for the 30° expansion and thus dissipates faster. 


‘This was mentioned i in connection with another item. The fact 


_— connection with another item. The fae 
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‘great. of the at a point in the fluid corresponds some- 

what to the magnitude of the friction or shear stress. If this intensity is high, 
7 > if it is quite localized, the energy loss due to friction is certain to be large. 
< It was n — oted that ¢ the maximum intensity of the turbulence for the 180° expan- 


‘Viglen. wilt i is about 1.3 times greater than 


TABLE Vantatton that for the 30° expansion for any 
ALONG ConpulIts, IN FEET OF 


4 (Expansion from 2.75-In 


. Pipe to 4.75-In. Pi given rate of flow. 
In, Pipe) The aforementioned items seem 


ot tob be the principal information 
| can be « obtained from an examina- 
15° E tion of the E,,-curves and E,-curves. 
These quantitative data give some 
information relating to the changes 
0.051 


0. the mean kinetic and ‘turbulence 
energy” that take pla vce in the flow 


To complete the picture of the 


‘conversion "process for flow 

through expanding co! conduits is nee- 

0.860 0250 essary to have data on the variation 


‘sured for the 7.5°, 15°, 30°, and 180° 


location a the section is given es the: dis- expansions: in the manner described 

tance in = from the of the in Section ‘The pressure varis- 
shown in in Fig. and Table 3. 
The v very ; appreciable drop in the wall pr pressure ‘at the start of the larger 


duit for the sudden angen, as shown in in Fig. 5(6), 3 is quite characteristi 


on 


the stort of the expansions w whieh: seem to increase with in 
: angle; thus, for the 7.5° expansion this distance i is about 4.0 ft, for 1S it is 
3.3 ft, for 30° it is 2.8 ft, and for 180° it is 2.5 ft. te 
‘In order to indicate the change in the total energy, exclusive of the turbu- 
lence energy, 2 along the conduits, the quantity (E> + Em) W was calculated and 
is plotted i in Fig. 6. The value of E, was obtained by taking the values of 
4 pressure head, h, and ‘multiplying by 62.4Q. It is worth noting that - the 
‘aoe that occurs in this total energy in the 180° expansion is much more 
; > a gradual than that in the 30° expansion. In other ‘words, as was mentioned 
m - previously, the energy conversion occurs more gradually i in the sudden expat 
i m sion than in the 30° expansion due to the fact that the relatively small-scale 
_ I meromaes that is generated at the edges of the high velocity stream i is not # 
effective in “ ‘eating into” the stream and slowing it down. . The diffusive 


of the turbulence depends on both its intensity and seale and ev’ and 


a _ that the 30° expansion appears to have as much total turbulence energy as the a 
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turbulence energy are minor in comparison to the changes in (Em + E,). 
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Q=0.147 CuFtperSec 


.75"' Diameter 
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6.—V ARIATION OF Kinetic Pius POTENTIAL ENERGY OF MEAN FLOW © 


= the turbulence generated in the m more gradual 30° expansion is of es “_ 
and thus gets into the main stream more quickly 
Bis fa If the values of E; are added to » (Em + E,) it will be seen that changes in 


i = other 1 words, the total « energy loss 0 occurring is several times the maximus 
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then, is the exact réle play ed by the turbulence? It appears that the: major 


portion of the kinetic | energy of mean flow is lost directly due to the a 
high, and quite variable, shear stresses that are set up in the fluid on the ~ 
: side of the high- velocity jet. The turbulence created is a kind of | by-product _ 


_ which, of of course, uses up some of the kinetic energy of the mean flow. 4 _ A; 
re a by- product of the ‘major energy-loss process, the turbulence, — 


the less, has an important influence on slowing down the high-velocity water 


and i in shaping the velocity distribution | curve. 


‘From m the values of (Em E,) it Was possible to ¢ calculate the actual loss. 


3 of a. This was done by taking the difference between the values of 
+ (Em - E,) at the start of the expansion and at about twenty to thirty diameters 


the expansion and subtracting from the normal friction in the 
conduit between the two points. In the conical expansions the friction loss 


was calculated by summing the loss for increments along the expansion. - The 
values of Br energy loss, expressed in feet of w rater, are shown in n Table 4. 


"Ako, the efficiency of the energy conversion is given, which is nes as: ane 


” 
‘in which | (he 1) is the 1 rise in pressure head (after correction for. friction loss), 


and the denominator is the difference in the mean kinetic energy. manuel 
_ TABLE 4. ~—Enanor Loss, h, AND EFFICIENCY, n, OF VARIOUS EXPANSIONS 


PROM 2.75-IN. PIPE TO 4.75- In. 


30° E 


XPANSION 1 15° ExPaNsION 
Flow (cu ft} 


XPANSION 180° 4 
sec) 


7 
0.027 
| & | (0.296 


e energy as vans in Table 4, do appear to ) vary con- 
 sistently with rate of flow . Checking the loss figures — given. for the 180° 
expansion: with the values obtained from the | Borda formula indicates t that f : 
- the higher rates of flow the figures given are appreciably larger. ‘The r reas 
+t these apparent discrepancies is that the kinetie energy of the mean flow i in 
‘ ‘the small pipe is considerably greater that that which Ww ould be obtained using 
the mean velocity i in the cross section. #&F or instance, for the flow of 0.210 cu : a. 
it per sec, » the ratio of the actual to approximate total kinetic energy in the- 
2.75-in. pipe is 1. A2; for the flow of 0.147 cu ft per sec, ¢, it is 1.11; and, for 0. 082 a 
cu it per sec, it is 1.02. . The \ variation in this ratio accounts: for the fact that: 
the loss: does not vary as the square of the rate of flow. For 0.082 cu ft per 
sec, the loss in the sudden expansion checks that obtained from the Borda 
equation very well, and the ‘Teason is apparent. It is thus obvious that the 
energy loss. in an expansion will depend on the velocity distribution at the 
inlet side, and this d distribution must be taken into account if the : 
is to be correctly evaluated. . Itis noted that the ‘efficiency | of energy conver- 
sio on in —_ 30° expansion is not much better than in the 180° expansion. = 
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oo Analysis” of experimental data obtained in this investigation leads to the 


ef follow ing conclusions relating to the energy conversion and dissipation processes _ 


23 The largest. energy occurring in the flow through any expansion 


are those due to the a in velocity and change in shape of the mean 
velocity distribution curve, 


2 The maximum total turbulenee energy i is a small part of the total energy 


change taking place. ‘However, at any particular section the turbulence: energy 

aaa may bea considerable part of of the total kinetic energy. | For the 30° expansion 
and 180° expansion, the maximum value, of the ratio of turbulence energy to 

kinetic energy, Em, rises to 0.50. ‘For the 7.5° expansion it is only 0.17. 


~ 
The maximum ratio the ‘root-mean-square value of the axial 


7 velocity fluctuation, u?, and the mean center velocity occurs in the 30° expan- 


i and is 0.50. The transverse velocity fluctuation i is always appreciably 


lessthan the 


ut 4. The rate of loss of nee occurs more rapidly in the 30° expansion than 
e = tlt ‘The principal loss of energy occurs at the outside of the high- ~velocity 
_ stream as it passes into the expansion due to the very high shear developed. 
: This loss is represented by a direct conversion into heat of the mean kinetic | 
- - energy of the flow. The turbulence produced is but a by- product of of this energy 
conversion 1 process. Of | course, the turbulence created also | represents lost 


; energy, since eventually i it is into | heat. 


"process are influenced considerably “ng the distribution of the mean velocity 

at the inlet to the expansion. — High ratios of center velocity to the mean 

» velocity in the cross section produce poor efficiencies and large energy losses. 


Since the experimental work reported herein eaten. over a period of six 


years, it is natural that many individuals should have contributed their efforts. 
a project w as sponsored by the Hydraulic Research Committee of ‘the: 
i. Society, and financial aid was provided by the Engineering Foundation. — The 


os work was done in the hydraulics laboratory of the Iowa Institute of Hydraulic 


Mavis, M. Am. Soe. C. E., formerly head of the of 

# oe = Mechanics and Hydraulics at the University of Iowa, at Iowa City, is Te | 

_ sponsible for starting this investigation, and for his advice, suggestions, ‘and 

ae Re ~ encouragement the writer is very grateful. Graduate students who assisted | 

at various times with the laboratory and analytical work were: E. R. Van 

ee and C. L. Pien, , Assoc. : Members, Am. Soe. . C. E., A. R. Luecker, 


J. M. Robertson, and Claude Lomax, Jr., Juniors, Am. Soc. C. E., and W. M. 
achter. The National Youth ikdininistration provided a ‘considerable nun- 
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an- n empirical golumn formula should be one that expresses theoretically 
bly possible relations. The prime requisites of | any column formula or formulas are 
the average stress be a continuous, monotonically "decreasing function of 
ol Ss the ratio of slenderness and that it approach the ‘Euler > value as the ratio | of 
city "slenderness becomes infinite. ‘Another ‘desirable requisite is that the first 
ed. derivative of the average stress with res pect to the ratio of slenderness be con- 
D 

etic tinuous. These requirements impose ‘certain restrictions on the constants in 

Tey il f d th tricti id od fc abe 
empirica ormu as, and these restrictions are consi or a num er of we 

7 known and less well-known column formulas. . 
eu 


rsion ra, Column formulas area convenience, not a necessity, for r relating the strength 


of a-column to its dimensions. is. It we would be quite possible to o plot the results of » 
_ tate, , determine a curve from the plot, and select design values from this curve — 
or 0 tildes of values read from the curve. W hen a formula i is s used, 
“ever, as it usually i is, there i is some virtue in having it of correct form, if only 


all ieee! the designer may have to extrapolate beyond the values on which the 
yf six 


is based. for example, at large ratios | of slenderness a formula 
- covering all ratios of slenderness leads to values of average stress that are known © 
tobe impossible, the formula will lead at less large ratios of slenderness to values — , 
that are likely to be impossible. — Therefore it is desirable not to preclude 
extrapolation. Is it not at least slightly better to know that extrapolated values Bs 


“may be right than to know that they probably are wrong? 
j on 1 elastic theory for the design of centrally loaded columns in the plastic range. 


‘There have been numerous and elaborate attempts to devise formulas based ~ 


- Before. the development of the double-modulus theory toward the end of the ps 
century these attempts were ‘consonant with good engineering practice. 


Today, however, it should be realized how futile such attempts must be. _ I It 
aa _ Nors.—Written comments are invited for immediate publication; ali insure publication the last al 


., National Bureau of Standards, Washington, D. Cc. att 
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— would be possible t to design wei me: ins of the double- modulus ‘theory, given the 

4 compressive stress-strain curve of the material; but the best testing technique 
- results in wider discrepancies between theoretical and experimental values in 
“ae plastic range than i in the elastic range, and most engineers would undoubt- 


2 Be prefer to use an empirical formula based on tests. An empirical formula, 
how ever, one that theoretic: ally possible Besides | 


structure, implicitly or explidtly, and one or more constants related to 
material. ould be very hard to justify n more than tw vo constants: of the 

AB to the for m, the } prime requisites. of : any y column formula or formulas are 
that the average stress be a continuous, monotonically ‘decreasing function “ 
‘i ratio of slenderness. and that it approach the Euler value as the ratio of 
slenderness becomes infinite. Another desirable requisite is that first, 


derivative of the average stress with respect to the ratio of slenderness be- 

continuous, 

It is proposed to present a number of formulas i in n their x most general forms 


satisfy the requirements of the preceding paragraph. or comparative 
_ purposes the formulas will all be written in terms of the two. non- dimensional 
variables, (the ‘reduced ratio of and (the reduced 


(1b) 


_ as w well as in ial of the eus customary variables— the average e stress P/ A, and te 


“cessive points of inflection. For supported = tee fixed ends, 


3; for one end free and ‘tee other end fixed, k = ete. In Eas. 1, Fis the 


‘modulus of elasticity of the material, and 8 is an arbitrary constant of the 


_ dimensions of stress. . In analyzing column-test data, S should be taken as a 
at compressive secant yield strength. . (The secant yield strength i is the stress at 


pe iy 


a intersection with the stress-strain curve of aline through the origin having 
slope m E, where < m 1. ) SUGT? ped over 


The most. general forms of some well- -know n and other, less well- -known 
= formulas will now be given. The dimensionless constants 8, 7, and n, and the 
constants B and C that « appear are arbitrary, subject to the conditions stated in 
case. They are constants related to ‘the material. The load, P, is 1 to be 
the maximum load, a ‘and after assigning values to the 


to be taken as Ww working lo load a factor of safety. le 
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a ‘The only formula co: containing but fends constant that will be con- 
is the Rankine formula. . This formula is used extensively. 
justification for it is that at low and at high ratios of slenderness it fits test 

- results well, and at intermediate ratios of slenderness, w here the scatter in 
~ test results is relatively large, the formula gives | lower values of the average 7 
- stress than the mean of the average stresses found in tests of most ‘materials. 


neta cee applies for all values of ratio of slenderness. 
I. Toe RANKINE FORMULA 


ormula is expr 


= 0, Eqs. 3 3 reduce to the Euler formula, 


= 


ttle 


pay 


of ‘the resulting line Gt ne results!) through the 
1. Alternativ ely, since 


8 may be taken as the average of the v values of — — 2. 


® 


There are several formulas containing two arbitrary eat ati one of 
F them, and often more than one, can be made to fit almost any column-test data 5 
f - that it may be desired to fit. a First will be considered formulas covering the 


range of ratios of slenderness. 

. THE SECANT 


secant formula | is 


1e 
ne 
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he | 
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COLUMN FORMULAS 
see? 


tone 
Th 


e constants B and may be determined 1 by plotting sec against 


Bo 


pres = and C are the slope and intercept, ee of what may turn out io 


Ill. Formuna 


A column duet ‘to M. M. G. Aarflot? is: 


‘When =B 0, or C=1, Eqs. 9 reduce to Euler’s 7 
_ The constants B and Cc may | be determined from a series 0 of observed values 
~ Mo against 
+00 


and 2C and B a- - are the the slope and intercept, the straight 
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2A Mer Column Formula,” by ‘Martin G.. Bulletin, Am. Sen. of Swedish Engrs., Vol. 26, 


 §**Keskeisesti Sauvan Nurjahtamiseste Ylapuolella,” by A Arvo A. 
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‘When 6 = B 
and n= ‘they to Rankine’s formula. 
Tti is easy to determine Ww whether the formula i is a good fit to test data and to 


and it is only necessary to plot observed values of l—NW¢0 against o on log- 
-arithmie paper. If the plotted points fall on a straight line, the constants may 


be from the and intercepts of this line. 


Provost BY ‘THE Warren 
“a 


1 by | tee Ww Titer is: 


+i 
J rmula, and, when 


p= y= B= 0, Eqs. 13 reduce to Euler’ 
8B > 0, _B > 0, and =C= 0, they: reduce to Rankine’ s formula. 


The 


Ao 
‘ and, if the resulting plot i is approximately a straight line, y and B are its ~<a 


0, 


— 
. COLUMN FORMULAS 1569 
— 
— 
= 
ay 
a= * ¥ a 
A formula propose =} 
— 
q 
may be ermined by plotting — A“ against A., 
2 
areierence to them if any reader knows of one.) 
Another suggestion by the writer is: a 


Pa Papers pers 


4 
a of the axis of to o= = line 
= 1 betw een Xo = 0 and A = 1, the reduced Euler curve ‘beyond 
As in the case of Eqs. ll itis easy sy to determine whether Eq. 15a. represents 
a a good fi fit to test and to determine the constants if it does, for 


Eqs. 15 are derived by the the double-modulus theory fr from a stress-strain curve. 


p 


in} which Ka nia 
A third suggestion by the writer is: 


Ww hen B = 18 Euler’s formula. ‘Like Eq. 15a, as n as n ©, 
S: | Eq. 18a 8a approaches. a broken curve consisting of the axis of o down to o 1, the 


"Tine ¢= 1 between and and the reduced Euler curve beyond 


Asi in the case of 11 the Kgs. 18 may 


No 


he next two as (Eqs. 20 and an two ‘constants 
each, but the formulas are intended for use only above the limiting average 


Parazouic-T TYPE OR 


owl g 


nB 
n+2 


_  4**Column Curves and Stress-strain Diagrams,” by William 1 Journal of Research, National 
a 


‘Bureau of Standards, Vol. 9, October, 1932, 


stresses indicated. | ‘The curves they represent are tangent to the Euler curve, 
for stresses below the point of — the Euler formula, 4, should 
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or 


W hen n = = 1, Bes. 20 become the. straight- line formula; and when n = 2, they 
become the Johnson parabolic formula, 


A formula due to H. 


= 


Gay 


“ant any fort formula can b be plotted t to as s large a scale as desired and values read —¥ 


the graph, or it m may be tabulated with increments in | the —— as — 
as 


The column formulas in Sections II to V, inclusive , apply for all values of the 
Tatio of slenderness, and one or all of them can almost be guaranteed to fit 
adequately any reasonable column-test data it may be desired to fit, = 8 


till fraga in om striifvors knickning,” by ‘Kreiger, Teknisk Tidskrift, och Vatten- 


value. of the po of hayend which the Euler formula must be 


When y = C = Eas. 20 reduce to the Johnson parabola. 
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AERODYNAMIC (STABILITY 


senting the stiffness om suspension bridges by the summation of exp expressions that 
characterize | the ‘stiffness of of the e component parts—cable e and girder—of the 
center span. — ~The cable is assumed to be a string under tension, rigidly held at 7 
its ends; the girder, to be a simple beam of uniform section. « ‘The accuracy of 
this approximation when determining the frequencies of a practically” 
stiffened bridge is well indicated | by the agreement between the calculated and > 
observed values for the Tacoma Narrows Bridge. However, the equations 
- given for the condition of three spans coupled do not define the values of m * 
corresponding to values of nin Table 1. For bridges with stiffening girders of 
small rigidity, the exceedingly simple 17, which was presented by ‘Joseph 
Melan in 1906,” gives useful | results. frequencies of the Tacoma Narrows 
ee differed by only 6% up to the observed n = =9, from those obtained by 


‘The K- we of Eq. 6 are characteristic of deformation when the algebraic 

_ sum of the areas of the two segments below and above the initial condition of _ 
“the cable equals zero, or when the deformation occurs without change in H, bs 
without stretch of cable, without mov ements ‘of the towers, and without in- 
volving the side spans. For these reasons, Eqs. 16 and 22 cannot reliably be 
used for estimating live-load deflections and rotations which, , besides being - et 
influenced by the enumerated changes, are also affected by the ratio of live a 
load to dead load. close an in Example 2 and Example 4 between 
the results obtained from these equations and those obtained from the deflection _ 


theory i is accidental as shown in Table wat 


_ _Nors.—This paper by D. B. Steinman was published in November, 1943, Proceedings. Discussion — 
on this paper has appeared in pee my as follows: February, 1944, by Edward Adams Richardson; 
March, 1944, by R. i. Bernhard; May, 1944, if M. Robertson, and C. H. Gronquist; June, 1944, _ 
Farquharson; and September, 1944, by H. Reissner, Birdsall, Hardy Cross, 
Smith, Hermann riedlaender, and George § . Vincent. 


| 
 ™Cons. Engr, New York,N.Y. 
Received by the Secretary August 31,1944. $= 


_ % “Handbuch der Ingenieurwissenschaften,” by Joseph Melan, Vol. 5, 3d Ed., 1906 
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‘The: application of E ugs. 16 and 2 22 to: a ges in addition to 10 those 
in ‘Table 12 resulted in deviations, up to a maximum of 48 %, from the values 


obtained by deflection theory devi ations also rove: that th 

~ author’s physical | inter values (as the loading from tower 


to midspan in n pounds per square foot which 


ATION OF Sea. 6: AND 22 TO _Varrous BripGes 


(a) Eq. 16 Comparep (b) Ea. 22 ComMPARED WITH 


by: 
“Meas 
| “Gon | erease |B crease 


= 
(116 
Golden Gate ate... 26 2,000 


A 
George | 100| 0. “4: 


++ . 


oS 
£2 


tl ++ 


0. 0. 


Eg, 22, T = 10,000 pert 


“deflection of stiffened suspension bridges) is eopennest. The percentages given 
ré f th rpr n. 


= The basic cable term of Eq. 6 5 , however, can be used for the 


tials. of qt quarter- point deflections of “the unstiffened cable, due to the 
tee mentioned, if it is considered to be inextensible and fixed at the towers. 


For this ¢ condition, 70.2 =f & .% or p = 100, considering the 


_ prevailing dead loads, — is between een 12 and 150 and, even in the extreme case, 


average value, 7 70.28 = 1.5 isa close approximation 2%) of the 
_ To connect the stiffness and aerodynamic behavior of suspension bridges, 


. author evaluates Eq. 6 for tw enty bridges, divides these values by the 


width, and lists the bridges in order of the magnitude of these values (Table 2). 

a From | experience the author states that stability requires these K per foot values 

a to be not less than 15 for bridges stiffened by plate girders and not less than 6 fo: 


-truss-stiffened bridges. Tn this n manner, an empirical connection is is establishe: 


between the =—-values and the oscillations of suspension bridges 
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ay The use of the basic — of Eq. 6 to express the e stiffness of of a suspension 

bridge able and its stiffening girder has been practicable in suspension bridge 


he design. The concept of extending these terms by the corresponding first 
- ened coefficients and dividing their s sum by the width of the bridge is new; 7 
the significance asribed to the result 
int 
nf i, for the --ratio of 100, and for the minimum values for K and R can best 
= | be discussed by the application of these rules to several suspension bridges 
| nating the span range of structures actually | “built. 
_ Table 13 demonstrates that, for r spans longer than 1,000 ft, the requirements 


TABLE 18. —APPLICATION 0 oF THE AuTHOR’s MINIMUM REQUIREMENTS FOR 


weer “Truss 520,000 | 2% 

Tacoma Narrows. gt : 
» the result i in plate girders of impractical eal dimensions. For the Golden Gate Bridge, > 
von. © trusses would be required which would nearly double the cost of the structure. 
- This bridge has been exposed for s seven _ to the most unfavorable wind - 
wail and safe 
winds of from | 60 per hr to 78 per damped out when 4 
p the wind velocity decreased. Therefore, the Golden Gate Bridge may be 
assumed to. represent the limit of satisfactory dynamic behavior, contradicting 
the general v alidity of the author’s dynamic stability criteria. 
idges, ‘The e requirement that Plate-girder stiffened bridges be 150% stiffer th an 
9). trus-stiffened their use. ‘The Tequirement that that = = 100 
n 6 for has the same effect because it is destructive t to advantageous | Ftatios. s. N or 
ni could plate girders be used for spans longer than 1,000 ft if — = 100 and an 


eh had to be maintained |. Nevertheless, the author concentrates on the 
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investigation of this section whieh, i in view of experience and his recommenda. 


' The author’ 8 s analysis of 1 the aerodynamic e excitation of vertical a1 and sain 
-oaslations is based on the characteristic drag, lift, and torque curves obtained 
ons section models of the Tacoma Narrows B Bridge. By idealizing the laboratory 

“results, the author obtains a family of lift graphs corrected for drag (Fig. 5). 


The introduction of the K-values into the frequency, as well as into the energy, 


~ relations of vertical and torsional oscillations has the effect of establishing the 


= between s’ (the slope of the corrected lift graph) or % (the slope of 
the torque graph)—both for a = 0—K, V, B, and the aerodynamic amplifica- 
tion ratio 6. The arbitrary function, 8, of the phase’ difference across the 

Sn section is also introduced for vertical oscillations to explain the behavior of 
‘the Tacoma Narrows Bridge in the low-velocity Tange. . The values of é,, 

x derived from the moment of inertia of the girders equated with 5, ‘result i in 
Eqs. 44 to 46 for the critical velocities of vertical oscillations and in Eqs. 62 

_ to 64 for the critical velocities of oscillations, 
mx. The author | proves the validity of these equations - examples w which give 
values corresponding to the observations on the Tacoma Narrows Bridge. 
Applying the equations to other bridges, however, the results are not con 


vincing the Thousand. Islands, Bronx. -Whitestone, and Golden | Gate 


> bridges ( ——— = 0.139 q Eq. : 45 b with n=9 yields 48, 85, and 70 miles per 
: hr, and Eq. 46b with n = 2 yields §, 1.7, and 2 miles j per hr, respectively r. The 
Tacoma Narrows Bridge subjected to tos a wind velocity of 42 miles per hr changed 
from a ‘a vertical oscillation with n = 9 to a violent torsional oscillation with 
n= = 2 The maximum critical velocities determined by the author’s equations 
7 were experienced by the enumerated bridges without any such effects. — 7 For 
this | reason, , the practical meaning of Eqs. 45 and 46 is not apparent. __ = 
ate 62, 63, and 64, which show i increasing safety | from torsional ane 


“For more than five years | the Bronx- Whitestone B Bride 
, but at no wind 
bint oscillations occur. . Using Eq. 62a, the. critical velocity 


_ of torsional oscillations of this bridge i is 18.2 miles per he: and, using Eq. - 64a, 
18.5 miles per hr. that everything is unchanged except 


The evaluation of observations on the Tacoma Narrows Bridge and on its 


ae section models failed to give . generally * valid relations even for the solid H-section. 


"el The author’s analysis is more interesting from the point of view of procedure 


than of practical applicability | to design. fact i is apparent from “Part 
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ted Bry. 
icati results of the presen 
i measures do not follow from the application of the 


“y analysis. The c condition 5 6, = 0 is a basic relation in the critical velocity a 

a equations (see heading, “Part 2.—A Aerodynamic Excitation of Vertical Oscilla- 14 _ 
nal tions: Structural Damping’’) w here 5, is essentially a function of R, the pro- i. 


portion of K due to I (the moment of inertia « of the girder sections). “Tn build-_ a ar’ 
ory ff ing up the resistance of the structure to aerodynamic action with stay ropes, — _ 
5). the author assigns a new ‘meaning to 4, by explaining (see heading, “Part 4 
‘gy, § Methods of | Securing Rigidity and d Aerodynamic Stability: ‘Use o of Longitudinal cs 

the Stays”’) the mechanical interference of stays with oscillations; and then writes 

e of (but without derivation) Eq. 65 for the he amount that K is augmented by the 
ica- B stays. Eq. 65 is similar to the cable term, and not to the girder term, of Eq. 6, 
the § § being a function of a unit stress and the initial stay Tope sag. ” ‘Hanes, the E he 
r of addition of this quantity to the girder term in the equation for K, to show a 

desired degree | of increase in is not warranted by the author’s ‘explanations 
t in Bor by the fact that pre- stressed stay ropes stiffen a suspension bridge to a_ 
. 62 B certain extent. _ With the recommended minimum K = = 300, the cost - of the a 
Tacoma Bridge, stiffened by plate girders, would be more 
give doubled (see Table 13); whereas by 2 application of 16 sq in. of wire ropes strand, 
dge. in the form of longitudinal diagonal stays, the K- value of the bridge that 1 failed 


con: | would be raised to 332 and R would be increased to 0.816. In other words, a 
Sate structural damping well above that of the Williamsburg Bridge would be 


2 secured at insignificant expense. © _ This result, which i is of great economic im- 
tess portance, demands a more thorough proof than the explanations i in the pa paper. 7 


‘The Experience with stays reference to the Brooklyn Bridge, which was 
nged stiffened by s six trusses, without —“ having values of K = 190 and R=0. 22, 

tions author’ s reference to the of transverse diagonal stays 


“For top struts, with the out-of-phase movement of the cables, does not prove that page 


 -«BEq. 66a represents quantitatively the dynamical effectiveness of this arrange- ths! 
tions, | ment (Fig. 9). _ The makeshift appearance of this detail (which has been used ie ; 
avior ff 00 ore- -transport suspension bridges) does not befit the permanent and monu- _ 


mental character of a vehicular bridge structure. 
wind The e author’ s statement that the aerodynamic forces causing the vertical and 
locity torsional oscillations of the H-section can be eliminated or reduced to. negligible 
. 64a, Proportions by leaving open spaces near the girders or by open-floor —— a 
the tion (see heading, “Part 4. —Methods of Securing Rigidity and Aerodynamic 
Stability Use of Open-Floor Construction’ is contradicted by experience and 


vind. tunnel tests. The Bronx-Whitestone Bridge was constructed in the 


xe the All steel, including -reinforeing grid with floor was 


erected and sidewalk areas open, leaving two lateral v vents with 

i total width of 13 ft ~ At this stage, — 1. 

cedure (2) Roadway concrete was poured ond sidewalk areas were still open. — At 

1€ 
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ALOG ‘ON A AERODYNAMIC STABILITY Discussions: 
@ Sidewalk concrete was poured and the vents in stage were 
7 In none of these stages, nor during the following five years, did the bridge eins a. 
torsional oscillations. largest vertical oscillations occurred, however, 
—_— stage (1) w hen, according to Eq. . 60, the aerodynamic lift was reduced : 
60%, and, in addition, the multiplier Bw as reduced 40%. + 
- Homer M. Hadley, Assoc. M. Am. Soc. C. E., = , describing e experiments with P 
the 100-ft scale | models of the Tacoma Narrows Bridge, states:77 
7 pep. “The steelwork of the original deck w as reproduced to ay it. the 
marginal stiffening girders, floor beams, stringers and wind trusses—but_ he 
the floor : slab and sidewalks were omitted. This model lacked 20 canal { lis 
Bt of the stabilizing eight of the prototype PA — = 0.67 | * ** It was” the: 
tested with and without midspan stays. Without stays *** it 
quite closely reproduced ‘vertical undulations of the prototype 


ae nor opening the horizontal member of the solid H-section 
iio the aerodynamic behavior of the: section as Eqs. 43 and 60 indicate; and 


a 


the dimensional ratio — cannot be as a measure of sensitiveness to 


~The structural arrangements, by the author as methods for 


securing aerodynamic stability, have been used on actual bridges and in pro 


posed designs for many years. It is not proper to call the H-section “conven. Mert 

tional.” _ Many existing suspension bri bridges have shallow, deck plate- -girder, +. 

suspended structures, and it is purely that bridges using 
this form have not been built. Such suspended s structures were proposed EI 

for the -3,000- ft ‘Wilmington Bridge, in ‘Delaware, « and also for the | 4,700-ft as fo 

Battery-Brooklyn Bridge in New York, | B.. ¥. , in 1939, because th this design =" 

satisfied the requirements of those bridges better than the H-section. — All stay and 4 
systems (center, longitudinal, and transverse, Fi igs. 8 and 9) discussed by the az 

author were suggested by others, or were tested in practice long ago, and then BY 

now they have proved themselves mostly inefficient. In 1926, open n side- be 

Walk ‘Spaces were used in the design of the proposed bridge over the First it 

- Narrows, Vancouver, B.C., Canada, by Henry C. Tammen, M. Am. Soc. C. E. te 

: Also, the 1 ,050- ft Otto Beit Bridge, built in 1938, over the Zambezi River in in 


1odesia, Africa, is with vents of width 2? ~ 6 ft. Outside 


i - sidewalks have been used on bridges too numerous to list. Ni Static or structural 
considerations suggested all these arrangements. Their effectiveness for 


securing aerodynamic stability is is not experience or 


wake relations from the Tigidity of a a basic str ucture 
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-empirical relations, the 
properties will be expressed by the cable and girder stiffness 


parameters: and =F ET’ , and, by their ratios, the dimensionless 

= — = and = arameter 71, refers to the 


‘jeden structure. The writer utilized these parameters by « equating the ms-values_ 

of the Golden Gate and Bronx- ‘Whitestone bridges, thus deriving the value of of 
[required to stiffen the latter structure i in 1941. Ata a conference on November - 
3, 1941, Shortridge Hardesty, M. Am. Soc. CE, , who was asked to pass on — 

the sufficiency of the design, outlined a dependable deflection ‘method on which | 

he based his recommendations. For two basic structures, Golden Gate and 0 a 

Lions Gate, the two methods gave identical results. | Generel utilisation a 

these parameters will be effected by means of a representation in which the 

and values of are Teferred to curves by 


‘ = Tis = » 

Eq. 122 contains the aerodynamic and ‘damping characteristics 


of the basic structure and indicates, for another structure of similar aero-— 
‘dynamic characteristics, the corresponding amount of structural — 
that will assure stability, in the form of a requirement for the moment. of 
‘inertia of the g girders. 7 For the r reasons 1s stated in connection with Table 13, the _ 


Golden Gate Bridge was selected as the basic structure. Introducing ‘its 

ro into Eq. 122, a criterion of satisfactory apne behavior is obtained 


and the required moment of inertia of the stiffening trusses becomes: 2 
the 
hen By considering m3z and solving Eq. (123 Tz — 


de 
), = 0. 00089 att?) - ~ 0. 0009 .(124a) 3 


(ls 0.00499 (in. *ft?) ~ 0. 005 (12)? Gn? . 


ental able 14 contains the numerical of ™, Ts, and 7 ™ 


The a,-values are plotted in ‘Fig. 22. by A circle is placed around the plots of 
the bridges that moved. Fig. 22 indicates that oscillations occurred in the full 
Tange of cable stiffness variation of bridges It shows that cable 
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4 | Thousand Islands*....... | 28.024 | 140.1 «157.0 
7 | 1,059 | 02703 | 5163 | 191 | 289 
| Mount Hope... 200 | | s392 | 463 | 470 
11 | Triborough...... (1,380 | (0.0576 ad 234 
12 | Manhattan...... ‘147000 | (0.0515 0.917 «180 ba 
16 | Philadelphia............| 1,750 | 0.0615 
17 | 1,850 0.1347 | 12622 | 937 | 707 
18 Present®: 2,300 «| 0.1455 | 181.795 | 1,249.5 | 1,0184_ = 
al 1940. 2,800 0.3262 | 1,044,090 | 3,200.8 | 5,848.6 
George W ashington...... 3'500 — “4 
All bridges are truss type except 2, 21 are plate-girder pes. Item } 
combines both types. Estimated from Data supplied by Glenn B. Woodruf, NM. 


stiffness does not govern dynamic - stability and demonstrates that weight is 
not only the most expensive but also is ineffective i in iii aerodynamic 


umbers from Table 14 


05 1000 200 25 80 


‘Fro, 22.—CaBLE oFr SUSPENSION Barozs, AS BY PARAMETER mn 
_ (Full Circle Denotes Objectionable, Dotted Circle Denotes Not Objectionable Movements) 


123 is ‘Tepresented a curve passing points in F the = 
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fats hi) Center Span Length, 1, in Thousands of Feet Bs 
Fie. 23.—Dynamic OF SUSPENSION BrIpGEs, as BY 
DIMENSIONLESS PARAMETERS 73 AND 74 
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| 
to the first of the parabola (Eq. 24a), representing 
th e ms-values. In In Fig. they” correspond to the second interpretation 
1246), representing the ,-values. The plots of those bridges with s: atis- 
factory resistance to aerodyn 1amic action are near or below, and the plots of 


~ those with objectionable n movement are considerably ¢ above, the curve, the di- 


vergence being nearly proportional to the character of the oscillations. — — When 
of the lay are with the ordinates of the 


re 


TABLE 1: 15. MoMENT OF INERTIA. AND 


(Excess = +; Deficiency = -) 


Eas. s. 123 and 1: Bridge Eqs. 123 1245 


‘Manhattan 4.65 Manhattan 
Philadelphia 


+ 
© 


jos) 


Transbay. . 

Tacoma (proposed). . + 2. Triborough . 

Isle d’Orleans Tacoma Narrows (proposed). . 

Grand Mére. Ambassador 

4 Ambassador i 
-Mount Hope 
Lions Gate 
Bronx-Whitestone 
Lempa Riv 
‘Triborough. . 


DO 


Bronx-Whitestone (proposed) . 
Davenport 

St. Johns........ 

Mount Hope.. 

Lempa River 

Grand 


+ 


Maumee River. 
aldo-Hancock . 


BR 


Bronx-W 
Thousand Islands 


Bronx-Whitestone (present)... . 
Tacoma Narrows (1940).... . 


Tacoma Narrows (1940)...... 


$ 


a 

| oT. Table. 15 ‘lists t the bridges” in the sequence of their dynamic stiffness, as 
_ determined by the foregoing simple relations and the ratios of the required 

~ moment of oningg of the stiffening girders to the moments of inertia in the 

actual structures. Examination of Table 15 indicates some evident incon- 
-sistencies it in both Eqs. 124a and 124b. _ However, the behavior of of the actual 
‘structure is the integrated result of action and resistance since neither the : acting 
sn nor the damping effects are shown quantitatively by the behavior. For 

a this reason a consistent arrangement of the bridges in the order of their dynamic 
stiffness is unknow n. Nevertheless, Fig. 23 disproves visibly “the high 

_ sistency and validity” of the stability criteria presented in the paper. For 
the Bronx-Whitestone Bridge, the author’s minimum requirements (Table 13) 

- ‘result i in trusses of 140% larger moment of inertia than i is sufficient for stability: 
contrariwise for the Deet Isle and Thousand Islands bridges: they result in 
ie trusses, the moments of inertia of which are deficient 60% and 380%, respec: 


tively, as indicated by Fig. 23. (Bridges that 2 are stable without satisfyi ing 
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this criterion have been discussed by the writer elsewhere.’*) attempt, 
how ever, to derive aerodynamic stability from purely empirical considerations — 

should be r regar rded as a temporary expedient for estimating purposes and should 
; @ not take the place | of | specific ex experimental research for each structure to be 


tog 


“design: for a aer maiheeingts: atability i is not app: went from. the results reported i in 7 

the paper, the importance of the variation of the characteristics of the suspen-_ 

bridge system as they influence aerodynamic stability will be 


"This ex examing ation reveals the followi ing: 


1) Cable sag (Vf) determines the fundamental frequency; 


Cable stiffness, deter mines the structural 
“nes Consequently the “requirements of static stiffness can often be met 
ally without t trusses. Fig. 22 shows” that within the limits of its practical 
variation, cable stiffness cannot assure ‘aerody namic ¢ stability. Ro 

(3) Trusses, within economically | practical limits, contribute little to the 
structural elastic stiffness but they can be the principal contributors to struc- 
tural damping. ‘They represent, “however, neither a progressive, nor alwa ays p 
an economical, ‘solution i in the design for aerodynamic stability _e 
es (4) Stay s sy rstems may be. developed into efficient contributors to structural i 
‘damping. Dynamic models overestimate the damping efficiency of stay 
systems grossly. > The Bronx-Whitestone model indicated high damping effect 


for stays” which, installed on the prototype, have found insignificant. 


4 


-Wind- tunnel models equipped with stays should supply basic data for the 
of efficient stay systems . Stay systems will then gain 
in the economical design for aerodynamic stability, 
(5) ‘Suspended structures which are ‘aerodynamically stable, with or with: 
oat efficiently provided additional structural damping, constitute the seconom- 

ical solution of the | design problem of aerodynamic stability. In other w rords, 
the structure whose satisfactory static rigidity has" been achieved the most 


as [economically will also be aerodynamically stable , or can be made so the most 


n- 2Ots«éLhee author’s s equations are useless for the so solution of the practical design. 
An "problem, because neither the variation of aerodynamic damping with arbitrary 
ng "modification of the shape of the suspended str ucture, nor the variation of 


‘or structural: damping resulting from any ‘practical ‘arrangement, can be caleu- 

sis lated by them. In fact, they give i inconsistent Tesults even for shapes similar >. 7 

is to the one from which ‘ they were deriv red, as proved b by the behavior of. existing 


13) ‘The author’s conclusions recommending that secured by in- 
creasing weight and depth must be rejected. The fundamental design n principle | 


of structural efficiency, expressed by lowering the ratio of | the w eight of a 
in the structure to the live load the structure is able to. carry, should nev ever 
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assure dyn amic stability cy sis of the ‘cost items ‘of long- -span wenpension 
bridges clearly indicates the great economic significance of the reduction of 
weight. girders also affect the cost unfav orably. he determination 
of stability, therefore, is is primarily the task of aerodynamic | research which 
should be be > connected directly wit ith the design of suspension bridges until the 
experiments permit a clear theoretical interpretation. 
Corrections for Ti ransactions: In Nov ember, 1943, Proceedings, on page 1384, 


12 and 15 following Eq. 56, change “aerodynamic’’ to and 
minimum” t to “critical,” 4 
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AMBRICAN SOCIETY OF CIVIL | ENGINEERS 


STORAGE AND THE ‘UNIT. 


: 
Don JOHNSTONE, Assoc. Am. Soc. ( C. —The difficulties of win- 
‘nowing out the grains of wheat in technical Neneahaen ei are certainly not decreased 
inventing new Meanings fc for old phrases. Engineering terms should mean 
the same me thing today as yesterday. if new concepts are introduced, new 
terms should be invented. . The author, in applying the term “unit hydro 
graph”: or ‘ ‘unitgraph” toa hydrograph which includes ground-water flow, is y 
misusing a basic definition - Similarly, in applying the term ‘ ‘concentration = 
time” to the time interval ‘ ‘between cessation of ‘runoff-producing rainfall and 
the minimum indicated value of K” he is tying another time-honored phrase 
to a new concept. _ Mr. Bernard’ 8 definition | of “concentration time (“a time ; 
interval at the end of Ww which all parts : of the watershed are are contributing to the 
bry at the point ¢ of observation”) gives that particular combination of words 
a job to. do and should ‘pre-empt it from other “uses.% Other examples of 
titicizable terminology include the use of ‘ ‘runoff” for ‘ rainfall” in the third 
and fifth paragraphs of the “ “Introduction” and i in the third paragraph of the 
section headed “Valley” Storage ‘the ‘Unit ‘Hydrograph,” followed by 
its use in an entirely ¢ different sense in the e second paragraph of the next section; 
and the use of “the channel” in that same section without indicating whether — 
reference i is to the main stream alone or to any ' “channel” in the w ratershed. 
Passing” to other “matters, attention is invited to the s statement in the 
“Introduction” determination of unitgraphs from “multiple ‘periods 
of runoff [sic] * ° els. | indeterminate, as there are always more equations _ 
q than unknowns and a multiple of possible solutions ’ Presumably the author a : 
meant “nore unknowns than equations,’ but, "waiving that point, the state- 


ment is still incorrect. By) unitgraph theory, surface runoff in the zth period —_ 


, Norz.—This paper by C. O. Clark was published in November, 1943, Proceedings. Discussion on a 

this paper has appeared in Proceedings, as follows: February, 1944, by James 8. Sweet, and Otto H. Meyer; 
March, 1944, by L. K. Sherman, Gordon R. Williams, and Ray K. 'Linsley, Jr.; May, 1944, 5 F. Brater, a 
L. C. Crawford, ‘Robert. E. Kennedy, and Victor H. Cochrane; June, 1944, ‘by Franklin Snyder; and | 
% Asst. Prof., Civ. Eng., Ohio State. Univ., Columbus, Ohio (on military leave). 

«Received by the Secretary November 11,1944. 
Vol. 1 “An to Flow,” by Merrill Bernard, Transactions, Am. Soc. C. E. 
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result refin 
+e are rainfalls for the couk 
petiods. and (a), 1), — n + 1) are the ordinates of the 


unitgraph for the ath, (x — 1)th, ---, (x —nt + 1)th periods. flow mea- 


and ‘rainfall measurements were > precise, if the flow (ground- flow, 
Ww ater flow) w ere know n, if the net rainfalls for all periods » were determinable, ‘rang’ 
if the stream behaved exactly accordance with unitgraph is ne 
‘theory, it w ould be possible to derive the unitgraph fom the hydrograph of itself 
multiple- period storm | by successive subtractions. Actually, successive Rath 
subtractions yield nothing of value, as 1 none of the four “‘if’s’” is exactly i in & rainf 
map 
~ accordance with fact and as the p process involves the accumulation of error &B intro 
from ‘step to 10 step. M... he concept is set forth here, how ever, ‘to demonstrate mom 
“that: there are _ actually : as many equations available as use. are unknowns is sat 
(so f far as unitgraph theory is concerned), and that there should therefore be the ‘ 


no possibility. of multiple s solutions from a given set of data if a proper er technique ‘surfa 


of derivation is applied. « One such technique has been outlined by William ee 


. Collins,” Jun Am. Soc. C. E. detai 
It follows (see ‘ ‘Tntroduction”) that. “the practical requirement of finding the r 
an isolated flood hydrograph resulting from a uniform unit rainfall,” which the is ob 
author states “ provides : an almost insurmountable obstacle to the derivation 
- 4 of the unitgraph for very | small [?] drainage areas and very | short unit periods the ¢ 
 Fintervals] of time,’ is no longer a ‘requirement at all. As a matter of fact, ing 0 
"3 a protracted period of heavy rainfal! i is” the ideal period { from which to derive contr 
the: unitgraph, because (a) 1 rainfa'! | is more likely to be nearly uniform over the JB done’ 
watershed, (b) infiltration ‘conditions are more nearly uniform throughout the cross 
period of | rainfs alll, and (c) the fact that rainfall actually extended over the they 
selected unit time interval is s assured. To ¢: carry this argument further would [& Mr. ] 
be » beyond th the e proper scope of discussion; the w vriter acknowledges that ¢ despite highe 
the incorrectness of the quoted thebeats there is still need for a method that # time 
will permit deriving unitgraphs to rainfall records. know 
7 ‘The argument under “Valley Storage,” down to the introduction of Eq. 4, JB techn 
is confusing. Particularly w veleome would be a justification of the caleul 
- that “reduction o of the i inflow rate necessarily would be accompanied by * bes TT 
an increase of outflow rate.” and its explanation are unconvincing, all ar 
and t the writer Ww ould like to see (here as elsew here) proof i in the form of actual i sump 
recorded data rather than i in the form of a rather bir theoretical analysis. 4 
The argument on the extension of storage theory to _ground- water flow 
appears to. “boil down” to this: It would be convenient if unitgraph principle Be mome 
were applicable to ground- -water flow; ; therefore they are /applical able. 9 _W ithout 
_ further proof, the author proceeds te apply them. There is something about [ie upstr 
neatly arranged tables of computations that seems to consecrate the result, a) 


and it. is likely that. computing | a number of ground- water flow hydrographs choice 
this method would confirm a in the Sugge: 


Distribution Graphs from Precipitation Occurring in More One Tie nit,” by 
William: T. Collins, Civil Engineering, 1939, p. 559. 
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JOHNSTONE ON STORAGE 


December, 1 944 
+ 


refinement a’ at all) seems scarcely Ww orth ‘the effort; 30% of curve (c) in ‘Fig. 8 
could be » replaced by the ‘customary sloping straight ie and the change in 


curve (d) Ww vould sc scarcely be noticeable. le, 


. — An even more basic criticism, however, is the inclusion of ground- water 
: flow, by w whatever means ‘computed, in the unitgraph. _ By including it, the — 
, tange of usefulness of the unitgraph is seriously impaired. Its major function 
1 is not— as emphasis in many papers would lead one to believe—to confirm — 
f itself repeatedly by reproducing second-rate floods with more or less accuracy. 
ef Rather, it is a tool to indicate what will probably happen \ when a really: first-rate 7 a 


rainfall comes along g; and it should be left to do that job: unhampered by the 
r “introduction of ground flow factors that will likely be inapplicable at the ne 
e _ moment. _ Suppose the maximum storm hits the v watershed when the ground | 


s is saturated—or frozen. true unitgraph ca can handle that situation 
the “pluviagraph”’ of Mr. Ber nard®); a a composite graph of surface and 


i JB surface runoff cannot t handle it—and it errs on the side of danger, Gc 


“fh Passing to the author’s example, there is much difficulty in following the 
- detail. Table 1 shows K = 15, whereas as the text states that K=9. Also, | 
ig the relationship between the plotted curves of | Fig. 8 and the data of Table 1 

is obscure. ‘Should not curve (d) be the plotting of Col. 

Hower ever, , the major shortcoming of the example i is the failure to explain 


ds the construction of Fig. 6. is is the heart of the entire matter—the draw- 
t, IB ing of lines on a map in such a manner as to define areas that make their 

ve contributions to “the channel’’ in successive intervals of time. How is it 

he done? What_ factors are taken into account? Why do the area boundaries: 


he J cross the ‘streams without bending stream- upward, when logic suggests that 
he [they should take the form of exaggerated contour lines? (See, for example, — 
id Mr. Bernard’s statement ‘the efficiency « of the main channel, through 
ite higher velocities, places headwat ater areas much closer to the outlet in terms of P 


at [time than extensive areas nearer but less « advantageously located. ”) With a é 
_ —— knowledge of how v to chart these | areas, the operator finds the remainder of al | 4 


technique to be a matter of measuring areas by planimeter and of running 


calculating machine; without it, he is at a loss where to start. Tass 
‘The justification for applying the same value of K to the contributions 
call ar areas is not apparent. It would ap appear to be reasonable only on the a as- 
sumption that storage conditions are the same in all areas and that ne no ‘storage 
exists in the ; main channel; but obviously the author does not make s any such ‘ 
“assumption. Rather, he takes for the value of K K the vi value existing at the | 


moment all surface runoff is presumed to have reached “the channel.’ 
It follows: that K for the downstream areas should be smaller, and for the 


upstream areas larger, than the average. 


ult, Me Since ‘delineating the areas on the map is subject to a wide range of personal a) 

pbs 3 choice, ; and : since the use of a single value of Ki is open to question, the w riter ; 

that BH Suggests that a somewhat. different approach | might yield equally satisfactory et 

nt” Tesults w ith less effort. «Using “ ‘concentration time” as defined by } Mr. Bernar 
t, base length of a any hydrograph of surface flow (including ‘the  unitgraph a as 


An dete to Determinate Stream Flow,” 


V 7 Merrill Bernard, Transactions, Am. Soc. 
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‘defined Mr. is to the duration of runoff- producing rainfall 
the concentrs ation time. This” statement includes the instantaneous 
-unitgraph, for which the base length is ‘simply the concentration time itself. 
es ‘oncentration time, by the Bernard definition, can be determined from an 
‘inspection of hy: drograp aphs, w with reasonable accuracy, as readily as concentration 

time by the Clark ‘definition. Having determined it, divide the watershed 
map into the ‘corresponding number - of time-interval areas and determine the 
—. area of each. . The tabulation of _ these percentages, without any 
fate operations, is in itself the instantaneous distribution graph of surface 
It cannot be said that the effect of storage is overlooked in this proposed 

~ method, for it is reflected in the fact that concentration time (the ] Bernard 
definition) takes i in the entire period of time required for a particle of w ater 
to travel from the hydraulically- most-remote. point to the point of observation. 
‘Basically, the analysis seems as sound as that of the author’s s more complex 


p procedure. — It would be > interesting | to determine a unitgraph for the . Appo- 
-ies7. River by this proposed method and test it, _ also, for the flood of f April 


(1937. U nfortunately, the necessary basic data are not available to the writer 
Also, an adequate technique for dividing the watershed map presumably 
remains to be dev eloped is This should be a fruitful field for further i inquiry. | 
“a In ‘making such a test, it would be necessary to assume a “reasonable” 
; rater flow, and the writer would b be content with a , straight lit line j Joining 

_ the point on the » dhgerved hydrograph at which rainfall begins, with the p point 
at which surface runoff ceases. “Fig. 9 indicates the latter 


_ Corrections to Transactions: See the errata published in _ Proceedings for 


February, 251, and Mareh, 1944, Pp. 409. arte 
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‘AMERICAN SOC SOCIETY OF CIVIL. — 


Founded November 5, 185200 us? 


GEOLOGY I HI¢ Il 


“ 


writer that the Society io oathahae another paper in the interest of applied 
geology. The y paper is specific in its application to highway e engineering, but 
the treatment is broad enough to apply to many fields in which civil engineers © 
are active. _ It Iti is of distinet value i in further calling the attention of engineers : 

and educators to the 1e application « of. geology | in engineering development ‘and— 7 
nek The late J. Cc . Branner, head of the Department of Geology at Stanford | 
‘University, Stanford University, Calif., , until his death in 1922, was a scientist | 
of broad vision, and his efforts and accomplishments i in the practical application 
of geology in the fields of petroleum and mining engineering are attested — 
‘through the accomplishments of his students after graduation. — The ¢ curriculum . 
he developed 1 was very broad, in that abridged | courses in civil, mechanical, 
chemical, and electrical engineering were required, with opportunities | to the 


‘individual for metallurgy, petrography, geochemistry, 


tie to Pie in work in a new field in which the professor was endeavoring to ot 
obtain interest, research, and study; namely, Engineering Geology. After 
being graduated | in 1913, the | writer was em employed under Professor Branner — 
and was subsequently recommended by him to other similar positions prior | to 
World War I. During these engagements, primarily directed to the nve 
B tigation of the geology and ground-water conditions of certain areas, it was i i = 

} found that employers had problems involving geology upon which the en- he 


| sineers solicited the writer’s services and advice, = 


Ne Note.—This paper by Marshall T. Huntting was published in December, 1943, Proceedings. Dis 

onion on this Pe has appeared in waa > as ielhowes February, 1944, by W. W. Crosby; March, | 
1944, by Carl B own; April, 1944, by F. H. Kellogg, Berlen C. Moneymaker, E. F. Bean, A. T. Bleck, 
Lyman W. Wood, Philip Keene, and Jacob Feld; May, 1944, by H. E. Marshall, K. B. Woods, and D. J. 
Belcher; June, 1944, by R. Woodward aioore. George E. Ekblaw, and Allen 8. Cary and Ade E. Jaskar 
September, 1944, y Karl Terzaghi, D P. Krynine, Robert F. Legget, and Irving B. 
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FORDES OW HIGHW AY ENGINE ERING 


The of the: Spring alley W ater ( any 


afforded opportunities to do w ork in connection with the 


§ vatertightness and stability of the abutments for the Calave eras Hydraulic: 


ristics of the materials in the fill. Similar were afforded 
i ‘og connection with the foundation investigation for, and the « construction of, 
the base of the Gibr alter Dam by the City of Ss inta Barbara, Calif., in 1915, 
when the writer was inv vestigating ground- water er problems. Professor Branner 
was consulted and, with the advice of consulting civil and hydrs aulic engineers | 
ees contacted during those | employments, the writer attempted the pre prepar: ation of 
a curriculum aiming at collegiate specialization in engineering geology. _ 
After Professor Branner died in M: arch, 1922, the attempt to institute a 
curriculum at Stanford ‘University emphasizing the practical application of 
pology in. the fe field of engineering development tion was 
> 
- applied engineering g geology which he est: ablished ir in associ: ation with Professor 
Branner and in the collection of material as his: experience has broadened. 
His files and library are indexed and cross-indexed under the headings of that 
‘outline, in all of which he has made practical a and technical application. of the 


science of geology and under which he has found much of va value in the engineer-- 


ing literature. The w riter presents that outline here in 1 acknowledgment « of 
: his debt ‘to Professor Branner and in the hope that it will be of | service ‘to 
engineering educators: 


at of Rock Formations. 
_ (a) Character and structure—jointing and fractures. 
(b) Weathering—effective depth and geochemical 
(c) Int Interstitial binder—primary and secondary cementation, vein 
filling, healing of fractures. 


@) Strength, specific gravity, watertightness, stability, solubility. 


Guar of Unconsolidated Formations. 
(a) Source of material and method of deposition. Ay je 
Erosion, transportation, and silting, 
ess (c) Geochemical change and physical change in transportation and 
= deposition, solution, oxidation, carbonization, hydration, 
abrasion, and wear; further changes to be expected. __ 
(d) Structure, porosity, and permeability. 
Water content—drainage and recharge. 
(3) Engineering Seismology. 
ay (a) Recognition and location of faults. — a 
(6) Determination of activity of faults: of 
 (c) Affect of vibrations on engineering structures. __ 
@d) ays, pipe lines, and structures in reference to 
movement, and head of ground water. 
@) Chemical and mechanical properties of movement 


pressure. 
(ce) of water 
as wells, spring 
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a dl. (d) Safe yield and expected life of ground-water titi 


(e) Spreading of water for ground-water recharge. 
ped (f) Drainage of saturated ground. 


(5) Landslides and Earth Movements. 
Occurrence, investigation, and correction. | 
(6) Soil Technology, igri 
(a) Occurrence, and classification of soils. 
(2) Development of soils and geochemic ral changes, from | origin in — 
a rock fragments or products through stages in development. — 
Soil mechanics—the mechanical properties of soils as foundation 
material and for use in engineering structures. 
atural Structural Materials and Methods Adapted in Construction. 7 

Earth-fill- construction, compaction, stability, | 
(b) Hydraulic fill construction, effect of water and transportation. 
(c) Rock-fill construction, durability, insolubility, crushing strength. 


(d) Subgrade for roads, highways, and air ports. 


 @ Geology of dam sites, requirements of structures. 
Foundations for bridges. 
(d) Foundations for buildings : and other structures. 
(a) Investigation, location, design of support requir 
on if the engineering student is given an opportunity to obtain a basic know! l- 
edge i in the for egoing subjects, he w ill enter his professional career well equipped 
to understand natural conditions and the use of natural materials, and, what 
is more important, , with sufficient, know ledge to be aware of his own limitations — 
and to know wherein the services of a specialist a1 are required. _ This lack is re- 


sponsible for a _. amen of failures’ of structures. . The geology student, 


tion that will bring ‘nisl satisfaction and opportunities for service. In ihe 
past, the engineer has made the mistake of assuming that a man who has 


acquired a college degree in geology “knows all the answers. ae | geologist: 

“cannot see underground any better than the next man, but he | can recognize 

certain conditions from surface indications « and, with adequate ¢ experience, will — 

learn how to gain subsurface knowledge and apply basic know 


‘Under the head of ‘ “Problems,” the author lists twelve items in which the 


‘geologist’ services might be pr ofitably engaged. He then discusses these 
ftom the standpoint of experienced engineering geologist. . The writer 
agrees with the first premise, but the testing of material goes beyond the 
‘province of the geologist i in the « opinion n of the highway engineer, and it is true’ 


that the geologist requires additional specialized training and 1 knowledge in ni 


applying his basic training and knowle dge to problems (1), (2), a and (3). ' The | 
geologist must | have more than a basic training and know ledge of ee 


and ground-w water. hydrology, too often lacking in in the collegiate training, in 
order to be of value in connection with problems (4), (6), (8), and (10). ten 7 


- Training and experience can be acquired just as Ww ell, and in as short a . time, 7 


on ay, development work and as in the field o i 
geology, Ww 
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FORBES ON HIGHWAY ENGINEERING Discussions 
-profession—namely, “mining engineering and petroleum engineering. It has 
always appeared to the writer that the failure to apply geological training and. 
“= ledge properly in highway “engineering was the cause of considerable 
a Beginning in 1928 and d continuing well into the depression years, the oil 
a ‘companies s of the W est found that the surface geological structures were prac- 
- tically all determined , overproduction | was causing difficulties in the industry, 
and geological staffs were ‘drastically cut. The older geologists were absorbed 


in administrative and executive positions, but the younger men were dis- 


missed. Many of the latter came to the w vriter to discuss opportunities for 


employment and the broadening» of ‘their activities. - - The advice given was 
that hepsi should take the examination for j aneared r engineer in the St State Highw: ay 

es as an opportunity arose to. 

show that the application. of their special knowledge was of value. Several | 


men took the e examination, , got on the lists, but were not called after long 


It was felt that some preference i in the matter of hiring s should be given ww 


“i having such qualifications, and a letter was written to the Chief Engineer of 
the State Highway Department setting forth this matter and submitting an 
outline covering the application of this special know ledge in highway  engineer- 

_ ing. A conference \ was arranged at the State Capitol at w hich the subject was 
received with enthusiasm by the division engineers in attendance but with a 

lack of enthusiasm m by the executive head The writer was given to understand 


that such a setup vw would require a centralized control, which highw ay engineers 
were trying to avoid, or the retention of a consultant for the direction of the 


oar field men, w which was against their general policy. As far as is now known, 
benefits e derived from this effort. 
“no enefits were derived from this effor 


_ In 1928 the ) Writer ¥ was appointed to a committee of the San n Francisco 
The 


a (Calif. ) Section of the Society to study the foundation and subsidence question 
= 4 over the built-up area of San Francisco. The writer presented the viewpoint 


of the geologist, and these discussions, resulting in a report published sub- 
ie ess by the Section, opened up to the writer a field of application of 


geology and -study—namely, ‘foundations for 


ad The author a the problem of acting as a witness in litigation sil 


problem (12) under the heading ‘ “Problems”). Litigation affords the time and 
“money for extensive and d exhaustive research such as is never available in 


engineering development. and construction. Court presentation of the 


The writer has" sometimes observed d expert w itnesses 
mathematically. the problems of the underground. hen purely mathemat- 
- ical solutions are attempted and broken down under cross-examination, or | 
"matched by one e factually determined, it soon becomes apparent, that nature 


One: is greatly impressed with the degree of confusion the ‘ridiculow 
results achieved when values of reasonable variance from those selected | by the 
a 2 "proponent, as assumption or small-scale test result, are substituted in ina aims 7 
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A mathematical computation, involving factors sidiagiad through assumption, 
or ‘determined by small-size sa sample test, i is no better | than tl the value of that 
assumption or test. Its result might have been beta stated as as the proponents’ 


| mz study of geology reveals that geochemical change breaks down en size 


“ith comparative ‘rapidity. _ Under on one set of ‘conditions geochemical change 
occurs in natural material to. produce a ce certain product and under ' similar, but 
slightly varying x conditions, the change undergone by the same material yields 
y a different product. _ Moisture contents of materials, determined as percentages 
; of their oven-dried weights, do not take into consideration any of the water 
&§ i adsorbed or water that remains in fixed chemical combination. — There is no 
i constant law as to the quantity of water that may be taken up in chemical 
4 combination by the natural materials. For the e foregoing reasons it is difficult 


to accept ‘mathematical calculations based on sample tests. 


There’ is no ‘such: thing as an undisturbed sample. The act of 


‘ f implies. ‘disturbance from ‘its position of confinement a and support, and a new 
relationship to. air—its oxygen, carbon dioxide, ‘and moisture. The degree 

disturbance cai can be ascertained only tl through knowledge of the geologic history, 
a, the physical condition and ground-water condition of the wider area sampled. 7: 
That does not mean that sampling and testing should be abandoned. Boring” 


X and sam ling, and more exploration boring, are needed to give as complete 

a bie ledge of the subsurface conditions as possible; but the limitations of the 
| results ob obtained must be realized and used, within those limitations, to denote 


os action, reaction , trends, and ranges, as tools to work w ith in the exercise of 


a judgment. The writer uses shear tests, compression tests, physical, chemical, | . 
’ 
‘a mechanical, and moisture analyses i in his work for that purpose; but, in bearing- “ | 
capacity determinations in foundation investigations, he depends u upon positive 


cyclic loading of the undisturbed foundation material; and for settle- 


at ® ment he depends upon the experience recorded upon the same or similar 4 
foundation conditions determined through complete boring exploration. - 
d summarize , engineering education has not provided comprehensive * 
ed ‘curriculums i in n applied geology but has gone é “overboard” on ‘ ‘soil 1 mechanics.” 
> 


restriction, of appropriations such as occurred after Ww orld War ar I, _ many Ww 


be seeking n new employment. if they are forced into commercial lines, as 
geologists were during the ear ly 1930's, much training, knowledge, experience, 
and ability will be lost to the profession. 


» i The profession : as a whole does not recognize ‘ends appreciate the importance pore x 
of applied geology i in connection with engineering development a and construction : 


and it has overemphasized soil mechanics far} more than its importance warrants. 
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RE MIN ON HIGHWAY EN NEE 


te 5 Huntting has shown the importance of geological investigations in highw ay 
i. a. construction. In the past few decades considerable interest has been expressed 


9 4 in the value of geology. In various universities the introduction of courses in 


. A. EREMIN,” M. Am. Soc. C. E.4*—By illustrative examples, Mr. 


soil mechanics and geology has increased the curriculum requirements for 
- - Economic considerations fully justify geological study for highway struc- 
tures. , Various devices have been developed to insure safety in structures. ’ 
a Streams are regulated and protected by various s types of embankments. " ‘Bridge 
seats are designed to follow pier ‘displacement resulting fr from the movements of 
4 various geological formations. _ In tunnel construction a temporary lining is 
used to study the characteristic features of the soil penetrated by the tunnel. 
Ww ithout these. engineering devices, geological science is still unable to yield 
reliable results. Further progress in soil mechanics geology may be 
tained through the effort of engineers and geologists s by interpreting 


construction failures. tendency still or to hide the 
real causes of some important structural failures. It is. true that learning 
through failures is an expensive method; but an investigation of the failure of a 
structure designed according to standard assumptions is still the most eee: 
ind effective method of i improving - design and construction. - Many types sof 


. jeshdeete are repeated ¢ constantly b because no one studies the causes sufficiently. 1 


_ The final design of many large tunnels is the result of careful investigation of 
partial failures that occurred during construction. 


Marsna ALL T. Hunrrina,?® Esa. 2sa__Every thinker or writer must inevits ably 

be biased by own particular ped and experience. Hence, the writer 


may have placed relatively too much emphasis on . certain subjects. — The dis- 
cussions charging over emphasis or improperly | placed emphasis are subject to. 
the same analysis, of | course. One man charges that soil mechanics or soil 
science has been neglected ; others object to including the subject of s soils in a 
~ geological discussion. As was pointed out by Professor Terzaghi, soil mechanics 
is not a substitute for engineering geology. The opposite is also true. The 
‘two. should be ‘complementary, w ith soil tests and geological investigations 
combined, to obtain best results, 
_ The experience of at least five discussers shows that clay is not a major 
am constituent of the soils where frost heaving is worst. The formation of ice 
— : layers causes heaving. W herever the water table is higher than the frost line, 
ae - heaving will almost certainly | occur, and the more permeable the-soil the more 
_ water can be furnished to freeze into ice layers. _ Therefore, even gravels may 


on be heaving soils. In fact, the sand and silt soils where the water table is high 


will be the worst offenders for that re ason, and clay i is pre eferable because o of its 


nearly impervious “ insulating” ‘effect. 7 On the other hand, in areas where the 


<_ calif Associate Bridge Engr., Bridge eh. Div. of Highways, State Dept. of Public Works, Sacramento, 
es - % Private, U. S. Army, Student Company G, Student and Service Regiment, Ft. see Va.; —_ 
Geologist, State Div. of Geology, Dept. of Conserv ation and Development, Pullman, Wash. _——- 


2ba Received the Secretary November 9, 1944. 
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HUNTTING ON HIGHWAY ENGINEERING 


ie wi ater table is well below frost level and heaving i is — common the water must 
= rise by capills writy to reach the trouble zone. a In this case, | the Ww riter contends 
1 that clay, mixed with other soil constituents, adds to the capillarity of the soil 
a and makes a given soil a frost-heaving offender, whereas, under identical 
drainage and water-table conditions, the same soil without its clay content 
7 would not cause the water to rise high enough to allow heaving to occur. —To- 
prevent frost heaving, the examiner first must make sure the water table 


below: frost line. Then he must balance two costs against each other: er: 


(1) The | cost. of lowering the water table to below the depth from 


which w ater may rise by capillarity t to the frost level; and 


a i (2) The cost of breaking the capillary rise from the water t table to frost 
lev el by filling with gravel or by some other means. 
i 


le W here sands, silty sands, and silts without | clay ar are soils that heave, ~~ 
3 Ww riter suggests that this fact in itself may be sufficient evidence of a high water . 
table and thus may simplify preventive studies. by indicating directly the a 7 
for vetter drainage—no o other e being applicable. 
reasoning would not necesst sarily apply in the "ase of silts, which have —_ 7 
siderable capillarity even w hen. clay is present in very small percentages. 
_ Admittedly, the writer’s statement (see heading, ‘ ‘Fr rost-Heave Problems” 
‘It has been shown that i in many, if not most, instances of frost- heave rh — 7 


clay i is a major constituent of the ‘subgrade “material” —is in error. 
[It should be amended by the addition of the clause “where the w ater tal table "4 a 


below the frost line.’ Also_ much too vague 


‘tine of research would bea a ‘study of the relation betw een soil particle size » and 
capillarity, using samples taken from heaving soils. AC comparison of the 
effects of the several individual classes of clay minerals would be important in 


— 


“such a astudy. A soil with a certain percentage « of clay of one mineral composi- — 

tion ‘might be expected to. give results different from those e given | by soils with 
the same percentage of clays composed of other minerals. 


ie How the soundness of a bridge superstructure ¢: can be independent of of founda- 


tion ‘materials, as stated by Mr. _ Feld, is difficult to conceive. My A number of 
2 bridges have failed because piers and superstructure were not properly designed 

the existing foundation conditions. In the discussion of bridge foundations, 

the writer did not wish to ‘ “indicate that taking and testing soil samples is 

“The geologist’s sphere of “usefulness in mechanical sampling lies in 

a examining samples; in helping to interpret the samples in terms of dimen- 

'B sions, structures, water- bearing characteristics, and load-bearing capacities 

“jf of subsurface formations: and in combining his interpretation of subsur- . 
he face conditions with features of general and local or form the basis _ = 
“we for recommendations for further sampling and testing.” 


entirely the work of stateme nt (see he: Competen ney 

of Materials for Bridge | Foundations’ ams 
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The actual w and most ; of the 

+ interpretation of physical 1 tests—should be done by soil mechanicians and or 

- it is interesting to paraphrase Mr. Feld’s instance of the | contractor w om 

found himself i in the class of an ex- -contractor because he relied ona geological 

prediction that a certain hill was composed of glacial till, when excavation found 

Se rock instead: A patient died as the result of using drugs prescribed by a 

doctor 1 who made an incorrect diagnosis’ of the patient’s ailment. Perhaps 

the doctor was ‘not very ‘competent to diagnose that type of case, or perhaps 

ii his examination was not complete enough. a In any event, such an example is 

scarcely s a basis for condemning | all doctors or the medical profession as a whole; 

nor ¢ does ‘it indicate that it would have been preferable e for either the patient 


himself « or the druggist ‘to diagnose the case and prescribe the medicine ™ 


= 


| 


but solid rock w w as not in the case Mr. Feld 
himself cited i in the preceding paragraph of his discussion. In most instances, 

| excavation of unconsolidated materials would be preferable; but the q question 
t will be found and in what quantities in 

a a given excavation. The geologic history of the area will contain the answer, 
19 and the better the investigator knows that history the better equipped he is to 
ae In reply to Professor Legget’s question a as to the use of the term ‘ “traprock,” : 
admittedly the term should be avoided, and the w riter’s only excuse is that, 
in an earlier, unpublished, more detailed version of the paper, a table was 
- included from a governmental publication giving the a average results of physical 
tests on thousands of samples of various: rock types. One e of the rock types 
listed wan “traprock. the | published | version of the p ‘paper this table was: 


purpose e of ie paper, “ extrusive ig igneous 
_ The “within limits” which Professor Legget questions should refer to ) the 


Re limits of the examining geologist’s ability and the limits of geology i in general, 
not to the limits of the depth and shape of the deposits. 
oe Moneymaker partly answers his own question: | ys “How should a a state 
highway department provide for, its geologic work?” Mr. Crosby offers 


further: clarification. answer depends upon the amount ‘of geological 


investigation ‘needed, the organization of the highway department, and the 

geological personnel and facilities available and economically practicable. 

a These conditions v vary widely from one state to another; each state has its own 


_ The geologist | whet is to work with engineers when needed should certainly 


be an engineering geologist who devotes his full time and interest to that, 
enka _ He cannot be primarily a petrologist or primarily a physiographer, 

for example, since ‘geology like engineering is. a highly diversified | subject 


geologists have to such an extent that they have little 
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December, 1944 
_ know ledge of, and less interest in, the geological fields n not included in their 
lines of specialization. The engineering geologist must be competent 

most of the special fields and should have a wide range of interests and a el 
“outlook. His willingness own limitations is. 4 


to call for th the assistance | of other specialists. 
oe There are many geologists who think of their science as an end in itself. . 
Obviously, an is to the field of engineering, 
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"SURGES IN PA 
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W. DW waRDs,? } M. Am. Soc. C. E., AND Epwarp Soucek 

Soc. C. E.1°*—It is interesting to note that one discussion expresses doubt 
-Tregar ding the ability of the model to repr oduce s surges that would be caused by i 
operation of the proposed larger lock, whereas another expresses as' astonishment aa 
it was consider ed necessary to conduct verification tests. 
. be Mr. Weggel states, “Te st No. III and Test No. - are of special interest be- a 


cause they the r results of consecutive lock. ations; 
da 


presented i in Fig. 


a In order to explain more fully, reference i is pen Fig. 4(a). The initial 


lock oper ation depressed the water-surface elevation immediately.  Thes second 
_ lock operation was performed 88 min later when the water surface was again 
_ approaching x the low elevation. Withdrawing more water from the canal at this 
y time lowered the water surface and increased the magnitude of the surge. — It 
- should be noted that 88 min after the second | lock operation the water surface 


was low again bt but was not below the elevation caused by the initial w vithdrawal. 


q The surge as augmented by the second operation, therefore, Tepresents the worst 
condition likely to be pr duced by ation of the ex isting locks. 
> i _ Reference i is made next to Fig. 4(b). _ By performing the second lock op opera- 
- tio n 44 min after t the first, when the w ater-surface elevation was approaching a 


maximum, ., the surface was prevented from rising to the elevation attained in 
Test No. III and the surge was neutralized. The choice of intervals between 
lock operations w was not accidental, as implied by Mr. L eypoldt, but was b based 


i i: “one calculation and on n study « of the water-sur face profiles of Tests Nos. I and IL 


 Nors.—T his paper by F. W. Edwards and Edward Soucek was published in . January, 1944, ‘Pro- 
< ceedings. Discussion on this paper has appeared in Proceedings, as follows: April, 1944, by ‘Harold A. 
ie | Weasel; May, 1944, by Harry Leypoldt; and June, 1944, by C. J. Posey. 


1a Associate Prof. of Civ. Eng., Pennsylvania State College, State a 
30 Hydr. Engr., Office, Chf. of Engrs, Washington, D.C. 


ae Received by the Secretary October 11, 1944. wine 
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December, 1944. EDWARDS AND SOUCEK ON CANAL SURGES 


‘Tests } Nos. III and V substantiate Mr. Weggel’s belief that, ‘“Not always is 7 


it the amount of change 6f discharge that determines the height of wave or 
depth of trough.” 7 The data regarding waves on the Ohio River are interesting 


but do not support the opinion that a model could not oa 
identical to those observed in the prototype. 


writers agree fully with Professor Posey’s s viewpoint that comparisons of 

models and prototypes ‘should be undertaken whenever practicable. Although | 
it was expected that the model \ would represent the prototype correctly, it | w as — 

“1 4 considered good « engineering practice to check results by every reasonable _ 


i The entire discussion by Mr. ‘Leypoldt is at vi ariance with the beliefs of the 


w profiles of Fig. 2 do not show the actual | 
ie They are seriously influenced by the findings of Professor King, who states 


: definitely that* the water is moving in an open canal when water is added = 
oe removed. In the prototype, as in all other seiches, the wave form an 
a nearly sinusoidal, instead of the odd assortment shown in 2.7 


7 In. spite of Mr. Leypoldt’s belief that the wave form should be s sinusoidal, weal 7 
i oh, as represented i in Fig. 2, show this to be untrue in the present case. — a 
_ The writers stated that the ‘profiles | should be considered pictorial rath 
than exact, but this qualification \ was made because a profile representing nine 
miles of omal is not precise when based on relatively few observation points. 
_ The st statement was not made in order to age-4 a fictitious | representation of the 
data. ae wave forms represented i in Fig. 2 2 were influenced only by the ob- 
, theory, or preconceived notion on the sub- . 
ject. The agreement f betw een the water-surface profiles and the field observa- 
tions for Test No. II es can be seen in Fig. «6. In this figure the observed 
water -surface elevations from Fig. 3. and the profiles originally sho shown i in Fig. 2 by 


described are undoubtedly to the of the wave. 
The relationships | pertaining to the height of surges formed in an open channel 
= by sudden changes in discharge are based only on the conservation of matter 
4 (implicit i in the law of continuity) and Newton’s Laws of Motion. appli- 
aid cation of the equations i is affected but not invalidated by the initial zero — 
eu pa flow. ti is well know: n that the initial shape of a wave ‘in a channel is 


function of the manner in which it is . generated; whether the wave pence dll 


| in _‘Tegularity in form or te tends to break up as it travels along the channel depends : 
ect eu upon conditions beyond the scope of the paper. Whether or not the wave form — _ 
Wan’ odd assortment’”’ probably depends on the point of view. To the writers, 
I any apparent oddity « can be due | only to an attempt to ) classify the waveforms 


in accordance with a preconceived ‘notion—an error which Mr. Leypoldt 


ranslatory. Waves in Open Channels,” Horace Ww. King, Civil June, "1933, pp. 


21. 
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e writers disagree with Mr. Leypoldt’s contention that the data —— 


» ak seiche of two or more nodes, one at station 1772 and another at station 1456. | ; stu 

Since no significant oscillation occurs beyond station 1478, it is clear that nate ‘ ma 

Ja occurs at station 1456, but here the agreement ends. If a node were located at | «gen 

> a station 1772, the following conditions would exist: (a) T The \ water surface atthe [Ris 

© Elevations Obtained from Figure3,TestI] q Iti 

Lock 

“had 

Bro 


an 70 Minutes add 

| 


"85 Minutes indi 


Fic. PROFILES AT 5-Mix 


node would oscillate less than at stations on either side, (b) 1 the : amplitude » Ww ould of t 
oe increase e progressively from the node to ) adjacent loops, and | (c) the concurrent - core 
s oscillations on the two sides of then node would differ in sign. r _ The water-surfact -Tepl 
ca profiles: of F ig. 2 show that none of these necessary conditions prevail fors Be ally 


‘Mr. Leypoldt believes the to ‘Colonel Brown’s work‘ to be that 
irrelevant. He states that “Gatun Lake proper” does not influence the part< 


4}  4"*Flow - Water in Tidal Canals,” ” by Earl 1. Brown, Transactions, Am. Soc. C. E., Vol. 96 (1982) = 
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‘the canal i dieonnstien’ in the s same manner that lakes do in Colonel Brown’ s 


studies. Itis clear that 1 Mr. Leypoldt means the main bi body of water ¢ approxi- 
mately six miles from station 1478. However, Gatun Lake is considered 
"generally and properly to extend from Gamboa, the point where the Chagres 
_ River flows into the canal at station 1478. The second paragraph of the paper 
and the notation of Fig. 2 testify that this designation was used by the writers. 
iti is Gatun Lake at station 1478 to which the work of Colonel Brown is relevant 
and: which does act in a manner identical to that of the lakes in his studies. . 
is Gatun Lake at station 1478 which was } represented i in the model by a forebay 
in which the water surface was maintained at a constant level. if the wa waves © 
had not been reflected from station 1478 in the manner described by Colonel 
Brown, the » model could have agreed with the ] prototype « only k by chance, the 


ww. © 


chance being the existence of a node de at that si station. 5 
ve Since Mr. Leypoldt’ s seiche theory does include a ‘athe near station 1478, 
iti is necessary to go one step further to show that Colonel Brown’s theory of 
_ negative reflections does apply. _ This step is proof that significant oscillations © 
do not occur beyond station 1478. a Mr. Leypoldt suggests that the existence of 
additional nodes could be verified « easily on the prototype by : additional gages 
‘in the southern | part of Gatun Lake. Preliminary reconnaissance prior to the ‘ 
S _tests showed that significant oscillations did not occur beyond station 1478, but 
to remove all possible doubt, gage a at station 1367 4-50 was established and 
operated during the first two tests. 
The following quotation from the official report on the prototype tests! 
“pli the abandonment of this station: ‘‘As had been expected, the surges — 
beyond Gamboa were small so this station was discontinued after the second © 
test.” This gage located near the point where Mr. analysis 
‘The “narrow, shallow Ww waterway except for the dredged cuts,’ 
Chagres ‘River each have mean width: of more than one-half mi mile. The 
quiescence of these bodies of water. during the surge tests shows conclusively 
that the bodies are large enough to reflect waves negatively in the manner 


N 


‘described by Colonel Brown . The creation of significant ‘oscillations in these 


waters by means of impulses no greater than those caused by filling the locks 


isphysically impossible. 
From Mr. Leypoldt’s comments concerning Gatun Lake proper, from his — 
‘opinion that no conclusions regarding model operations are valid from these “. 
"periments, from his statements that the width does not enter the problem and 
- that roughness i is unimportant, ‘it is not evident that he understands t that the — 
model extended only to station 1478. model did not e1 encompass the length 
of the oscillating system which he believes to exist in the prototype. Ac-— 
"cordingly the remark that it would be astonishing if the. ‘model were not a | good 
--Teplica of the prototype is especially inappropriate. 7 Although it is true that — 
“all parts of the model could have been reduced proportionately i in width on the 


_ basis of grevitetionsl and inertial considerations alone, it is by no means evident 


” “Surge Tests—Gaillard Cut, ‘Report on Procedure and Results,’’ by Edward Soucek, Third Locks — 
Project; The Panama Canal, January 16, 1940, p.3 (unpublished), 
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E roughness vee s the nt number of oscillations, w rhich Mr. Leypoldt concedes, it 
must also have a progressive effect on their r magnitude. . Therefore, neither the 
width nor the roughness may be altered without careful of the 
7 7 results. Since, in a model of the type described, many problems in 1 addition to 
— ' the period of oscillation are investigated—for example, the flow pattern at the 
junction of the ~pass yass for the proposed new lock—Mr. Leypoldt’s sugg sestions 
future experiments have no value. does not appear desirable to 
adopt the restricted meaning of the term ‘ ‘surge,’ as advocated by Mr. Leypoldt, 
a paper written for engineers. 
Except to state that the specific suggestions offered by Mr. Leypoldt for 
— solving the surge problems in the e Panama C Canal are considered impracticable, 
‘no comment on this question can be made a at this time. . The ‘paper wi was 
! specifically confi weed to an explanation of the surges a and to verification of the 
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L.J. M. Soc. C. E The and airfield 
“pavements: will much from this excellent § Symposium, espe- 

cially if he compares the data presented by Brigadier- General Stratton with © 
the specifications for concrete ‘pavements available in the 1930’s. This” 


discussion will be confined mainly to the joint: spacing and design of joints, 


In most cases longitudinal expansion joints may be dispensed with even in 
airfields 300 ft wide. Concrete pavements expand lengthwise from 
eauses—moisture variation and increase in temperature. ‘The. Arlington (Va.) 
tests of the Bureau of Public Roads (5a)*” show that, due t to moisture variation 
alone, concrete ] pavements lengthen 0.0003 in. . per ai in 4.5 years after con- Z 
‘struction. — It is probable that th the expansion may be as | much as 0.0005 in. 


per in. after 10 ary s. The expansion due to an increase in ™ temperature — 


set ‘the a average of the awe een 1 top and 
bottom | surface) when the sun is shining on the pavement in June or July. 
If this interval is 50° F, for exam] Se the expansion from this interval will | 
be | approximately 0.00025 in. per in., or a . total expansion of 0.0005 + 0.0002 5 
= 0.00075 in. per in.; that is, 1 35 in. at each side shoulder in a pavement 
300 ft wide. shoulders will yield, and the only considerable resistar 
comes from the friction of the ground and the w eight of the slab. Assume the 
slab to be 12 in. thick, w weighing 150 Ib per sq ft, and the coefficient of friction 


—This Symposium was published in January, 1944, Proceedings. Discussion 
has appeared in Proceedings, as follows: February, 1944, by James B. Newman, Jr.; March, 1944, by : 4 
Thomas E. Stanton; April, 1944, by W. E. Howland, and David 8. Jenkins; May, 1944, by W illiam E. 
Raymond L. Irwin, G. G. Hibbert Hill, Jacob Feld, and Robert Horonjeff; June, 
1944, by Clarence S. Jarvis, D. D. McGuire, E. D. Parmer, "A. Casagrande, and Harry E. Cotton; Sep- 
tember, 1944, by D. O. Murray, -‘Thomas B. Pringle, John R. Haswell, W. J. Turnbull and William - 
Jervis, Ralph R. Proctor, D. M. Burmister, Ralph Hansen, H. V. Pittman, G. R. ee and H. M. 
Williams; and November, 1944, by Ira B. Mullis, and C. F. Izzard. 


Civ, Engr. and Contr., Chicago; ag tt 
Received by the Secretary October 9, 
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‘to be 2 (5). The total resistance in a slab: 12 in. wide 150 ra and 12 in, 
thick will then be 150 X 2 X 150 = 45 000 ‘Ib This force divided by the 


be -mum stress in the center line of the 300-ft-wide airfield. Such a compression 
a. helps the: concrete as it diminishes a any tensile stresses that may occur. — Because 
of the time _ loading, it would shorten a pavement 150 ft long about 0.2 in. 


Boge: the ‘modulus of elasticity i is s assumed as 1. 33 x 10%, leaving a net increase 


section of the slab, 12 in. X 12 in., gives about 300 lb per sq in. as the maxi-- = | 


in width at each shoulder of about 1.13 in. A drop i in temperature « of 100° F 
will shorten the 150-ft length about 0.0 0004 in. per in., or 0.72 in. _ This is 
expected | tol be distributed among t ‘the ten to fifteen longitudinal joints, So that 
there is a gap ‘at each joint of from 0.048 in. to 0.072 in. Of cours, if the 
“3 dowels are not properly he held in place, absolutely parallel to the direction of the | 
Be esigragerds (for example, if if they should depart from the parallel as much as 
1 in. in 16 in.), the dow els may may freeze this movement in the joints and concen- 
trate ‘it in one joint. On the other hand, they may cause a longitudinal crack 
or at least a high tensile s stress of the or order of 300 lb per sq in. near the center 
of the p pavement. Ww hen the - designer has no confidence in the kind of load- 
transfer device he specifies or in the kind of workmanship he ‘may obtain and 
when he fears that the dowels will not be e parallel, then he is quite justified i in 
providing : a longitudinal expansion joint such as that shown in Fig. 2. 
e size of the pavement sii between t the various longitudinal and 
“transverse joints i is shown in Fig. 2 to vary between 10 ft by. 12.5 ft and 15 ft 
: by r 25 ft. It is generally believed that this size should be limited by the 
“stresses generated by moisture and temperature variation. Another important. 
point has been overlooked by General Stratton, however. individual slabs’ 
are subject to. rocking and vibration due to the heavy concentr ated a | 


_ loads at the edges and ends -_ These W heel loads must be resisted by the w veight 
of the slab acting with a moment arm from its center of gravity to the _i 


~ 


“shock. “Well -dow weled joints s will reduce this 1 greatly, 
dummy joints without dowels have no place in airfields. «dtl is w — known (5) 
"another that joints i in pav vements should have 
load-transfer devices. No reliance can be placed on the interlocking power r of. 
concrete aggregate a few years a after the pavement | is built (5). Hence, the 
greater the wheel load, the greater will be the size of » the individual dab. : 
Possibly a size of 12.5 ft by 18 ftshould bea minimum. 
_ The fear that pavements may be injured by irregular cracking, due to 
warping from moisture and temperature variation, when the individual slabs 
are. 10 ft by 25 ft or larger is based on the. absence of an accepted theory con- 
cerning such a actions and, of course, on ever yday experience. In desperation 
some « engineers ers have advented slabs 10 ft by 10 ft-which are e entirely: unsuitable 
for airfields. Slabs 25 ft long crack because of faulty load- transfer 


devices, unequal fr ost heaving of poor subgrades, or excessive Ww heel loads. 
Approximate of Stresses Produced by —Temperature warp w 
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es _ Roads (56), who found that the highest temperature differential, A, between 


MENSCH ON MILITARY AIRFIELDS 
top and bottom surfaces occurs in the spring and early summer when the sun 


is shining on the tela value, in _—- Fahrenheit, is ‘given’ approxi- 


A=3.75d... a 


as = 3.75 Xd X 0.0000048 = 
For example, a temperature differential of 33° F Ww ion idk for 
A _ A weight slab would deform into the shape of a sphere of radius : 


0 000018 d 4 


he length of these of ‘4 


55,555 in. The ne of Public Roads measured the the warping 


ofa a 9-in. oe 20 ft long, ata temperature interval of 23° F he radius a 
this case would — 55,555 X 33 = 79, 700 in. To simplify the calculation 


and at the same time to arrive at safe and acceptable results, assume ‘that the 
Ww warped | surface | is a cylinder with its axis along the transverse center line of the | 


Subgrade 


> 


slab, ‘thus forming the ‘arch ACB (Fig. 52(a)). The r rise of this arch 


, depending on the soil 
the condition of the near the ond the moisture warping. 
‘The line ADB (Fig. 52(a)) has been drawn equidistant from line EFG to — 
nen the deflection of the slab with the édeal condition shown by line ABC. 
The deflection CD is caused by the loading shown in Fig. 52(b); that is, by a 
“uniform dead 113 Ib per sq. ft ; acting over the entire span of the arch 


| 
.(19), — 
— 
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reactions (shaded areas) — concentr rated at points 


ame deflection caused by the uniform dead load is easily found to be: 


The deflection of the end in relation to to the center p point C caused : 


by the load * = at point H is most easily found. by the moment ar method: 


688 w Lt 


‘he ‘difference— 
8X 16ET 


088 ET the deflection CD. Let b= 240 in; w = 18, 

= 0.75 L; E=4 xX 10°; and I = 729. ‘Then the the deflection CD = 0.0578 j 


The vertical distance betw een points E and. F was measured las 035 in. (5c). 
‘This: is equal to DN; and the total rise of the arch, CN = = 0.0578 + 0.035. 
0928 ‘in., as “compared with the 0.09 In. found from 


temperature interval. mB other w ords, agreement between theory : and pr: vetice 


is very good. - The bending moment at the center of this arch, 1 ft. wide, e, is 


found to be: X = 33,80 800 0 in-lb, leading, to. flexural 

7: 

_ of 209° Ib per ‘sq in. with tension at the bottom, as compared with a 
7 flexural ultimate strength i in paving concrete of rar ‘ely less than 600 Ib per sq ir in. 
_ Furthermore, this stress should be relieved considerably by the action in a_ 

-_ transverse direction, which the w riter has neglected, and by the continuity of 
good jo joints equipped with p proper load-transfer devices. ‘The Bureau of Public 
Roads_ estimated the largest temperature stress in this slab based on 
someter readings at 350 Ib | per sq in., but admitted that there might be some 

_ reservations concerning the assumed modulus of elasticity. ; For a differential 

' of 33° F, a similar calculation ion and an assumption of L' = 0.8 L would result 

in a flexural stress of 250 per sq in. 

due to “moisture variation was investigated for the ‘Bureau 

_ Public Roads by ‘W.K. Hatt, M. Am Soe. C.E. (61). He found the maximum — 
differential change of length betwe een the top and bottom surfaces of a wet 
 g subgrade on on a dry summer day ' to be 0.000176 in. per in. in a slab 7 in. . thick 7 


7 ina in a beam 8 in. thick, and (it may be assumed) also for a 9- in. slab. A weight- | 


less slab w saab tellewee into a sphere, », convex downw vard (see line ACB, Fig. 
53(a)), with a radius r= _— 51,200 in. and a ‘rise (line NC) of 


0.14 in. for L = 240 in. Assuming again a » eylindrie: al warping except. 


that the e dead Sinks of the cylinder is sania ‘ted by the subgrade it in n the middle 
w 


of the arch, these upward reactions may be considered t to be nearly 
4 uniformly distributed | on a length L' = 0.6 L without significant error. The 


— deflection CD of the slab due to its uniform ¥ weight of w per square foot and 


ti the reactions = per linear foot onthe length L’ may be found in a manner 


=) similar to that described for temperature warping: The deflection due to dead 
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load = , when w = and the deflection | due to the 

flection CD resulting fr from time loading = (1 — 0. 408) - ssa py = 0.0082 

in, in w hich E= x x 10, The measured vertical rise from the transver 


joint edge to the center of he deh one: be estimated (5d) at less than 0.04 in 


giving the rise NC : = 0.0982. Ea 0.04 = 0.1382 in. Ww ‘hich is practically the was 


as 0. 14 in. _ computed theoretically - from length changes. The maximum 


bending moment per linear of width, M L = = 26,000 in- lb, 
leading to ‘a flexural stress of 168 Ib per sq in. with ‘uaadiap on top, whereas 
the ma aximum stress from temperature | Variation was found to be 209 lb per 


sq in, at the bottom. * Fortunately. these stresses are not a additive, except for 


the small effect of temperature warping at night. aaa nei 
‘The: foregoing speculations seem to indicate that w arping stresses increase — 


“somew hat less than the square of the length of the individual slabs, and, as the 


deflections i increase V with the fourth power of the slab length, ‘soon a stage w ill _ 


‘reached where the points D in Figs. 52(a) and 53(a) will rendh the subgrade, 
Rano will probably occur in 10-ft-wide slabs when the length is about 25 ft: 

any increase of length will probably n not increase the warping stresses to any 
extent. Much lower stresses occur in 1 transverse directions at the center, and 
| this fact should be taken into account when it is decided to use reinforcement 


steel to to strengthen the slabs at the bottom. Reinfor orcing bars are advantageous 


if placed ‘at the bottom in a Jongitudina al direction n along the edges of the slab 
: and in a transverse direction at the ends. — Where top reinforcing i is selected, 


_ the bars should be of high arbon steel and preferably %) in. in diameter to 


avoid ‘surface ¢ cracking. 7 Such bars have been used by the writer with great 


pepe for Dé wvements designed Stratton 


“choice of of positive holding soy resilient fillers for expansion joints. “7 


Bureau ¢ of Public Roads has shown (5) that q-in. dowel bars are not of much 
value heavy loads, even 12 in. on This has" also been 


= 7 become — by the w riter) while the expansion joints which have been — 


e become close. a The writer has shown that dowels heavier than 34 in. are needed 
y for heavy w mn loads, especially such as occur on airfield pavements (62) and 
e that it i is of the utmost importance that the dowels be held positively, by other — 
‘means than w wire and stake am. in the correct p eon po to the 


s — 
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& 
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loads, dowels of 1s i aimi 
have to be specially reinforced to prevent the | 


edges: of the slabs from cracking due to the e ntrated bearing pressures _ 


(8) “The Structural Design of Concrete neni by L. W. Teller and | 
. Earl C. Sutherland , Public Roads, Vol. 16, No. 8, ‘October, 1935, pp. | 
145-158; Vol. No. 9, November, 1935, pp. 197; Vol. 16, No. ‘10, 


ecember, 1935, pp. 201- -221; Vol. 17, No. ‘September, 1936, pp. 

148-171; Vol. 17, No. 8, October, 1936, pp. 175-192. (a) November, 

1935, Fig. November, 1935, pp. 169- 197 November 

Fig. 26. (d) November, 1935, Fig. 37 

(61) “The Effect of “Moisture on Concrete,’ by W. K. Hatt, Transactions 

Soe. C. E., Vol. 89 (1926), p. 
(62) “Joints Concrete L. J. ‘Mensch, Chica 
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CHARACTERISTIC REDUNDANTS USED FOR | 
_ ANALYZING STATICALLY INDETERMINATE 


RALPH Ww. ‘STEWART, ‘MYLE J. 


W. Srewarr,” M. Am. Soc. C. E. moments in the 


a full page or more of and form-— 
containing three X-values, two 6-values, one y 


_d-value, one k- -value, and nine Q-values. These 
moments can be obtained from two ‘o readings | at 
one setting of a . slide rule and a little e easy men- 
tal arithmetic by sketching a traverse of the elas- o 


tie curves es of this frame (Fig. 8). This" traverse 
diagram shows at a glance that three sevenths of | | 


the e sum | of the moments in either | column goes to 
the top « of the column and four sevenths to the 
bottom . From symmetry and the fact that | 
sum of all column moments must ‘equal ‘the 
eral force multiplied by the frame height: oo 


Moment a at A= 3/7 x 6 P = = = 2.5 714 


Moment at C = 4/7 X6P=2.4286P 


(check) = 6.0000P 


This « one-minute solution uses two constants of flexure, one moment balance, a 


gt 

and one traverse closure. The > author’s solution re requires three simultaneous 


_ Notrse.—This paper by John B. Wilbur, Assoc. M. Am. Soc. C. E., was published in June, 1944, Pro- 


- gedings Discussion on this paper has appeared in in Proceedings, as follows: September, 1944, by Charles A. 


10 Bridge Engr., City « of Los Angeles, Los Calif. 
10a Received by the Secretary | October 2, 1944. = “ty 
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‘not of value as a practical method of computing stresses. _ 
J. Houiey, Jr.,"" Esa. A new method the solution of 

indeterminate str ‘uctures is suggested in this paper. In addition, an entirely 

ug new line of ré reasoning is ; introduced vy which may yield valuable concepts and 
methods beyond those illustrated by the author. _ Practical application of the 
method, along the lines suggested, should prove it more advantageous than 
= conventional methods for some special pr ‘oblems, and, ‘in n addition, from such 
a application may develop even more potent analy tical methods, For this 


_ reason, it is to | be hoped that some members of the profession \ will be encour aged 


to apply the method to practical problems. 


The method of characteristic redundants for the solution of statically. 


indeterminate trusses offers « a reduction in the effort of solving simultaneous 
equations. This gain results. from. (a) the fact that only one set of simultaneous 
Vas equations need be solved, regardless of the number of different loadings for 
_ which the ‘structure is to be analy zed, and (b) the fact that there is one less 


pes ‘equation in the simultaneous set that must be on than. there would be if a a 


antages of the method are ‘that it requires the solution of a 
for the determination of characteristic forces, and that 
it also requires a substantially larger amount of te tabular work ‘for the | determi- : 
nation of the participation factors. — With these relative adv: antages and 
disadvantages in mind one might use ~~ following criteria in deciding w nether 


the author’s method would be worth while in a particular case: 


if the structure is indeterminate to a small degree, the will 


not prove advantageous, even if there are many loadings to be investigated. — 


2. If the structure is indeter minate to a large degree but only “one, or at 

‘most a very few, separate loadings must be considered, the method will not 
proves advantageous. 


loadings 3 to be the method may 


tT The writer realizes that the foregoing criteria are very -erude, but if if the | 


analyst follows some such reasoning before star ting he may avoid an unfair test 
of the author’s method and may also save himself considerable effort. Of — 


— the experience resulting from actual analyses ¢ can refine these, 


riteria to a point w here they will g give ¢ critical combinations of the. degree- ol 
cy and number-of- loadings. 


Although the w iter can visualize no ‘0 case in w hich i it would be practicable 


to use characteristic redundants to determine deflections, it may be. of academic 


interest to point. out that this could a be done. If If the deflection of a point : 


7 
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at which there i is! no load i is desired, or if the deflection at a loading point ina a 
- direction other than that of the actual load is desired, the characteristic set _ 


_ ould include a characteristic load at the point in question, in the direction of : 
the desired deflection, as well as a characteristic load at each loading point | in _ 
the direction of the actual load bit The characteristic load at the deflection point 


will be removed by the participation factors in the superposition process which ae 
_ results i in the actual loading - The he w riter believes, how ever, that such a pro- — 


“cedure Ww vould be far more laborious than any of the conventional methods of 

determining deflections. This i is in no w ay a criticism of the p paper, which (as 


its title indicates) a method for analyzing indeterminate structures 
(not a method for finding deflections) and, for this purpose, the writer — s 


M. 12 Esq.!—The skill displayed i in producing this intricate 


and interesting solution to a complex problem i is deserving of high | praise. The 


aim of the author, which 1 was to find what he calls “the heart of the — : 
indeterminate solution” (see “Synopsis’ ”) and by means of which the analysis 


of any indeter minate f frame i is reduced to the solution of one system | of simul-_ 

taneous equ: ations for all loadings, is a laudable one. © Unfortunately, the in 


tion obtained has not simplified the analysis of such frames. _ On the contrary, 

for most cases , it has added to the complexity of the analysis. At seems quite | = 

within the realm of possibility, how ever, , that some future coven of - 


method may serve to simplify the steps. 
it would appear that the W ilbur method i is in much thes same ‘category as the 


“relaxation, method,’ proposed by Prof. R. Southw As an extension 


“of the principles of moment distribution to the solution of multiply redundant 
“trussed frames, Professor Southw ell’ ’s method was 3 first developed to ‘simplify 


the analysis | of. dirigible frames. In it, , stresses are determined by a series of 

approximations, but, unlike the method of moment distribution, , thea approxima- 

t | tions do not converge readily. In fact, considerable skill and patience are re- 

st quired i in solving even some of the simpler frames. In spite of this, the method | 

a | seems to hold promise and has the important adv antage that it eliminates, 

il completely, the : solution of simultaneous equations. 7 It may be used to analyze 

= problem w hich normally involves the solution of simultaneous equations. — 
As w “ith the method proposed by Professor W ilbur, it has its g greatest : advantage © 


” hen used to analyze one frame for several different, separate ty] pes of loading. © 
Brskty," 4 Assoc. } C. An attempt is made in this 
paper to apply to statically Sidtedintetiinale structures a method of calculation — 
_ which is highly successful in various fields of applied mathematics. | The 


fundamentals of this method, as presented t by the author i in Sections | 1 to i 


12 Asst. Prof., , Civ. Eng., Univ. | of Toronto, Toronto, ‘Canada. 
12a Receiv. ved the Secretary Nover ember 2,1944. 
18 ‘*Relaxation Methods in Engineering Science,”’ by | Southwell, Press, Oxford 1 1940. 


4 Engr., Consoladated Vultee Aircraft Corp., San Diego, Calif. 
‘a Received by the Secretary November 8, 1944. 
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and 16 and 17, can be found i in various good mathematical treatises -15,16,17,18 the op 


material in Section 16, however, lacks precision inasmuch as Professor Wilbur les 
does not apply to his ; problem the classical properties of definite quadratic th 


forms which constitute the basis of the method. This attempt to transpose ‘sir 
- into the field of statically indeterminate structures a method which is essential in sti 
other fields i is not s successful, due to the fact that the problems to be solved in the els 
ease of statically indeterminate structures are considerably simpler than those a 


ot which the method has been developed. . Only in some special cases, most. th 
_ of which are already well ll known, does t the method present definite advantages. 
- Similar developments. were presented in 1907 by S. Miiller!® and i in 1944 W. M. 
Kincaid and V. Morkovin” presented a straightforward method of ‘orthogonal-_ 
ination” limited to the case of rigid fr frames. 
* P To define the type of problem examined by the author, it must be remem- | 
- bered th that. the Strain energy of of a system is a quadratic function of redundants, 
and can be 1 represented by ellipse (two redundants), a a -quadric (three re- 
-dundants), o or a a generalized quadric i in a _ hyperspace (more than three redun- 
dants), these curves or surfaces corresponding to a determinate value of the 7 
strainenergy, 
quadratic fenction | is minimum at the center of the quadrie, the o co- 
“ordinates which represent the actual values of the redundants. solution 


of a system of linear equations deduced from the condition « of minimum str “ 


sa energy is tantamount to the determination of that center. re 


2s 


_ Actually, the author examines two methods of determining that center. a 
: ‘The first is the study of the principal » axes of the quadric, which i is of interest tl 
en these principal axes are needed as such— —for instance, in problems of —( 
vibration, in which these axes represent the periods of vibration. When these 0 
- quantities have no special interest, this method is much too complicated, since : nd 
it requires the solution of an algebraic e equation, _ the degree of which is equal a 
to the number of redundants. The author has devoted more space to this ba 
method than it deserves. 
The second method, Ww vhich is called orthogonalization of 
linear equations, is the study of an arbitrary set of “conjugate axes of the same ¢ 
uadric—a much simpler problem than the ‘first. “ Once these axes are de- n 


unknown 

method, to be an improvement, must be shorter than the conventional = Be 

- in nvolving tl the ¢ direct solution of linear equations. © It substitutes for the solution a 

 ofn linear equations with n unknowns the solution of the sequence, n - —1 equa- i 
| = 


4 

fined, it is extremely simple to find the ¢ coordinates of the center, which are 
"given by independent linear equations, containing one unkno 


-15‘*A Course of Modern Analysis,”’ by E. T. Whittaker and G. N. Watson, Cambridge, England, 1 
--—-s-:16 "Die Differential- und Integral gleichungen der Mechanik und Physic,” by P. Frank and meee von 
Mises, F. Vieweg & Sohn, Braunschweig, 1930; American lithographed edition, Mary 8. Rosenberg, New 
_-___47*Methoden der Mathematischen Physik,” by R. Courant and D. Hilbert, Julius ‘Springer, 1 Berlin, 
2d Ed., Vol. I, 1931, ChaptersI=ndII, pp.1-50. 
3a ‘*Mathematical Methods in Engineering,” by Theodor von Kérmén and Maurice. A. Biot, McGraw- 
—- Book Co., Inc., New York, N.Y., 1940, Chapter V, pp. 162-216; Chapter VIII, pp. 323-364. 
ees 19 ‘Zur Berechnung mehrfach statisch unbestimmter Tragwerke,” by S. Miller, Zentralblatt der Bauver- 
20**An Application of Orthogonal Moments to Problems in Statically Indeterminate Structures, by 
M. and Vv. Quarterly of Applied Mathematics, Vol. I, No. 4, , January, 1944, PP. 
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BESKIN ON CHARACTERISTIC REDUNDANTS 
_tions, n—-2 equations, ... 1 equation. A comparison n of the number of 
. operations shows that the method is advantageous only for four unknowns or 
less. | For these small numbers of unknowns, the additional work required by 

the preparation of the ‘method offsets the theoretical advantage, except when > 
simple formulas can be established for all problems related to a given type of 
structure (for example, the centroid and principal axes of fa beam s¢ or the the 
elastic centerinarigidframe)) = 
__ Moreover, the most efficient provision fo for solving simultaneous equations is i _ 
the Gauss method, a remarkably simple application of which y was presented by — 

P. D. Crout®in 1941. The Gauss method is a systematic method of eliminating © 
unknowns,” which is at the same time a method of orthogonalization. This 
method i is so expeditious that other methods are never justified except, when > 
special physical interpretations make them easier to understand, when 
further her computations require so some intermediate steps | which the Gauss n method 
by-passes. few examples are offered to demonstrate the special conditions 
under w hich | orthogonalization along a definite physical pattern is useful and — 
also to indicate when research involving the aXes is justifiable 

The problem involving combined | bending and direct load is a case _ 
point. When moments are taken with respect to the centroid, the property of 
-orthogonalization will appear ; that is, moments th then en produce only rotation and - 

“direct loads produce only translation. 
“ In the general case of bending, the principal axes of inertia have the property a 
orthogonality—namely, that bending around an axis produces 


around that axis. In this case, pr principal directions are , defined, thus r requiring 
the solution of an equation 1 of second degree, although, in the previous condition — 
(combined bending and direct load), principal directions are not defined because 
of their lack of physical significance. In combined bending and direct load, 
an orthogonal system of the thi 
A more complex case of orthogonalization, introduced in the analysis | of 
: open cylindrical shells (aircraft structures and cylindrical roofs of the Zeiss- 
-Dywidag type), is that of torsion bending. — When these structures are sub- 
- jected to twisting forces, transmission of torque produces axial stresses, which 
_ can also be produced by externally applied direct forces with no resultant and 
“no resulting couple. — These direct forces constitute a mechanical entity 


which the writer has called a ‘bicouple.”23 The distribution of stresses under 
arbitrary forces is determined by a method of orthogonalization, which defines, 
‘ successively, the direct resultant, two couples at the oe centroid (along principal — 
axes), and a bicouple i in reference to the shear center. The 

is of the fourth order in this case. 

oA well- known case of the third order i is that of hina: This. problem | is 
_ generally treated by the use of an elastic center, but : sometimes it is presented 
j as a case of orthogonalization (complex bending and direct load), and is eal 


at “‘A Short Method for Evaluating Determinants and Solving Systems of Linear Equations with Real 
: aa Complex Coefficients,” by Prescott D. Crout, Transactions, A.I.E.E., Vol. 60, 1941, pp. 1235-1240. 


eo _ *‘Methoden der Mathematischen ‘Physik, ” by R. Courant and D. Hilbert, Julius Springer, Sette, 


Stress of Open Cylindrical Membéanee,’ ” by Leon Beskin, Journal of the Aeronautical 


Sciences, ‘Vol. No. 4 4, October, 1944, PP. 343-355. 
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™ _ The writer has also examined the case of the two-story f frame, leading to an 
_ orthogonalization of the sixth order, * which reduces to a fourth- order and a 


second-order problem, when the frame i is symmetrical, 


__ Incomplete orthogonalization is also used. The three- -moment equation | 
for continuous beams is an n example: All moments at supports are e orthogonal, 
_ except moments at three consecutive supports; the two-moment equation (con- 
point. method”) is an | improvement on the three-moment theorem, since” 
only two consecutive moments are not orthogonal to each other. | The method 
can be applied a _ second time, and new moments can be perry which are 


entirely orthogonal to one another. 


ih all these problems, o  orthogonalization is is  justitied because it has a physical 


“reduces: the work to solve without 

_ introducing additional complications in the preparation of these equations. ns. 

: — It is obtained by a a rational choice of unknowns w hich, “once established, does 
>. not require either any additional computation (eentineanen beams), or very little 
additional computation mee, elastic center), every time the method is 

Sometimes, also, the of unknowns facilitates handling equations 
os by reducing the values of the rectangle terms (such as 54. in Eqs. 2) as compared 

with those on the first diagonal dss, Then, the solution of the 

- simultaneous ‘equations i is less sensitive to small errors in the determination of 


6 and in the ‘subsequent numerical This makes it 
The redundants 


are chosen i in n such 2 a ‘manner r that the regions of the structure on which | they act 
ere are reduced to a minimum. 7 For instance, considering a struc- - 
Bh aseah on elastic supports of various rigidities, it is inadvisable to use as an un- d 


known: the reaction of a very flexible support, which greatly increases : all the 


‘cae a of the simultaneous: equations SO that the coupling, represented by 
is nearly equal to 1. 


The foregoing comments show the broad a advantages and limitations | 
_ orthogonalization. ‘ Details of the method presented by the a author will next be 


d. 


has Sections 1 to 13 are devoted to a method of 


correspond to zero 4 The author applies this method to condition 


— which all masses and mass moments of inertia are equal to . _ Then the 
- coefficients « are the reciprocals of the. squares of the frequencies w a 


i= 


“4 “Rings Decked Airplanes,” | by L. Beskin, London, (publica 

2 % **Moments in Restrained. and Continuous Beams by ‘the Method of Conjugate Points,” by L. H. . 
Nishkian and B. Steinman, an, Transactions, Am. Soc. C. E., Vol. 90 — “206. 
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an | ‘Eq. ‘76a shows that the set of conditions , expressed by Eqs. 1, is unduly 


a8 restrictive, and that the actual set should be : 
& 


in which m; is a generalized mass at point j (mass if Aj; is a translation and 7 
is a force; ar sail mass moment of inertia, if Aj is a rotation and P51 is a couple). 


4 Then, n solutions are found in exactly the same manner as the author’ Ss 


solutions, the latter corresponding tom;=1. Itis obvious that, when couples — 
and forces are used simultaneously (Section 13), the m; = leads to 
coeflicients: « which have indefinite units, since they are for. translations 
and forces and for rotations and couples. Thus, the author’s solution 
3 is ‘not homogeneous and changes with the choice of the unit length. _ For’ 


inst: ance, a cantilever loaded at the end by : a load P and ac couple C has a deflec- c- 


is 

ct Tf K with zero dimension are chosen n—for instance, 
o's - 
n- 
he 
by | 
be +8 


Remarks Concerning Solution of Eqs. of the statements in § 
16 must be examined. For ‘example, concerning ‘the first paragraph of that. 
section, it is well known that the strain energy of an elastic solid is always a - 
definite Positive quadratic function, either of the forces or of the displacements. 5 
‘This means that all coefficients x ave veal and, moreover, positive. | In other 

words, all vibration modes are real. Thus, an elastic system of the type a 7 
amined is always stable. This statement is true when the strain energy is of 
form 3 A). Then, the effect of the change ¢ of configuration under load 
is. neglected, which is not the case when problems of elastic stability are ex- 


amined. — It is understood that such cases are not considered by the author, 


since he uses Eqs. 2 to define displacements. 
For the case of roots Professor Wilbur’ examination is incomple 
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to the ilies by. eliminating other thus: 


Eq. 81 has — roots only 


which requires that that ai. 2 = 0. | 
: = The physical significance of Eq. 82 is th that the tw 0 ) corresponding frequencies. 


be equal. This will happen in ‘any system, if masses are chosen in 


a = manner, an will occur for a definite | configuration when the dimensions 


For instance, in Fig. 9, let a n a mass be placed at the elastic center, which : 


Or 


Se 


with the extreme end, E, of the beam. The deflections are: 


V2 
“CRT 3y2 


If the radius of gyration of the mass is (* = 0.486 the 
two frequencies are equal. mes the author’ s method, the two values of x coin- 
when # = A, ,making 34 [2 = 2.06. This is independent of the 
unit ‘Of length. Thus, if the | computation is in 
“<4 inches, the cantilever equal to 2.06 in. has tl the 
exceptional property that its coefficients kK are 
equal. if the computation is in yards, : a cantil- 
ever 74 in. Jong] has the same exceptional property. 
‘It is obvious that this is a characteristic that has 
9 nothing to do with the structure, but merely with 
the choice of units. This effect can be eliminated 
by) choosing other units. In other words, the ‘condition imposed by ‘multiple 
— roots is generally of no significance for the problem e examined by th the author, 


despite the fact that it has a physical significance in | the problem o of natural : 


a ™ When consistent units are used, multiple r roots may have a physical sig- 


nificance. For instance, , if the ‘problem of general bending and axial force is a 


— 


W 
e 


examined in connection with a beam having | equal moments of inertia in two 

4 directions, there : are no principal directions ; or, more exactly, all sets of mu- 

tually perpendicular ‘axes can be assumed to be e principal directions. It ; is 

simple and logical, in that case, to choose the principal axes. of references as’ 


— axes. examines two so solutions (see Eqs. 59 and 60)—P Pas 


al axes, The 1e logical choice in that case is ‘Pay 0; and 
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0, Pee. Then, instead of defining orthogonal systems by such complex formulas 
as orthogonal system can be defined 


Pr. = Pro: : 


sacs has i in mind a special definition of independent components of displace- 


ments. ‘The latter are defined conventionally a as the differences between the 
total number of displacements and the number of equations of statics. More- 
i : over, for the example show! n in Fig. 4, the 1 author indicates | (see third paragraph, 
f Section 13) that there are three iabegondent components of distortion, and it as 
| is well known that a built-in frame is redundant to the third order. — Thus, the 
| example is not consistent | with the first sentence of Section 13, the Meaning of _ 
| which is obscure. | 
; Alternate Basis.—The approach used in Section 17 is much simpler than 
4 that used i in the remainder of 1 the paper, and (as has been | indicated) i is the cor-— 
| rect method for nearly all the problems considered here. ‘ It is correct to state 
¥ that a judicious choice of arbitrary values simplifies ‘the determination of 
dependent values ; but ‘simplification i is generally obtained not not by using values 
—Oor +1, as stated by the author, but by finding orthogonal group: groups from physical 
considerations, such as symmetry, antisymmetry, and polar multiaxial sym- 
metry. — In Section 14, the author examines such cases of simplification, ‘but 


does not connect them with Section 17. 


It is useful to clarify the connection between the two methods, seemingly 

independent in in the e author’ S presentation. This connection is very simple 1 when 

| _ the author’s nonhomogeneous special relations (Eqs. 1) are replaced by Eq. 76b. 
_ Then, if a system > is orthogonalized in an arbitrary manner by the second 
_ method, there will be a system of values m; which will satisfy n systems of 
Eq. 760, written successively for the n independent sets of values obtained | by 

the method of orthogonalization. — in] Figs. 3 and 7, the values 1/m; can be 
chosen equal to 1, 3 0, and the values w? can be taken as 1/100, 0, 1/100. In = 
this case, there are indeterminate quantities | 0/0, which can be defined by con- 


-tinuity from orthogonal sets differing from the solution selected, by infinitely Z 


_ _ These considerations show why the second method i is simpler than the first: 


‘The system of masses ™; can have ¢ arbitrary y values in the second method, 


whereas these masses are e equal to 1 in the first. _ This determination of masses’ 


the solution of an equation of n, otherwise purely 


jae have been eliminated for a given set of loads (second rhs of the 


equations), the supplementary work necessary to find the redundancies in other = 
loading conditions is extremely slight. _ When many loading conditions are con- 

templated (more than) the order of redundancy), the use of unit conditions | for 
= second members reduces the work in great proportions, and, if Mr. Crout’s” 


method i is applied, this supplementary w 
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paper have recently y attracted much attention, and the quite complete treatment 
: te author has given is most welcome. i It is hoped that in his closure he will 
explain the relation of his results to those of A. A. A. Kalinske and C. L. Pien,® 
Assoc. Members, Am. Soe. C. E., which seem to be founded on an » otedy ad 
different approach, but which also lead to laws of logarithmic vertical distri- 
bution of bot! both longitudinal velocity and sediment concentration. It may be 
_ stated, however, that G. I. Taylor,**i in his‘ ‘modified vorticity transport theory,” . 
deduces that velocity distribution for two- dimensional flow i is eycloidal (except 
for the laminar boundary laye er at the bottom). . Although the cycloidal cl 
_ tion is entirely different from the author’s logarithmic one, the numerical values _ 
differ: hardly enough to be “distinguished experimentally, The writer 
smaused himself by deriving a simpler distribution by thinking of the laminar 
two-dimensional flow of oil of high viscosity | with a thin layer of water between 
it and thet bottom. . The velocity profile, of course, proves to be the parts 0 of tw o | 1 


parabolas, the lower one in the water straight line, and the 
. one with its vertex at the upper surface of the oil. This gives an ap- | 


. proximate picture of two -dimensional turbulent flow, the water representing the 


laminar boundary layer, and the oil the turbulent part with its high (mechanical) 
viscosity. It is only : approximate because the mechanical viscosity is is not uni- 


form, but is a function of the depth. 


Neither this rough approximation nor the more accurate cycloidal 


logarithmic equations f fit the case where the flow i is not ‘two- dimensional, how- 
ever—as, for example, in narrow rectangular channels or rivers of 


Nore. —This paper by Vito A. Vanoni, Assoc. M. Am. Soc. C. E., was published in June, 1944, . Pro- 
Associate Prof. of Mechanics, Ohio State Univ., , Columbus, Ohio. 

%5**Experiments on Eddy-Diffusion and Suspended- Material in Open Channels,” by 
A. _A. Kalinske and C. L. Pien, Transarwons, Am. Geophysical Union, Pt. II, 1943, pp. 530-536. 
- 36**Flow in Pipes and Between Parallel Planes,” by G. I. Taylor, Proceedings, Royal Soc. of anaitinds, 
"Series A, Vol. 159, 1937, pp. 496- 506. 
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depth. — As i is well known, in these cases the maximum velocity « occurs some 
distance below the surface, , and Eq. 15a could not possibly apply. The only 

satisfactory explanation of this depression, which the writer has seen, Was 
developed independently in 1882 by the late P. Stearns,?” M. Am. Soe. C. E., 

America and Max in Germany, although Mr. Stearns 

“that the explanation had been suggested to him some six years earlier by the 7 
late Hiram F. Mills, M. Am. Soc. C. EK. . The explanation i is that in straight 
channels there is a double spiral flow as show! n in Fig. 16, , superimposed o on the 


forward anes which carries the slower moving water from the bottom up along 
the sides and out along the top to the center. ‘At the center it mixes with what - 
would otherwise be the fastest m moving water, thus slowing the rapid water to a 
velocity somewhat less than that in the region just below. ‘name unaeil 


_ Those \ W ho i inquire as to what causes these spirals ae a very amusin : 


situation. A. follows 0 out thinking « of Mr. Stearns and state a 
Ry rietional losses at , the sides reduce the energy and thus the head of 
the water in their neighborhood, with the result that the surface level at the 
: sides is lower than near the center of the stream, and the cross- -sectional 
profile of the water surface is a curve concave to the bed. " Owing to — 
a superelevation of the water near the center and its tendency to find its own - 
level, transverse currents are set up which travel downwards near the center 
of the stream; outwards along the bottom to either bank, upwards along the 
sides, and, for permanence of régimé, inwards along an and near the surface.” 


Mr Moller, 38 however, takes just the opposite view: 


“Each filament of water that ‘approaches the bank * * * suffers re- 
tardation, and hence requires greater cross section, whereas the conditions 
in the center of the stream are reversed, water filaments there being ac- 7 
. celerated and requiring less cross section. Both actions lead to the produe- : 

tion of a slight transverse drop in eng water surface * * * and in the | 


Pe formation of of of spiral on either side.” 


Mr. Miller then shows an iHustration with the water surface highest at the 
sides; but all three agree the are as. shows n in Fig. 16. 


37**On the Current-Meter, Together a Reason Why .he Maximum Mag of Water Flowing in 
3 Ch: annels Is Below the Surface,” by F. F. P. Stearns, Transactions, Am. Soc. C. E., Vol. XII (1883), pp. 
-'8 “Studien tiber die Bewegung des Wassers in Fliissen mit Bezugnahme auf die Ausbildung des 
Flussprofiles,”” by Max Moller, Zeitschrift fiir Bauwesen, 1883, p. 201. A summary in English is given in 
Hydraulic Laboratory Practice, edited by J. R. Freeman, A.S.M.E. 1929, p. 70. 
$9 **Hydraulics and Its Applications,” by A. H. Gibson, D. Van Nostrand Co. 3d Ed., 1925, p. 21, 
following the article ‘‘On the Depression of the Filament of Maximum Velocity in a Stream Flowing T rough 
an Open Channel,’ "by A. H. Gibson, eRe, ae Soc. of nee Series A, wa 82, 1909, pp. 149-159. 
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at the ieaiiiaes sala: but, since the velocity of the spirals is very small, the cor 
additional head due to curvature would be . extremely minute, and at any rate 
‘it would be an effect, rather than a a cause, of the Spirals. — The writer has the 
same feeling i in regard to the author’s ‘eteaiedien (10), and believes that the 
double spiral flow is a cause of the unequal distribution of suspended load, 


rather than a result o of it. 


es cause of these double spirals is still obscure to the writer. Prandtl 
- has shown why there will be secondary flow toward the corners and away from 
the sides 1 in n the case of turbulent flow in rectangular pipes. — Similar reasoning 
"applied toa rectangular chann el might explain spirals somewhat as shown i in 
| ‘Fig. WV, _ but this is a pure flight of fancy on the part of the writer and he has 
> no sonal that it actually exists. a ‘There is just the possibility that such spirals 
Ton have been observed by Messrs. ‘Stearns, Miller, and. Gibson, , and have 


been mistaken for sais form shown i in Fig. ie _— it does not seem very probable. 


Average Mean Ratio of equiv: alent 

of values ‘concentration, diameter, roughness to 

7 (mm) of bottom sand 
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The other point. whieh. the writer wishes to discuss is | the ‘reduction i in 1 re- 


sistance to flow whicl hich author fi finds accompanies the presence of sediment. 


(26a) 


n which c is the Chézy coefficient, R is the hydraulic radius, ¢ is the equivalent — 
ee of | the | channel bottom (in the same units as R), and ¢, is the hel 


4 Uber die “‘Turbulenz, by L. Prandtl, Proceedings, 2d International Cong. for Applied 
Mechanics, Zurich, 1927, p. 62. The essence is given in English i in ‘‘Fluid Mechanics for Hydraulic Engi- 
neers,”’ by Hunter Rouse, McGraw-Hill Book Co., Inc., New York and London, 1938, pp. 267-268. 
_41‘*Flow in a Channel of Definite Roughness, by R Ralph Ww. . Powell, Am. Soc. C. _E. 
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correction coefficient, proposed by G. i Keulegan® to correct the Chézy 
“coefficient to a corresponding channel of infinite width; that is, to eliminate the 

effect of the side-walls. (The value of c, is 14.18 times the value of B given by 

Mr. Keulegan. ) I In the author’s experiments the side-w alls were smooth and 


far away f the s sections studied, sO that Cs may be as 


ralues of w ere computed for all the data of Table. and ‘from these the 
value of « €, the « equivalent roughness, § given by Eq. 26b, was determined. These 
values w vere then divided by the diameter (in feet) of the sand forming the 
bottom roughness. For water without sediment it might be expected that this 
ratio would be 1, but the scale of roughness used in Eqs. 26 is not quite the same 
as that of Nikuradse, values of the writer's s e being. somewhat smaller than 
- Nikuradse’ s sand size. All r runs with approximately the same concentration 
were grouped together and the ratios of € to size of bottom sand were averaged. 
wean 5 shows the results, arranged in ascending order of sediment size. The 


general y W veighted average for each size of sediment, irrespective of concentration, 


Ratio of equivalent 
wi. 1, millimeters 
60.32 
0.42 


Although it that the roughness i is reduced by 
the presence | of sediment, there is great variation in the effect and there are even 
a few cases in which there seem to be an increase (runs 3, 5, and 7, Table 5). 
z= also seems to be no relationship between the amount of sediment and the “ 
—_— run 17 with only 0.048 g per liter giving the largest reduction of any 
concentration of the 0.160 mm sediment. _ The only consistent variation is that. 


| shown by the a average for each size, which indicates that the reduction is greater 


for the smallest size and least for the largest. _ This might be interpreted as _ 

- meaning that the reduction in resistance isa function of the number of particles — 

present rather than of their total mass. It should perhaps be pointed out that 

the variation is not quite as great as might | be thought at first. if the bottom 

sand 1 were 0.47 | mm in diameter and the depth 0.50 ft in both « cases, , the v: values" 
- of c corresponding to the ratios of 0.34 and 1.07 (the smallest and largest — 


in Table 5 5 for 0.160 mm eee are 119.2 and 9 99. 3, &\ variation of of only 9.1% 


- The writer has not tried to touch on some of the most interesting points in | 
this valuable paper, which can better be discussed by those more familiar with sl 


> details of this subject, , but has tried to indicate s a few v directions i in which he 
hopes that further investigation will be 


—42*Laws of Turbulent Flow in Open by Garbis H. of Research, 
Bureau of Standards, Vol. 21 (December, 9. 
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E. R. Van Driest, 43 Assoc. M. C. E. *—Experimental investi- 


water where the of is great. paper 
t ca the trials and tribulations of an experimenter (1) i in building apparatus 7 

- for the proper c control of conditions and (2) in bringing theory ¢ and data together. " 
_ It In the presentation of the theory, it it might add ‘to the paper to point « out, ——? 
ei fact that the problem of suspension of sediment by flowi ng water is but. one «| chat 


of the three or turbulence transfer problems usually encountered in = dat 
( © 


attention is Fun heat transfer, the formula for w hich i is 


in ia in addition to the notation of the paper, q is the rate of heat transfer 
_ acToss a unit a area normal to the y- axis; ch is the specific heat of the fluid; T i is 
the temperature; and is the heat transfer coefficient. Through thea assumption 
that Ga = Ex m (the basic assumption of the analogy. between heat 
; transfer and skin friction), great progress has been made in heat transfer for 
“flow in pipes. - The third transfer problem—that of sediment transfer—has for 
its basic relation Eq. 4, _ The paper proceeds to show that progress is now being 
made in the sedimentation field when e, is assumed. equal tO €m Or when the — 
mechanism for the transfer of sediment is assumed the same as the mechanism _ 
for the transfer | of momentum. Since these three phenomena are of the same 


nature, they ¢ can be derived from | one general transfer formula, — 


. in which Qs is the mean rate of transfer of any property @ per unit volume across 
. # a unit area perpendicular to toy. By merely substituting any of the foregoing - 

properties per unit volume—that is, momentum pu, heat content pc, T, or 

sediment concentration C—into Eq. 28, the aforementioned equations will 


Indeed, any other properties may be assumed to be transferable. in 


‘The author shows some concern about the fact that the velocity distribution 


2 data do not fall entirely. on a straight line when plotted on semilog garithmic 
paper. _ Furthermore, he states that this departure is due to viscous action, 
but finds. that such an assumption does not explain the results. The truth of 
the matter is that Nikuradse’s data do not fall exactly on a straight line 

either. either do his. data for smooth pipes® fall precisely on a tne 


Asst. Prof., Mech. Eng., Mass. Inst. Tech., Cambridge, Mass. 
Received by the Secretary November 4, 1944.00 


Rag yg in Rauhen Rohren,” by J. Nikuradse, Forachungsheft 361, Vereines Deutscher 
Ingenieure, Vol. 4, 1933, Figs. 15a, 15b, 15c, 15d, 15e, and 15f, p 9 


8 “Clesetemissigkeiten der turbulenten Strémung in Glatten Rohren,” by J. Nikuradse, Forschungsheft - 


6, Vereines Deutscher Ingenieure, Vol. 3, 1932, pp. 1-36. 
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when one plots the data sepa ion for each run. All Nikuradse’s s data for | 


smooth and rough pipes, except those for the lowest three Reynolds numbers 
in the smooth pipes when the - viscous layer does show its presence, have a — 


characteristic S-shape. line plot does appear toe exist when all the 


are on ~ log = \-scales, in*® which ¢ is the linear dimension 


= 
g roughness; ever, even here the prevailing S-pattern | can be 


observed for each run. That the viscous layer does not affect the author’ Ss 
data is evidenced by a simple computation w hich | Will show: that the position — 


of the d data taken nearest 1 the bottom — is at least ten times as far from the 
‘bottom of the roughness grains as the edge of the laminar boundary layer. 
which is at the most one half as thick as the sand Toughness grains. = _ Measuring " 
the distance y y from the bottom of the grains makes the measurements consistent — 7 
with those of Nikuradse who measured his distances from the wall or, rather, 7 
ig from the base of the : sand grains, and this ¢ does have a decided effect upon the 7 


plotting of the data v very close to the wall. In Nikuradse’s plots for rindividual 


runs, the straight line drawn for k =0. 40. passes through the data as a mean a 


line. It appears that t such a line for k = 0. -40 has the same mean elation to 


the data of Fig. 8 for no sediment added, even after the data are shifted to 
conform y vith measurements from the bottom of the grains. 
_ The data of Fig. 8 seem to prove quite conclusively 1 that the value « of k i 


decreased by the presence » of sediment, and it appears logical that this should — 
be so, because of the dampening effect of the sediment. _ - However, the absolute 
effects of size and quantity are still uncertain, thus showing the need for: further 
controlled experimentation. 7 If the concentration has a decided effect 1 upon nk, 7 
then it. should be logical to assume that the velocity distribution curve should : 
be distorted considerably because of the large variation of concentration over — 
the depth. . This effect does not appear to exist. 
An intriguing observation i is that the exponent 21 suddenly changes from the Oo 
chee 
theoretical value when the : size, and not the amount of ‘sediment, i is changed. : 
The only explanation that the writer can give for this disagreement follows - 
from the apparent circulation o of sediment. ‘Since the value ¢ of € in any vertical a 
line is more uniform near the wall than in the central part of the channel, the 
sediment transfer theory requires the sediment concentration to be more uni- 
form: vertically near the wall th: than farther in the channel thereby causing the 
concentration to be greater at the water surface near ' the wall than el 
channel. — ‘This variation along the surface then causes a a horizontal flow of 
sediment toward the center of the channel. In this process, however, all the 
sediment at the bottom corner is : picked up, leaving : a sediment gradient he 
the bottom decreasing toward the corner, which in turn results in a flow of 
sediment along the bottom toward the ‘corner. The « entire process pa 
results in a circulation of sediment, not necessarily water, which may develop : 
into a complex circulation pattern. | It follows that heavier or larger a 
“should not be diffused too far into the flow and therefore that such sediment ;, 


in Rauhen Rohren,” by J. Nikuradse, Forschungsheft 361, Vereines 
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should be more uniformly distributed horizontally i in n the central part of the 
_ channel. — On the other hand, the e lighter or ‘smaller grains w ould be diffused 
farther into the flow, thus easily ¢ causing higher concentrations near the surface 


the ‘theoretical \ values. It is ‘concluded, ‘therefore, that 
: traverses throughout ° the stream are necess sary to interpret, fully, the results 
4 As mentioned by the author, sediment appears at the stream surface eyen 
though the momentum transfer coefficient there must become zero if the 
logarithmic v elocity distribution law is assumed to hold all the way ay ‘jo the 
surface. This s inconsistency promptly | disappears when it is noted that the 
logarithmic law does not hold at the surface . the velocity gradient there is to 
be x zero because of a presumably zero shear. 4 Although « one would expect. the 
transfer coefficient to approach zero. at the surface as a result of a zero 
approaching velocity gradient, yet some agitation should be expected because 
of the migration of eddies from the adjacent lower layer of turbulent fluid. 


rivers, _A more simple case of such di diffusion i is demonstrated d by the movement 

of eddies from a ‘point of disturbance at the bottom of a jar r of water into the 

overlying quiet fluid, 


; = writer wishes to re-emphasize the need and importance of such basic 


experimental work in sediment transportation. The ‘Peper 
considerable progress in this field. 
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Evidence of this diffusion can be observed readily on the surface of turbulent J 
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By G. GREULICH, AND RAYMOND P. PENNOYER 


M. Soc. C. E. new ‘methods, and an 


pansion of some old ones, for the precise mathematical evaluation of cortal ain 
phenomena relating to cofferdam design are offered in this paper. — The phe-. 
nomena have been noted before, and empirical allowances have been made for 
their design and construction. The new concepts of analysis presented in the 
paper should be useful in ‘providing additional cross checks for the adequacy 
of a cofferdam design. However, the use of precise mathematical evaluations | 
can result in misleading, erroneous, and dangerous conclusions in the hands of 
_ inexperienced individuals, w hen dealing with such uncertain quantities as the 
widely” divergent ‘materials in ordinary soils whose properties resulting 
- physical characteristics and effects can vary widely due to either man-made or 
Deflection of Cellular Cofferdams.— —In cellular cofferdams, the nature of “oy 
deflection depends largely not only | on the elastic properties of the fill in the 
cells, but also on the method of depositing the fill. 4 For example, fill deposited 
in cells inundated i inw ater will not produce the immediate tamping or Seeiel 
effect obtained by ordinary methods of dry filling. Tightening the interlocks 
removing loose play, as well as by stretching the. entire cell, aids materially 
in reducing the deflection. The analogy is similar to the comparison of | an n 
empty headless barrel or keg | with dried- out staves and loose hoops to one 
with tight hoops a and staves in solid contact. mi ‘Even without fill, the tight barrel 
_ will stand considerable lateral pressure at the top without appreciable deflec- 
— tion when the bottom is fixed. The loose barrel, of course, will rack far out — 
‘ of position. Where cells have been built and filled i in deep water, the deflec- 
tions under load are less for subsequent unwaterings than for the first unwater-_ ae ~ 
ing because the consolidation and settling effect from the first drainage of the ve 
_ fill in the cells is still felt during later cycles of filling and emptying the coffer- ns da 
-Nors.—This paper by Karl Terzaghi was published in September, 1944, Proceedings. as 
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dam. The deflection is also affected appreciably by the character of backfill. 
material and the method of its deposition. 
Permeability Cellular Cofferdams and of Their Sheet-Pile Enclosures 
- | With reference to to the permeability of sheet piling n near the b base of cells, tension 
_ in interlocks also y prevails for a short distance below the base of fill. ~The: softer. 7 
= the foundation material, the greater the depth of tension will be 
Resistance Against | Siding; - Cellular Cofferdams on Rock.—Sheet piling is 4 
: “usually “toed in” on rock by driving 80" that it penetrates the soft upper layers 
“of the rock. ~ here the rock i is too hard to be penetrated by the sheet piling, 
mechanical methods such as dowels grouted into the rock have been proposed 


: to prevent sliding. The method of grouting dowels into rock is highly objec- 
tionable due to ‘intermittent ‘supports and thigh stress concentrations | at the 
- points of s support which w ould tend t to rupture the sheet- pile interlocks. : 

Another ‘proposed method, which presents some very difficult construction 
problems, is the concreting of s sheet piles i in n trenches cut into the the rock. Bio. ks 
Prior to the early 1930’ s, cofferdam design methods were based largely on 
“rule on precedent. The customarily accepted Pennoyer 
(10)(11)* > was offered as a rationalized, ‘simplified procedure, and 
as a perfected, all-inclusive scientific approach to the problem. 
— AS Professor Terzaghi states (see text concerning Eqs. 1 and 2), the Pens 
~ ¢ noyer design . method disregards. some factors that 1 may reduce the margin of 
stability. When these factors are omitted, greater sliding resistance is indi- 
cated than is really the case when the ; ostven! shear distribution and height of 
‘saturation are taken into consideration. _ Inasmuch as there seem to be no 
‘sliding failures on record, the Pennoyer method has : apparently been adequate, 
even though the factor of safety was not as high as computed. — _ A more precise 
approximation ean easily be obtained by applying the author’s method. 
- However, as Professor Terazghi states (see paragraph introducing Fig. 3), 
there is at present no reliable method of determining i in advance the state and 
properties of the soil in the fill; and quantitative data on the ‘stiffening effect 
_ of the sheet-pile v walls are not: available. Therefore, during the present state 

of the art precise mathematical ‘analy best be considered only 
dq Resistance Against Internal Shear; Cellular Cofferdams on Rock.—The author 
_is in error in stating (see paragraph introducing Eq. 9) that possible failure by 
- ames along the vertical neutral plane in the cells has never received any atten- 
Colonel Pennoyer has called attention to the investigation of internal 
_ shear) resistance and its relation to the safety of the cofferdam (10a)(10b). 
Resistance Against Overturning; Cellular Cofferdams o on Rock. —The pre 
ceding remarks on resistance against sliding are also applicable to ‘the f factor 

of safety : against overturning. _ Since the factor of safety of 3.8 suggested by 
the method of Colonel Pennoyer wh hen applied to the selected example (see 
text following Eq. 20) is based on premises that may not , actually exist, the 
factor with the assumed values, as asserted by Professor Terzaghi, is still 
about 2, which should be adequate. The writer presumes that in making ing this 


4> Numerals in parentheses, thus: (10), refer to corresponding items in the Bibiography (see Appendix 


I of the paper), and at the end “ discussion in this issue. _ 
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assertion Professor Terzaghi 


shear are largely offset by the ieaetian: of soil compaction . and of friction in the _ 
* 7 ith respect to overturning, fortunately the e cellular type of cofferdam can 
be widened at ‘practically. no additional cost. Extension of the width makes 
E the spacing greater between cross-walls so that the amount of steel is not 
— increased in circular cells and is increased very little on the straight diaphragm 
_ wall type. In any case, the extra fill is a minor item. Of course, the allow- 
‘ a able interlock tension may limit the width for a given height. In general, the 
_ cofferdam should be made as wide a: as Teasonable limits ts allow. This i — 
stability and also helps i in reducing p permeability. 
aii. he ‘increased stiffness of cellular cofferdams, as —_— with other types, 


; has been — in many ways, and lock friction (see second paragraph 


“cofferdam soon to increase the sliding resistance of one piece 
- piling against the other in 1 the cross-walls and that a large percentage of the | 
to shear may developed by the steel the earth fill is 
'There- 
fore, true may lie somewhere the value for earth fill alone _ 
4 


and the combined value, as determined in the paper. _ Professor Terzaghi has: 
proposed a method of approximating the combined values which appears to be 
. simpler and more ‘rational than any that has previously | been brought to the — 
writer’ s attention. The author’s method, however, , does take care of 
another factor introduced by the sheet piling. _ The added factor is is the direct 
tension introduced | in | the: sheet piling when the cells cant or lean inward. 
‘This movement crowds the interlocks together so that individual pieces can 


exert tremendous resistance to sliding or upward movement. Because of these 
phenomena, the diaphragm | or cross parts of steel sheet- walls 


down to the outside bottom o: 


_ the introduction of stay rods from the i inner top to the outer base of the cell. 


= fact that interlock tension failures, a as s observed by y Professor Terzaghi, 


preceding Eq. 23) se seems to indicate that pressures sand shear 

i resistance in the fill are considerably more than the Rankine active pressure. 

: This is probably due to the tamping effect of dropping fill in cells, as well as 
to. a squeezing and consolidating effect occurring as the partly filled cells weave 

b- stretch during filling and drainage of water. These factors all eeeeare oil 

to the ultimate stiffness and stability of the -cofferdam. sf 

Ene There could be considerable danger of interlocks failing by tension at aan 
inner toe of a cofferdam on rock if the piling were not. toed or bedded in the — 
tock. Only the most foolhardy individual would permit such construction, 


especially cells. is is generally recognized that driving on rock tends 


1 
— _ the stiffness of the cells. Colonel Pennoyer has mentioned this fact (100), et 
n — Heretofore, it has been difficult to evaluate the combined effect of the shear rT a 
resistance of the fill and of the steel interlocks. It is recognized (returning to 
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damage the lower ends of the 1e sheet piling an may produce incipient 
men 

locking of adjacent piles. For this reason, manufacturers always recommend, 
“ and sound practice should dictate, the secure embedment of the lower ends on 
Even n if the vertical reaction of the piles from overturning moment is 
theor etically y sufficient to hold the: piles in place as is repeatedly demonstrated 
_ by the shape of the bulge above the inner toe, the remarks regarding embedment_ 
_ of the lower ends still apply. Although, as Professor Terzaghi has emphasized _ 
(see third paragraph following Eq. 27), two of the three customary equations” 
may be inadequate and the third is open to question, the fact remains that com- = 
paratively little trouble has been experienced with cellular cofferdams. As the 
a author states, in the application of these > equations, very conservative values" 
have been used for soil constants. i He ex expresses concern 1 that an ambitious — 
2 : designer s some day might use soil constants determined by laboratory tests, with ; 
fatal results. Fortunately, -cofferdams usually represent the e ‘expenditure of 


considerable sums, and the final decisions regarding design and construction 
‘jon been in the hende of mature, experienced individuals. The writer knows 
of no serious failures in which customary practices were followed by engineers, — 
who generally allowed ample margins beyond the obtained by 
equations. 

ae more exact methods of design are desirable to to prevent future 2 economic 
waste. _ There are also so cases in which | cofferdams with dimensions : smaller 
4 than those indicated by the e usual equations have been used without. encounter-_ 


a ing ‘difficulties; s; although i in a few of such instances the factors of safety w ere 


ae low that disaster could have followed any minor disturbance of soil or water 


a level. In general, however, unanticipated forces or acts of nature have ‘caused | : 
disasters which conceivably could have regardless of any kind of 
Penetration of Sheet Piling; Cellular on Deep Sand Strata.— 
"Aside from consideration of uplift of sand due to percolation, it is believed 
— the depth for driving the sheet piling in sand need not be 3 H (see second | 
* paragraph following Eq. . 29b) but that 3 H would be -ample—because the sand 
‘in the cells can grip ir one side of the sities, , whereas, 1 when. saaane the piling 
“Sand Strata.—Sheet piles in pairs may be driven and then loaded or extracted 
by means of hydraulic jacks to determine the forces required t to shove down or 
a ae pull out the piles: (see second paragraph following Eq. 37). In similar soils, 
the values for a , given depth of penetration are , approximately uniform. — It is 
easier ‘to make an extraction or pulling test than a load test on the sheet piling, 
a8 anchor piles are not needed and the ‘Tequired pull may be exerted 


value might be much greater ‘then extraction value. The ev: cer 
of vertical load capacities on sheet piling by means of hammer formulas is a 
very uncertain practice. The resistance of the interlocks when driving in a 
is a greatly varying factor— —s0 great in some instances as to prevent 
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_ driving t to moderate depths. 4 ‘The values for the sheet piles i in . the soil must be a 
determined on piles. separate from a continuous wall. 
- Conclusions.—The conclusions are most excellently stated. However, _ * 
clusion 2 needs further clarification and possible expansion. With respect to 
con conclusion 5, the reduced pressures seem to be due to arching or or cohesive action, _ 
or to a combination of both, since the Rankine theory is conservative and 
deals with granular materials only. In conclusion 6, the depth of ee 


3 sheet piling can be dictated by requirements of watertightness which, in af 


> 


some cases, n may indicate a a greater depth than i is needed for load capacity. 
Cellular Versus Double-J Wall Cofferdams.— 

‘the superiority of the cellular, as compared to the double-wall, cofferdam are 

~ true only when full consideration has been given to adequate foundation sup- - 
‘ port and the reduction of permeability i in the supporting soil. Inasmuch as 
7 ~ double-wall cofferdams usually have inner, and often outer, berms, their service 
7 record has been about as good as that of the cellular type. There have been 
i failures of individual cells due to boils or undermining from percolation which 
not have occurred on a wide-bermed, double-wall installation, 


. He has provided | 


- With this paper at hand, the designer may readily differentiate between an 

ultraconservative design for a per 

a temporary | structure | which, “although it must be safe, may be allowed to 

excessive ina ‘temporary structure, it may be checked by one of the 

_ The author i is to be commended fc for the thorough, thoughtful, and orderly 

"presentation n of a comparison between old and entirely new ‘mathematical 
"approaches on this subject. He deserves s the w wholehearted thanks of th the engi- 

_ iheering profession for assembling into convenient reference form much valuable 
_ historical and technical data on cofferdams, as well as for providing a most 
complete Bibliography on this difficult, complicated, and abstruse —— 

about which considerable quantitative information and experience must still 


_ Acknowledgment i is andes to J. R. Woods, sales engineer, Carnegie-Illinois 


Steel | Corporation i in Pittsburgh, Pa. , for a constructive review of this discussion a 7 


RAYMOND P. Pennoyer, Assoc. M. Am. Soc. C. E. “—The basic computa-_ 

A tion to determine the necessary width of a cellular cofferdam, by bringing the _ 
_ intersection of of the resultant of the moments created by two opposing forces 


= lateral overturning force and the weight of the contents of the cells) a 
within the middle third of the base, is expressed by Eq. 8. The author has — 
stated that “This conclusion is open to serious objections.” The only : substi-_ 

- tute offered in the paper for this basic determination is a statement that ie 


— 


5 Col., Corps of Military Police, U. S. Army (Inactive), Pittsburgh, Pa. 
Received by the Secretary Ni November 13, 1944. 
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PENNOYER ON CELLULAR COFFERDAMS Discussions 

‘average width” of circular cells built on a rock foundation has | been 0. 85 H, 7 

in the past (see paragraph containing Eq. 27). All the remainder of the paper 
1” analyzes the theoretical internal stresses in a cell on the basis of this ratio. 
_ Overturning Moment; Width of Cells. —The writer believes that 0. 85 H has 
been more nearly the minimum width of completely excavated and unwatered 
cellular cofferdams than the average value. Furthermore, the proper and Safe 
width of cellular cofferdams from which Professor Terseghi has derived this: 
rane average has been determined, for years, by the very equation he 


it would seem axiomatic that, since the weight | of the material contained in — 


a cell determines its capability to resist the overturning moment, this material 

_ should receive . careful consideration. Broken stone, gravel, sand, sand and 
—_— etc. , have different weights, as have the undisturbed soils in which cells 
are built tow various elevations. | ‘The analyses i in the pa paper have all been based 
on n sand Ww eighing 11 110 Ib ‘per cu cu ft the cell. No mention been 


i 


weight of the contained soil (that is is, the submerged weight) ; nor does 
- Professor Ter rzaghi explain how he proposes to adjust the width of the cells to 
make the cells safe against overturning when ane are almost infinitely variable 
factors affecting the weight of the contents. 
2 a No method is proposed to determine the width of cell necessary to resist the 
- lateral or overturning forces safely. | The paper deals only with an analysis of a 
_cofferdam built in water on bare rock, with no submerged soil overburden. This 
— is extremely rare in practice and alw always requires that sheet piles be. 


ican 


driven into empty cells on the rock—a major or construction p pr roblem. The w width 
of a cofferdam in salt water should be different from that i in fresh water, even 
“under such conditions. — Likewi ise, there 1 must be an adjustment i in width to 
resist the lateral pressures of submerged soils, for various elevations, in combi- 
nation with layers of natural soil or fill above such submerged soil. There must 


also be an adjustment to resist other possible lateral forces due to construction 


- equipment, storage loads, railroads, derricks, | Wave action, ice, swiftly flowing | | 


water, etc. This basic. consideration has been. dismissed in the } paper by stating 
that the average | width of cellular cofferdams in the past has been 0.85 A, and 


further that this ratio may be be considered an acceptable basis for the duden of 
cofferdams i in clear water on rock. 


7, with its ‘supporting text, is the assumption that a 
athe cofferdam filled with soil has a factor of safety of 3 when the width is 
alculated to bring the resultant of the forces within the middle third of the 
ase. - This formula is in serious error since the soil within the cells must be 

assumed incapable of Tesisting any tension or uplift along the outside face (the 
heel) | of the cofferdam. e _ The width of a cell must be adjusted to bring the stress 
to zero at this point: (neither tension 1 nor ‘ compression) anc and, as a consequence, 
- the factor of safety against overturning is s actually unity, which means | that no 


factor of safety exists. Professor Terzaghi then discusses the possible, theo- 


‘retical distribution of | pressures under | the cells (Fig. 7) and decides that, since 
a sand- cofferdam having a of 0.85 actually has its 
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Internal Shear—An attempt was made in 1934 to mathematically 
as ~ occasional leaning or distortion of the cells because | of evident lack of internal el 
ed &f shear resistance of the ‘soil contained in the. cells’ (10). At this time, , the 
offered the following formula (10c) f for for vertical and horizontal shear, s, at 


9, giving ‘the total along 
so-called neutral plane of the cell here a given moment, results in exactly the same _ 
value as Eq. 39 for maximum shear. Several successful bulkheads have been - 
built which, upon analysis, were not more than 60% safe in internal shear, but 
which, at the same time, had a width sufficient to resist the estimated overturn- 


ing based on the | “middle third” method of determination. A 1 mathematical 
| te of analysis will never ‘serve as more than a guide based on ‘precedent, 7 


because the assumption that the soil fails, along a a plane of maximum shear, 

the vertical axis, similar to al beam, seems fallacious. 


alk W hat actually happens is that the entire mass of nated soil changes" 
- shape simultaneously. This action can best be described as a universal plastic — 
_ flow, or distortion, , which bears no relation to the strains : set up in an elastic 
homogeneous mass that can resist shear, compression, and tension—w ith 
‘Special e emphasis on tension. Ww hen a cell leans inward, as a result of the simul- 
taneous distortion, or flow, of all the contained material, the sheet piles i in the — 
intemal, or diaphragm walls, remain parallel: to each other and necessarily 
maintain their original length to lean into a parallelogram. | ‘Under this condi- 


: tion, all the interlocks in a » cross-wall, _ or interior wall, slip the same amount at 


~ the same time, irrespective of their position in relation to the fictitious neutral 


cae .s the writer stated previously ( (10), the r resistance to to lateral deformation of of 

the bulkhead, through friction in the interlocks of the sheet piles, a affords an 
7 indeterminate factor of safety i in the ¢ calculation of internal shear. This still 
seems true despite the author’s able mathematical analysis of this feature for a — 

hypothetical cofferdam with a width of 0.85 H . With | sufficient sand trapped | 
and | wedged into the interlocks, or sufficient penetration to develop frictional — oe 
against the sliding, ‘the shear resistance of a cross-wall might be 


equal to that of a continuous plate of steel over the entire cross-wall. «dit is 


likewise possible that the assumption of the frictional value of an interlock as . 
oe tenths ‘of of the actual tension may be entirely too high. A cellular bulk- 7 
head built in undisturbed virgin soil, followed by excavation ot on one side, isa _ 

case in point. There Tight be no tension on the interlocks before the cells 3 
— Furthermore, the paper states (see sentence following Eq. 21) that the fric- 
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interlock—wit 


= 


pore and ‘extracting, with emphasis. on n reducing the area of contact “olin 
surfaces to a minimum and on eliminating wedging action from forces tending | 
to separate the interlocks. There is no relation between the frictional value of : 
two pieces of steel in what amounts to line contact and the coefficient of friction 

a of s steel on steel times a tension force per unit of length, when no area of contact 

ean be assumed. In the discussion of the bulkhead in Fie. 11, Professor Ter- 
- gaghi has ‘expressed some doubt : as to the safety of using three tenths of the 


tension for the frictional value o of 2 an interlock in his : analysis. of the actual 


7 ‘Tension in Interlocks— —The author recommends that the e method of 
Te proximating the possible maximum interlock tension by the so-called Rankine 


ormula, u using ta i arrive at the conversion factor between the 


_vertical and horizontal pressure, ‘be abandoned in favor of using 0.4 when the 


cells are filled with sand and gravel. — He states that this method does - not 
create the illusion of resting ona scientific basis. % It i is seldom, however, that a 


entirely filled with sand i is encountered i in actual practice. 


al ‘Professor Terzaghi indicates that designers seldom, if ever, use a conversion 


factor of much less than 0.4 for the lateral pressure of the contents of a cell 
_ filled with sand or gravel, completely drained and founded on rock. | With a 
_ working stress of 4 000. Ib per lin i in. in tension for the interlock of the sheet t pile, 
- cellular cofferdams w vould be limited to a height of approximately 48 ft and a 


maximum radius of bret of about (24 ft. a With a working stress of 6,000 ‘lb per 


beyond these containing soils far less 
ideal than the assumed sand and gravel, this conversion factor of 0.4 does not 
agree | with actual practice. A conversion factor as low as 0.17 has occasionally | 
been used without any indication that the lateral pressures in the cells were 

underestimated. _ ” Obviously, the ‘relieving effect of wall friction, even to the 
point where a complete arch is formed as described subsequently, is an im- 


portant factor to be considered when establishing the conversion factor for use 


_ The paper offers no , suggestions for approximating the possible interlock 

tension: (1) When the cells are ‘partly. full of water and submerged soil (see. 

“Part Practical Considerations’’) ; (2) when the cells have a layer of clay 

oe: or other soils at the bottom or at various heights as part of the original s soil on 


which the bulkhead i is : built; or - (3) ‘when possible ‘surcharge loads are present. 


— 
this factor was assumed ana ogous to the 
“8 ie of friction of steel on steel, the statement would seem fallacious. Interlocks of &§ = 
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condition not discussed in the paper is that in which 

~ contained soil is not fill but i is composed of various types of undisturbed ‘soil 
into which the bulkhead is built before excavation on one side proceeds. — i 
In his prior paper (10d) (10e), the writer illustrated | a ease which ex- 
ternal and internal earth } pressures ; differ which included \ weep holes, es, submerged 
soil, and water levels at various heights. In the light of the practical problems 

ill illustrated and discussed, and considering the complexity of conditions 

countered in practice, the writer feels that Professor Terzaghi has oversimpli- 
- fied t the problem by assuming a cell width of 0.85 H and a conversion factor | of . 


hey Probably no one realizes more clearly than t the writer that, at best, the use 


of the Rankine conversion factor is a very - approximate assumption. . He will | 
continue to use it, however, until something better is available, especially si since 
water r and submerged soil seldom are the critical considerations i in overturning 
forces, weights of material, and bursting pressures. In the light of many 

precedents, , the use of this is a guide w which helps the engi- 

bro 

hood will be ‘completed. on the other hand. 
_ Mathematics is is only one consideration in reaching this decision, because the 
y verdict must take into account the strength of the sheet piling as manu- 
-factured, precedent, the variety of conditions into which the sheet piling is set up — 
and driven, and the experience, ability, and equipment of the constructer. 
‘anes can be ruptured by sand trapped in the bowls during the pile « driving; 


builders can fail to set and drive the sheet piling accurately; or the final closure ~ 


pile ¢ can be used to pull the remainder of the wall into its proper sr place. _ Ing all 
such cases, indeterminate stresses are created in the steel, or an interlock is 
“ruptured i in the process. Certainly, mathematics is futile under such conditions. 
‘Sand and gravel, or rock, may actually arch and form a complete bottom i in 
i’ above the base of the sheet piling. About 1932 the writer had occasion to : 
i¥ Basen a bulkhead composed of a series of diaphragm-type cells, about 35 ft by 
27 ft between cross-walls and 35 ft high, the sheet piling being driven into sand 7 
and gravel along a river bank. - With the bulkhead in place, the sand and > 
gravel on the river side were excavated to bare rock and used to fill the cells 7 
— to the top. Me As excavation proceeded, it was found that the sheet piling had cm 
been through: a layer boulders to about from and 


Water ater coming through the bulkhead had washed away the original soil theough 
the gay gaps in the e piling up to the top of the breaks, some 12 ft off the rock, leaving | 
_ three or four of the cells completely | empty to about. that. elevation. The ite 
ry were standing, however, and were resisting the maximum lateral load from — 


the river bank as intended. _ Above the 12-ft elevation, the sand and gravel 
had formed into about as perfect a dome arched between the side-walls, front- 


walls, and as could have been built in 2 brick 0 or stone. 
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days later, that a bastet actually had to be used to break the dome , So that the - 
An analogy may be found ina pail filled with water. fi Its resistance to ov over- 
nnine depends upon its w eight, height, and cross section. ~ All shear i is taken 
by the steel shell of the pail as water has no shear or frictional resistance » value. 
Other ‘material, such as $ sand and gravel, does | have a high shear and frictional 
value. and consequently | possesses : the inherent : ability to arch as described in the 
aforementioned example. _ Thus, the shear stresses in a steel bulkhead would 


ta be relieved by the soil or fill material, whereas in the pail of water the tots il 
* shear would have to be transmitted by the walls of the pail without relief. 
Pre... the foregoing theory is correct in principle, it would indicate that ‘per haps 
the “ ‘middle third” limitations on width are somnetvalbee and do have a factor 
of safety k because the “ bottom” i is actually able to resist some uplift at the heel. 
a This may be one reason w hy y mathematical analyses « of internal shear seem at 
_ wide variance with the actual resistance of cells to leaning or distortion in some 
eofferdams. . Calculations of interlock tension of moist (not submerged) soil 
_ thus become a study in the theory of pressures in bins. _ This may ace 
count for the fact that the maximum pressure § seems to be a at the zone of g greatest 
in a cell, ordinarily about halfwa ay up the height, although it has been 
observed t that this varies § slightly for different soils. — _ The maximum bi bulge may 
 oceur at the level of the first actual “bottom. ” Of course, such a theory cannot 
be used when there i is water in the cells but ‘the bursting pressure from the ‘soil: 
alone might safely be ignored under ideal conditions below the point of maxi- 
mum b bulge. It is difficult to o imagine a. a condition for which the pressure of the 
water in the soil at the bottom could be critical, especially if the bulkhead had 
wi Furthermore, if the foregoing theory i is correct, , the distribution of pressures: 7 
at the surface of the rock | (Fig. 7) would differ from the true conditions, because 
alll the pressure at the toe would have been transmitted to the steel by fi friction 
and arching. 7 ‘Observers have long noted that the overturning moment does. 
not create any consequential added bursting | pressure in the interlocks of the 
"sheet piling because no consequential | added deflection of the steel has been 
perceived. The arching theory suggested herein may account fot that phe- 
nomenon also, since interlocked sheet piling i is a fully supported column. Be 
In general, it has been difficult for the writer to be completely satisfied with 
? i any of the theories as to why cellular bulkheads function successfully nareegg 
 — the thought that arching of the fill actually forms a “bottom” 
“bottoms” to the “pail” or cell. _W hen a ‘materialj like? ; Plastic clay or is 
use used to fill the cells, the action becomes uncertain , variable, and completely 
different from the action caused by materials like rock, sand, or gravel. It may 


be that the incapacity of such plastic material to arch accounts for the difference. 


serving the sand dome in the broken cells as described herein, and as s yet pon 
~ found nothing to disprove his surmise. On the contrary, when such an actual 
_ bottom | is assumed, the guides | obtained | from the application of t the > basic: 
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_PENNOY ER ON CELLULAR COFFERDAMS 

oo Cellular Cofferdams 7 in , Clay-- —In his former paper (10), the writer did not 

— attempt to discuss cellular cofferdams founded on anything except rock, be- 
~ cause he considered there were sufficient intangibles involved in that problem. 
_ Even with Professor Terzaghi’ s contribution, it would seem that the proper 
be of cellular cofferdams in ‘clay (or even through overlying strata of clay 

or silt on top of stable strata), especially when submerged, is still largely de- 
on making mathematics coincide with what experience, observation, 
"stratum of unstable soil, is still a p aides largely dependent | upon on judgment 
_ based on precedent and experience. — Itis axiomatic th that, as soon as the sheet-_ 


pile structure — into soil, it becomes partly a cantilever structure. . Ifit 


= 


(Rock, Sand 
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Original Unstable soil 
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Assumed Bottom 


of G Struct 
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Counter Moment and Cantilever £ 
Portion of | Bulkhesd 


Bottom of 


Sheet Piling 

extends deep” enough and has insufficient width (we eight: of contents), it no 

longer i is a gravity structure but a , cantilever structure. " So far as the writer 

knows, no attempt he has been made, even theoretically, to design a , cellular sheet- 
- structure acting entirely as a cantilever. When all the intangible factors 
_are taken into consideration, the problem wl seem to be indeterminate. ‘i : 

‘Fig. 12 illustrates a design problem quite usual in practice, to the solution 

of w hich the paper. has rendered little of tangible value . In this A ea 
Addition to the notation of the paper, R is: the resultant of all applied forces, — 

vertical and horizontal; and X’ and ’ are moment arms as shown. | Successful ca 
"practice has s indicated that a can be considered as a structure 
down 1 to an assumed elevation, or “bottom.” The selection of this fictitious _ _ 
bottom i is a matter of judgment, all factors into consideration, including | 


the width of the cell. This elevation i is not selected the “coun- 


ssive re: ces 
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“moment the calculated width of the cell would indicate greater 


(heel) of the in mi thle case the structure would be assumed 
_ eapable of resisting the support stresses of a cantilever and there is no certainty 
“that the structure is capable of transmitting the internal stresses properly. 
That part of the structure below the assumed “bottom” is thus considered as a 
foundation only for the gravity structure above this elevation, transmitting the 
‘compression forces due to the overturning moment and weight of cells into the 
‘soil below. resistance as that of this this foundation— cantilever against 
overturning, tension at the heel, and shear resistance against lateral distortion— 
may be assumed to be an indeterminate, but added, factor of safety, 
_- The structure above t this elevation of the assumed “bottom” ’ is then analyzed 
for proper width, based on the width necessary to bring the intersection of the 
resultant of the forces represented by the Ww eight of wap contents of the cell and 


“of the cells must also be brought within : a proper Betiegs so as not to exceed the 
safe interlock tension in the sheet. is no a consideration 


method ap- 

‘ the proper width of cell is open serious sbjections’” see sen- 
tences following Eq. 8); but he has proposed no alternate solution. The prob- 
lem of how to determine the proper w vidth of a cell under the aforementioned 

- conditions, which are subject to variation on 1 every project, has been dismissed 
with the very doubtful statement that the : average width of such structures on 
"rock has been 0.85 H and w ith the further note that internal shear i is usually 

‘ the governing factor. — - Obviously, the shear resistance of the contents of a 
_ structure under the aforementioned c conditions could be much greater than the 
shear resistance ‘of the contents of as a structure founded o1 on rock, which does not 

> ae with | special reference to the internal w walls. Eq. | 32 giving the vertical 

; = height H for a given bearing value of clay is of f questionable practical 
value wad distribution of tlie ‘Pressure at the toe of the cofferdam, due 
ng, has been disregarded, the deeper 

= the sheet piling is driven, the greater the distribution will be. és Safe construc- 


tion, as it concerns s possible settlement at the toe, has always been assured when 
s the sheet piling was driven to high resistance against further penetration, 
: _ provided, of course, that there was no underlying unstable | stratum below the 
Bs After reaching the conclusion « expressed by Eq. 32, the author proceeds to 
analyze mathematically the internal stresses in cellular’ cofferdams on a com- 

pressible base, using the ‘ “customary” values of b = 0. 85 H, etc., which, i 
Part I, he stated were the av average ratios for sand-filled celis built upon va 
= J In practice, there i is no ‘relation between the two conditions; nor is em ratio of 
—6b = 0.85 H applicable, except by chance after the “middle ‘third” m hod de- 
- termines the proper width based on the assumed overturning moment rene the 


= 


e 


| 
— : 
4 
— > 
£ 
— 
B 
| 
| 
| 
orl | 


‘1946 PENNOYER ON CELLULAR COFFERDAMS 

™ 4 Also, the author has analyzed the division of the total overturning moment 
ty, on a cellular bulkhead (see Fig. 11) between the shear resistance of the fill and 
ie ‘the sheet piling. i He attributes the fact that the bulkhead began to deflect — 
ed a (in spite of theoretical analysis) to the probability that the real — 
ty lock tension (friction) might actually have been smaller than calculated. The 
i Oo . writer doubts whether the actual overturning moment could be determined as 
ae accurately : as outlined, since a few of the uncertain factors involved include the 
the assumption a as to lateral forces, the arm of the overturning moment or point of © 
the rotation, the fact that the structure I is partly a cantilever and partly a gravity © 
nst : - structure, the amount of cantilever support dependent o on the lean of the cells, 
a Mie the supporting value of the soil il (both | under and on the — of the ull 
Cellular Cofferdams on Deep Sand Strata.—Professor ‘Teranghi states 
the “Part II. Design of Cellular Cofferdams”) that a review of current methods of _ 
wed - designing cofferdams reveals an “amazing situation,” inferring that the seepage 


- of water, possible boils in the bottom, and loss of stability of the foundation | 
soil inside a cofferdam have not received consideration. _ The writer does not 
know of any phase in the design of cofferdams which has received more atten-— 
f tion from designers and contractors. In 1935, the writer presented a , discus- 


4 sion (29) « of this ee. This article offe ered some rule-of-thumb ratios between 


in single-wall when driven into various types of soil. These 
were based on long observation of successful practice in preventing bottoms 


from boiling seriously or ‘“‘blowing up.” The ratio given for water- bearing» 

~~ sand was a penetration of 30% to 40% of the unbalanced hydrostatic head, — 
ally 7 provided the cofferdam is properly pumped, as against the 67% mentioned i in 


‘the pe Also, , the author recommends a , penetration of 2 H into sand as being «an 
ad sary to ) prevent possible settlement at the toe | of the cofferdam, due te to vertical 
a ; load on the sheet piles. It is commonly known that sand is probably the s soil Se 4 
most resistant to 1 the penetration sheet piling and, under most conditions 
- encountered in practice, the writer doubts if sheet piling could be forced into 
the depth of sand, i indicated as necessary in the paper, jetting, without 
rupturing interlocks, ‘or without encountering absolute refusal. The writer has on iq 
test results before him giving safe bearing loads of sheet een hich show that 
Z. single pile tested to 47 tons w hen driven 25 | ft into sand and w when interlocked 


with two. other sheet piles tested to 120 tons. hen one considers that for 


tion, 

-the even a 60- -ft height of cofferdam the load on the sheet piles at the > toe could 

4 ; “hardly exe exceed 7 tons pe per foot width of sheet pile, 2 3 H, or about 40 ft, of penetra- 1 
Js to tion, seems at wide variance with the actual need. 


In the discussion of filters and gravel toes for taking care of pie Se | 
author has omitted to mention the usual methods of stabilizing the soil inside 


hn on the usual methods of stabilizin 

rock. eofferdams by pumping from sumps and the remarkable achievements of. the 


‘The discussion of breakwaters founded on sand includes a recommendation | 
that, the sand should be compacted to a depth of at least 1.5 H in some situa- ‘s 
? tions. - Under the conditions in Fig. 9(f), , Showing water going | over the top of * 

the breakwater, the critical consideration | for stability is the fact t 
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PENNOYER ON CELLULAR. -COFFERDAMS 
‘a turbulent w water scours the bottom on the harbor side until the brea akwater 

fails, wren protected by a stone berm. Several articles on breakwaters pub- 


> 


‘ 


to the design of adder on a more basis. 
. ; = He has as demonstrated | that few y problems have been 0 overlooked i in the past by the 
7 engineers an and contractors who have had experience with the practical considera- 


- tions necessary to design and build such structures properly and who have 


"i tried to reduce some e of the intangible factors to some semblance of a precise 


basis for analysis. _ The paper is largely confined to a mathematical analysis. 
— of a cellular r bulkhead ¢ of hypothetical : size and w reight of fill under a condition 
that is not common in n praetice. The : attempt to apply the analysis, howev er, 
7 serves to indicate a wide variance between theory and practice. — _ Dismissal of 
- the old “middle third” method | of approximating the proper | Ww idth « of a bulk- 
head under r the w ide variety of euay-enne met in practice by the substitution 
of the “average,” a “acceptable,” or ‘ ‘usual” Tatio of width equal to 0.85 H i is 
not only erroneous but dangerous. 


& In estimating the interlock tension, a factor of 0.4 for the com- in 
‘a putation of the lateral pre pressure ina ‘sand fill would be too high under s some | ta 
conditions and dangerous in others. No method is proposed whereby the s.. 
“critical interlock tension can be e determined or the 1 maximum 6 size “of cells fe: 
established to meet the wide variety of conditions encountered in practice. th 
The conclusion that consideration of the shear value of the soil in cellular me 
» cofferdams is more critical than the determination of the proper width and st 
maximum size of is fallacious, li 

(10) “Gravity Bulkheads and Cellular Cofferdams,”’ by Raymond P. Pennoyer, 
Lod Civil Engineering, Vol. A 1934, pp. 301- 305. (a) Eq. 18, p. 303. 
p. 305. (c) Eq. 12, p. 303. @ Fig. p. (e) Fig. 1, p. 301. 
(11) vU. Ss. Standard Steel Sheet Piling, Carnegie-Illinois Steel Corp. (Pitts- 4 
burgh, Pa., Chicago, Til. 1st Reprint 1936. 
(29) “Design of Steel Sheet-Piling Cofferdams ” by Raymond P PL Pennoyer a 
and George Hockensmith, Civil Engineering, Vol. 5, 1935, p. p. 19, 
“Steel Breakwaters,” by W. G. Caples, The Military Engineer, -March- 

en “Experience with Breakwaters at Milwaukee,” by D. A. D. Ogden, ibid.. A 
A (32) “Breakwaters on the Great Lakes,” by W. F. Heavey, ibid., November- x 


(33) 3) “Developments in Steel Breakwaters, WwW. G. Caples, ibid January- 

Corrections f for or Transactions: In September, 1 1944, "Proceedings, on page 


1016, change line 3 to read: ‘cylindrical cells of the type shown in a 
‘Fig. 1(a), by ares on the outer side of the ~The 
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PRESSURE ON THE LINING 


CIRCULAR ‘TUNNELS IN PLASTIC SOILS 
4 ‘Dis cussion | 


iscussion 


M. Bmp,’ Assoc. . Soc. E. paper is an excellent 
interpretation ¢ of field data collected, for the most “part, 
taken during and after the construction of the Lincoln Tunnel | under the 
Hudson River at New York City. In addition to analyzing these. se data, 
fessor Krynine has compared the pressures on the Lincoln Tunnel with | those on 
the Central Avenue Pipe Tunnel in Detroit. A few data referring to the 
Chicago Subway are also cited by the author. All three tunnels were con 
5 structed in plastic soil. This paper is a definite contribution to 10 the extremely 
a limited literature on earth p pressures on tunnels. The writer feels that 
i interpreting these f field data, Professor, Krynine has revealed a very important 
pattern that the earth  pocnneeen a on tunnels constructed in 1 plastic soils appear 
to follow. — A study of the fi field data for the Detroit Tunnel indicate S a pressure 
pattern ‘similar, in part, to that of the Lincoln Tunnel. a A study of the earth 
pressures” observed on other tunnels see seems to ‘substantiate ™ existence ofa 
pattern that the earth | pressures follow. 
7 in three parts. ‘ Condutan 1, 2, 4, 
“and 5 5 will be discussed i in Part i aie 3 in Part 2, and conclusions 6 onl 4 


_ This fact appears t to be scnanailaie it the field data are t taken at f face value. 
-The relief in pressure during construction and continuing for a certain short 


after construction has been | as brought out in conclusion 


7 __ Norg.—This paper by D. P. Krynine, M. . Am. Soc. C. E., was published in May, 1944, Proceedings. 
. Discussion on this paper has appeared in Proceedings, as follows: ptaagascige 1944, by M. A. Drucker. a 


xas, Austin, Tex. 
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PRESSURE IN PLASTIC SOILS Discussions 


semen there are two stages of pressure variation that must be considered in 

an analysis of data concerning pressure on the lining of a tunnel whether the 
tunnel be subaqueous or subterranean. Z The first stage is one in which j pressures 7 
exist due to the action ofa a disturbed soil mass tl that i is in a condition of f temporary — 
equilibrium. | During the second stage, ‘pressures exist due to the action of a 
soil mass which is undergoing readjustment in ‘Teaching a final c« condition of 
- equilibrium. The first stage takes place during construction and for a “od 
time following completion of construction. During this first stage, the vertical 

pressures on the lining are decreasing. Pressure- variation records for the 
Detroit Tunnel indicate that the vertical pressure on the lining decreased for a 
period of approximately three months. Similar for the — 


standard loading, decreased “systematically for the 253 days of 
“until it apparently a assumed a constant value” (see conclusion 4). Iti is evident 
that the equilibrium | of the earth mass was disturbed by driving the lining 
through the soil. — Also, the soil mass is in a tempors ary condition of equilibrium — 
during this stage since another variation in pressure, indicative of a readjust- 
_ ment of the soil mass, must occur before the final equilibrium is reached. — Since 
; long-time pressure observations were made at the Detroit T unnel (from De- 
_cember, 1930, to January, 1941), it is possible, by examining the pressure varia- 
7 tion during this period, to find evidence of the two stages. of pressure variation. 
a records indicate that a decrease i in | the vertical pressure on the lining pre- | 
vailed for a period of approximately three months (stage 1) } and thereafter the 
_ pressure was constantly i increasing for several years (approximately five years) 
4 until a final and permanent condition of equilibrium was reached (stage 2). 
ee Since, after this five-year period, the vertical pressure appeared to remain fairly — 
constant, it is not unlikely that a final condition of equilibrium was reached. 
_ The period of observation for the Lincoln Tunnel is short, and, ‘as the author 
cS states, no definite conclusions can be drawn from the information published to 
3 ee date. It is likely, however, that the period during which a steady decrease in a 
the vertical pressure was noted (253° days) corresponds to the first ‘Stage of 
= pressure variation and that a second stage during which the vertical p pressures — 
ao _ increase steadily will follow. d In the case of the Detroit Tunnel, the vertical 
pressure increased until it nearly equaled the static head at the elevation of the 
measured pressure. Probably, a similar. stage ‘occurred in the case of the 
Lincoln = unnel. iY T he author states, however, that a final stabilization of pres: 


- gure was reached at the end of the : stage of decreasing v« vertical p pressure. 7: This” 


seems unlikely tothe writer, 


| 
oOo 


7. 


Se: a i on As the author states in conclusions 1 and 2, the. decrease i in pressure is caused 
wee by vertical shearing stresses within the entth mass. . It must follow that an 
‘eel increase in pressure, as was noted on the lining of the Detroit Tunnel, is due to. 
diminishing or vanishing of these shearing stresses. During the of 

temporary equilibrium when the vertical pressures are decreasing, maxi- 
mum resistances in the mass are being g¢ mobilized to transfer v weight aw ay from 


Ske the tunnel,’ and, during the period of increasing vertical pressures culminating 


‘Earth Pressure on Tunnels,” by W. S. Housel, published as part of the Symposium on ‘‘E arth 
Pressure and Resistance of Plastic Clay,” Transactions, Am. E., Vol. 108 (1943), p. 1045. 
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December, 194, 4 ‘BIRD ON PRESSURE IN PLASTIC SOILS 


riter believes that soil mae, upon being disturbed, tends to to an arch 


"through the soil arch to the earth 1 mass ad Ijacent the tunnel. 7 This iden was 
presented by the late J. C. Meem, M. Am. Soc. C. E., in 1907. ie ‘The two ideas, — L 
one of arch action within the sc soil and the other the relief i in pressure through © ys 
vertical shearing stresses, are e merely two aspects of the same 1e phenomenon. 
_ For subaqueous tunnels, it may be considered impracticable to assume that a a 
definite relief in pressure will result from such arch action... However, 
material such as the organic clayey silt of the Hudson River and the blue « clay 
‘in n which both the Detroit and Chicago tunnels arc are built, a soil arch is probably 
formed. In the case of the Lincoln Tunnel, the creation of the bulge or the - 
te at the top boundary of the earth mass during the initial period of the: soil’s a 
disturbance appears to indicate that this arch action is initially in effect. —dAtis: a 
unlikely that the soil arch will endure for : a long period | of time; in other | words, 


down and back to the tunnel if not all, of the 
_ of the overburden (preconstruction pressure). _ The writer likens this breakdown 
“of the arch action in the soil mass to the action that would take place within “ * 
asphalt beam that was initially deflected upward and allowed to remain over its - 
= supports for an indefinite time. | Although the beam would maintain this" 
deflection for a certain period, eventually, due to its own weight and possibly — 
“due to vibrations from external causes, form of flow would cause the 
beam to return to its original horisontal position. 
The nature ¢ of the stresses set up to relieve the vertical pressure on the tunnel Bie 
is of great interest to the foundation engineer. Of much greater i impor-— 


tance, however, is the necessity of answering the following three questions: A 
How much decrease in pressure can be expected from vertical shearing 
: str resses, and for what period of time will this relief exist? ery a 
a 2. Why i is not this relief in weight pana, and during how long a period “4 
0 the shearing stresses diminish? 


3. Do the shearing stresses finally vanish or do they transmit, ar. ee 


“feta studies like ass for the Lincoln and Detroit tunnels and (2) an eam 3 
tion of the data from these field studies like that i in this paper. the apg . xe vy 


know fully how much decrease in pressure can n be expected, n more 


_proximations. The exact shape of these surfaces is and 
on the type and condition o of the soil mass in which t the maeene stresses occur. 


“The Bracing of Tae and Tunnels with Practical Formulas for Earth “Pressures,” ” by J. 
m, Am. Soc. C. -E, Vol. LX (June, 1908), p. 1. 
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“4 ‘effects of vibrations set up asa of externally : ‘Evena slight 
vibration might cause sufficient readjustment of the e soil eunthiies to cause a 
redistribution of stress. _ The period of time for the transfer of load back to the 

- tunnel is a variable dependent on local conditions. It may be a few months 

& but, as was the case for the Detroit Tunnel, it may be a number of years. That 

4 ‘some permanent or residual shearing stresses exist in the soil mass as a result of 

* 4 the dis disturbance ‘to its ; equilibrium i is indicated by the fact that the lateral pres- 

sures a against the tunnel lining of the Detroit Tunnel remained fairly constant 

and were ‘substantially less than the static head after final equilibrium was 

et reached. W. 8. Housel, M. Am. Soc. C. E., states* that an indication of a , 

permanent i change in in t the condition of equilibrium ir in the soil mass fr from its 
preconstruction state 1 may be found in the fact that the tunnel structure weighed 


= than the displaced earth and yet, in we equilibrium, “ore stresses of 


A similar was used the author i in determining the 


average static head. 4 ‘A comparison of the average ve vertical pressure and the 
"average static head should indicate that corrected average pressure sta-_ 


, is still considerably below the 
4 average static head value of 39. 2 Ib per 84 ig Probably some increase in the 
“average vertical pressure has occurred since this date. _ It is possible that the 
pressure will become approximately equal to the static head, as was yr case “7 


for the Detroit Tunnel, 5) wh ‘ 


stresses actin on =. mass.’ 
This decrease in pressure is indicative of the action of shearing stresses and only 
Br, by an increase in pressure will the gradual disappearance of these shearing 


stresses be indicated. 
eel Part 2 ‘2,—The auth author states under the heading, “External Forces Acting on 


the Lining and on the / Adjacent Silt Silt Mass: Interpretation of the Curve W, — 


ot od te “Tue to symmetry, there can be no horizontal shear and hence no verti- 
uy cal shear at point O (center line of lining). In other words, there can be no — 


relief in the weight of the overb t that 


_ plotted from a horizontal line, has a zero ordinate and a horizontal tangent — 


Fee. “The answer to the quéstion 
ree months ton year, approximately. | 

= lasted from three mo 
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December, 1944 BIRD ON PRESSURE IN PLASTIC SOILS 
“73 Unquestionably, aw to symmetry, there can be no vertical shear at sated 0, 
but, if there is to be no relief in the weight of the overburden at that point, the | 


te 


Ag from point O will be a ** Ax. The pressure at the base of the element is a 


pounds per square inch, which corresponds to the measured pa Shown in Table i 
“s riting the equation LF v = 0 for the forces 


on the element: 


p=W.—-Pi= 


a or in which the weight of the overburden expressed 
| in pounds equals h Az; and y = unit. weight 


To have no relief in weight at the center line, the rate of change of the 


‘shearing stress og must equal zero at point O. _ The auth author apparently makes _ 
this assumption in conclusion 3bys stating that the diagram of shearing stress 
has a horizontal tangent at the center line. writer admits that this may 
ue be the case. E Evidence of f such a a condition will come only through experimental = 


7 study. How ever, it does appear that the field data fail to support this. theory, - J 
: provided symmetry is assumed. A study of Fig. 7 shows that, on July 8, 1935, 7 
there is an overload at point O (Pa being greater than W.) and also that d down- 
ward vertical shearing stresses exist, although the values of these stresses have iB 
been omitted | from Fig. 7. The riter has di drawn the complete (W. — P4)-curve 
f for this date and has determined the downward vertical shearing stress at aed 
an n approximate fashion (see Fig. In making this rough computation 
of the vertical shearing 1g stress at point | O, the writer followed a method similar 
- to that used by the author in his determination of the shearing ‘stresses at the - 
other pc points on the periphery. The area under the (W. — P.)-curve between 
3 zero point (4.8 ft to the left of point O) and point O was computed on the e 


- 


= 


basis that the curve was parabolic. This area, represented _by 2/3 X 4. 
(1.6 144), equals 740 Ib—the accumulated relieved weight between the 


two points in question. The area of action of the ‘shearing stress is is 28 eq ft, 
if using the value of 28 ft for the height of the silt column above p point O in x 
| ig. The: approximate downward shearing stress equals ———— 38 28 ft r 26.4 Ib 
On July 8, 1935, the vertical downward pressure Pa 39. 


— 
7 mple, if an element of the soil mass of unit thickness is taken directly above _ _ 1 a 
| | | i 

d 

— 
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BIRD ON PRESSURE IN PLASTIC SOILS” 


ownward) was 


Bronce 26. 4 lb per sq “tt. ‘Fig. 7 shows that, on July 16, 1935, and on 

- August 8, 1935, the weight of te. overburden was relieved by | the action of 

1935, no relief in weight, and stresses, at point 0. 
Accepting the field data as exact, it must follow that on the first three dates 

symmetry was lacking and the rate of. change of the shearing stress —~ was not 
equal to zero at point 0. - if pressure. intensities had been taken at , corresponding | 
Scare on each side of the center line of the lining, it would have been —_ 

Bios determine if symmetry did exist. _ Reviewing Table 1, dow nward pressures, 
Py appear to have been taken only at points O, D, C, and B on the periphery; 
. - it is impossible to tell whether there was a symmetrical condition about 
the center line on the aforementioned four dates. The deformations arising 
the decrease i in the horizontal diameter of the lining that commenced as — 
. lining was being driven ‘Gaongn the silt m may have caused the symmetry to 


m 


q 


17. brews got vinh wl) + 


a 


Part 3 —The decrease in the horizontal commencing 

te after the erection of the lining i is an extremely interesting phenomenon. ‘ The 
undeniable implication a arising from, this occurrence is that the horizontal p pres- 
on the lining must exceed the vertical pressure. “Such a relation between 
the horizontal and vertical pressure is not commonly expected i in an undisturbed ' 
mass. has been to be 0.7 for soil of 
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ay 


the type of the Hudson River silt.!° © Itap appears logical that this phenomenon 
. ~ should occur only in soils that have a considerable ratio of the horizontal pres- 
_ sure to the vertical pressure in the mass. Since the field data for the Lincoln © 
_ Tunnel show that the ratio of horizontal pressure to vertical pressure is 
“than unity (0. 74 on August 29, 1935), an bes apices of ‘this phenomenon i is 
necessary writer feels that the au 


examination o of 9 9 showing the horizontal pressures on ‘the tunnel 

lining indicates to the writer that the general shape of the curve from point O 

to point to N) follows the straight- line law (triangular pressure distribu- 


tion) as used by Rankine i in his formula. would lead to the conclusion that 


- cause must be s be sought at the very beginning of the shoving—ie a probable 


lateral pressures point O to point Conclusion 7 is excellent. Itis 
that the Rankine formula should not be used for determining 
pressure below point B’ since the pressure on this section of the lining i is shown 
General Remarks —Great care e should be exercised i in interpreting field d data 
8 ‘such as those for the Lincoln Tunnel. . Fig. 7 shows that the maximum relief 
in weight, occurring at point B B on the periphery, ranges | from 2.7 lb per sq in. 
on. July 8 to 5.8 Ib per sq in. on n August 29. Also, at the: center line of the lining ae 
(point O) the “‘relief in weight” varies from an overload of 1.6 lb per sq in. on > 


July 8 toa load of O Ib per sq in. on August 29. In view of the fact that the ae 
accuracy of measurements made with the special silt gage and plugs was limited © ily 4) 

- to the range, 1 lb per sq in. to 2 lb per sq in., it is entirely possible that much of _ 

the difference in W, and Pz may have been by the inability of the gages 
to record true pressures . Certainly the overload at point OonJuly8couldbe . 
attributed to this cause. . Iti is within the realm of possibility that there was es ce | 
‘Telief i in pressure. . Although the writer acknowledges that this i is not likely, he : ona % 

merely wishes to state that a ‘slight i inaccuracy in the pressure-gage_ measure- 


; ments could easily account for : a difference i in pressure of 2. 7 Ib per sq in. In oa 


= cu ft as the density of the Hudson River water and a value of 104 Ib per ¢ cu 1 ft a as 
_ the weight of the saturated clayey silt. An error of as little as 4 lb per cu ft i in 
¢ estimating the density of the silt could produce a change it in pressure on the 
3 lining equal to some of the values considered in this p paper as relief in weight. fs 
+ To illustrate the effect of a variation in density of the silt on the lining pressure, — 
us a column of.the silt mass 20 ft high and 1: sq in. in cross section, with a volume | - = 

P 240 cu in. or 0. 14 cu ft, has been selected. OA variation in density of 4 lb per 


cu ft would cause a change i in pressure ( of 0.56 Ib | per sq in. at the base of thi 


3 column. - Such an error in detetmining the density of the silt mass is certainly 
possible, since there is doubt i in the writer’s mind that a homogeneous, uniform * 
ig material exists for 20 ft below the bed of the Hudson River. In addition, it 


. does not seem likely that the bed of the Hudson River is a level ground surface. — 


1932, p, 248. 


10 of in Unstable Ground,” by Robert G G. Hennes, Engineering, 
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BIRD PRESSURE IN PLASTIC SOILS Discussions 
‘Therefore, it is s difficult to —_— using - values for the height of the silt columns 4 
above t the tunnel lining expressed to the nearest 0.1 ft. — . 

‘ments s of this accuracy would be most difficult to make i in n the ne near fluid: silt. r 
A slight variation in the height of the overburden could cause substantial 
changes in the values of W, and could be responsible for some of the difference 
between Wa and P, hitherto attributed to relief in weight. Inaccuracies in 
gage readings, ¢ errors in estimating the density. of the silt mass, and errors ins 
determining the height of the silt column above the tunnel lining could produce 
results giving a pressure difference less than that w which a actually existed. . The 
only purpose in emphasizing these possible : sources : of error is to caution against 
too strict an interpretation of the field data. Unless these data are exact (a 

fact which any engineer familiar with field and | laboratory measurements will 

doubt), misle ading results may be obtained. criticism of the author's 

7 As valuable interpretation of the field data on the Lincoln Tunnel is intended. _ 

; a his excellent p paper is a 1 real contribution to the literature on earth p pressures 
on tunnels. — Many more similar studies are needed before definite conclusions | 


can be reached a as to the variations in pressure on on tunnel linings pera “) i 
_ Corrections for Transactions: In May, 1944, Proceedings, on page 6 632, change | 
line 21 to read ‘* * * order. = If two bodies in contact tend to move along. the 
plane of contact with respect to each o other, the force * * *”; on page 637, 
change line 22 to read a * * * lining a as in Fig. 1a), Pa and P, are forces 
computed f from ‘and, after after line 41, page 637, add the ‘following new 


Fig. 5 and i in all figures and discussion capital s desig- 


to its longitudinal axis. These forces are “represented 


and P, in ‘Fig. Unit p pressures on the lining, 
~ in vertical planes parallel to its s longitudinal axis are ordinates of the corre-- 


curves. The use of lower case letters for “unit pressures was thus 
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